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CHAPTER I 


INTRODUCTION 

Recycling of most of the consumables in manned spacecraft be- 
comes essential in extended missions. The most urgent task is to 
supply a continuous stream of breathable oxygen to the cabin to 
balance the consumption. This can be derived from the metabolically 
produced carbon dioxide. In recent years several systems have been 
developed for that purpose. They consist essentially of some device 
to remove the CO^ so as to keep its concentration in the cabin atmos- 
phere below a prescribed level, followed by a reactor to decompose the 
CO2 eventually into O2 which is recycled to the cabin. 

One particular method has been adopted by NASA as the most practi- 
cal and reliable method for missions lasting six months or more. The 
recovery of O2 from CO2 is achieved in three steps. Fir'st, CO2 is 
concentrated and mixed with H2 in an electrochemical cell (sometimes 
called a hydrogen depolarized cell or HDC); the mixture is then trans- 
ferred to a catalytic reactor where it is converted into water and 
CH^, the methane is then dumped into space; and finally the water vapor 
is electrolysed in a special electrolysis cell (called water vapor 
electrolysis cell or WVE) to O2 which is recycled to the cabin and H2 
which is transferred to the concentrator. 

Electrochemical CO2 concentration provides several advantages. 

The process is continuous and the equivalent weight of the system 
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needed to maintain the atmospheric level of CO2 concentration ( 0.25 
mn Hg) is considerably lower than the adsorption systems which do the 
same duty. CO2 compressors which are frequently needed in other sys- 
tems, e.g., in the vacuum desorption from molecular sieves, are 
totally eliminated. In addition, the concentration of CO2 protects 
the catalyst in the Sabatier reactor from contamination by the cabin 
air and supplies the CO2/H2 mixture in proper ratio to the reactor. 

The development of the device was undertaken in the past five 
years by two private contractors, Life Systems, Inc., and Hamilton 
Standard Division of United Technologies. The designs agree in 
principle but differ in detailed mechanical design and rperating 
conditions. One of the major problems encountered in designing the 
CO2 concentrator is the need to use a concentrated carbonate solution 
as a cell electrolyte. The air temperature and humidity determine 
the electrolyte concentration at steady-state. At a high air relative 
humidity the electrolyte becomes more dilute and increases in volume. 
This may lead to electrolyte run-off in the gas cavities and is avoided 
in the Hamilton Standard design by providing a reservoir to absorb the 
excess electrolyte at humid conditions. At a low air relative humidity 
the concentration of the electrolyte increases and may exceed the 
solubility limit. In that case precipitation would commence which may 
lead to the crossover of the H2 to the air stream causing a hazardous 
condition. This should be avoided by proper selection of the carbo- 
nate electrolyte and air humidity. K2CO2 was first used but was re- 
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placed by CS 2 CO 2 which has a higher solubility. Still, this did not 
allow for cell operation at air relative humidities below about 60%. 
Below this limit the cell performance decreased because of the loss 
of electrolyte and precipitation of CsHCO^ at the anode- It was 
necessary to use extensive humidity controllers at the air inlet which 
added to the total system weight and power requirements. Through its 
efforts to overcome this problem, Hamilton Standard introduced aqueous 
solutions of a new electrolyte, TMAC or tetramethyl ammonium carbonate, 
which have such low water vapor pressures that they can be used safely 
for air relative humidities as low as 35% without significant change in 
performance . 

Very recently, Hamilton Standard completed the design and fabrica- 
tion of a full-scale device for use in one-man air revitalization sys- 
tems and finished a 90-day test program to demonstrate the maintain- 
ability and durability of the device with a predicted life of two years 
(Huddleston and Aylward, 1975 a,b,c). The HOC and WVE cells are con- 
nected electrically in series so that the energy produced by the fuel 
cell reactions in the HOC cell is utilized in the water vapor electro- 
ly s is. 

The cost of air revitalization in a spacecraft is measured in 
millions of dollars and the final design and operation of the system 
should lead to the minimum cost for a certain duty. The economy of 
the system is based on two factors: minimum system equivalent weight 

and minimum power requirements. The weight appears in cost of pro- 
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pell ant needed to accelerate the spacecraft and it includes the weight 
of the HDC cells, WVE cells, Sabatier reactor and other auxiliary 
equipment. The power requirements are functions of the current 
densities used and the electrical efficiency of the cells. The total 
system cost depends greatly on the performance of the HDC cell. For 
example, a low CO2 transfer efficiency would result in additional 
cost of the power needed to make H2 and O2 from the H2O produced in 
the HDC cell. A poor H2/CO2 flow to the Sabatier reactor would lead 
to a low conversion efficiency which requires additional amounts of 
the expensive catalyst and would add to the weight of the system. The 
decrease in cell performance if precipitation occurs would require 
extra cells adding to the system weight. It is clear that an optimum 
design of the system is unattainable without a mathematical model that 
describes the performance of the HDC cell under all possible conditions. 

Lin and Winnick (1974) developed a fundamental model of the early 
versions of the device which used CS2CO2 electrolyte. The model was 
based on the steady-state transport equations and kinetics of the 
reactions which take place in the cell. Since then, several changes 
have taken place in the Hamilton Standard design; new electrodes and 
matrix materials were used and a new electrolyte was employed. It 
was necessary to develop a new model for the units which use TMAC 
electrolyte. The model should be able to simulate the new test data 
and predict the cell performance under all practical conditions. In 
that way a clear basis for design evaluation can be obtained and 



5 


optimization studies can be started. 

This work describes the development of this mathematical rrtodel . 
The starting point will be to review all information pertaining to 
the electrochemical CO2 concentration. This is followed by a brief 
description of the cell design, test data and their accuracy, and 
all useful information concerning the new electrolyte and modifica- 
tions of the cell design. Chapter IV summarizes the equations which 
form the basis of the steady-state mathematical model. The important 
model parameters are then evaluated and used to predict the cell per- 
formance over a wide range of operating conditions. Conclusions and 
recommendations of areas of possible improvements are finally pre- 
sented. 



CHAPTER II 


LITERATURE REVIEW 
Introduction 

In this chapter several topics are discussed which are inti- 
mately related to the electrochemical concentration of CO 2 . These 
include the chemical absorption of CO^ in buffer solutions, transport 
phenomena in electrolytic solutions, and theory of gas electrodes. 
This will be useful in the development of the mathematical model in a 
later chapter. 

Absorption of CO^ in Buffer Solutions 

The absorption of CO 2 in aqueous solutions can take place by two 
mechanisms depending on the pH of the solution. At a pH greater than 
10 the basic mechanism predominates 


CO 2 + OH" :> HCO" 

[ 2 . 1 ] 

HCO" + OH" ^ 003 ^" + H 2 O 

[ 2 . 2 ] 


The first step is slow and rate-determining while the second step is 
instantaneous (Kern, 1960). At a pH less than 8 the acidic mechanism 
predominates 


CO2 + H2O j H2CO3 


[2.3] 
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H2CO3 + OH" X HCO3 + H2O [2.4] 

Here the second step is instantaneous while the first step is rate- 
determining. At pH in the 8-10 range, both mechanisms are important. 
The rate equation for the first mechanism is 

d[C0j 

- — dt^= ki [CO2] [0H‘] [2.5] 

The rate constant k- was measured by Pinsent et al. (1956) over the 
range 0 to 40°c. The dependence of k^. on the ionic strength of the 
solution, y, could be represented by an equation of the form 

log k^. = g + u [2.6] 

where g is a constant dependent on the temperature and 5- is another 
constant dependent on the temperature and the ionic species in the 
sol ution. 

The rate equation for the second mechanism is 

d[C0,] 

at ^ ^iii 

The measurements of k.^.^. taken by several workers were summarized by 
Edsall (1968). The ionization constants of carbonic acid as obtained 
from cell measurements were reported also by Harned and Owen (1958). 

The problem of simultaneous diffusion and chemical reaction was 
treated rather extensively by Astarita (1967). He classified pro- 
cesses of chemical absorption into different regimes according to the 
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relative magnitudes of the characteristic diffusion and reaction 
times. For gas absorption with a first-order chemical reaction, the 

process is considered in the fast reaction regime if 

. . 

D|_/h2 » 1/k [2.8] 

v4iere is the molecular diffusivity of the absorbing gas in the 
liquid phase, hj^ is its mass transfer coefficient in the liquid 
phase, and k is the reaction rate constant. Under this condition 
the absorption rate per unit area of the interface, N, can be calcu- 
lated by the following equation derived from penetration theory: 

N = (C’^- C") [2.9] 

where C and C*" are the concentrations of the reactant at the gas- 
liquid interface (the gas solubility as ordinarily defined) and the 
equilibrium concentration, respectively. 

The physical solubility of a gas in a liquid containing electro- 
lytes which react with the gas cannot be determined from equilibrium 
measurements used for other gases. Only kinetic measurements can give 
an estimate of the solubility in that case (Nysing and Kramers, 1957). 
However, the solubility can be obtained approximately by the equation 
(Van Krevelen and Hoftijzer, 1948) 

log (cV°) = - c y [2,10] 

• 4 * 4 ’® 

where and are the solubilities in the electrolytic solution 
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and in water, respecti vely , and c is a coefficient which can be con- 
sidered as the sum of elements characteristic of the cations and anions 
present in the solution and of the dissolved gas 

c = + c_ + Cg [2.11] 

Transport Phenomena in Electrolytic Solutions 

Mass transfer in electrolytic solutions can be expressed by the 
Nernst-Planck equation (Newman, 1973) which, in one dimension, has 
the form 

dC. 

where N. is the flux of the species j in the x direction, D. is the 

J J 

ionic diffusion coefficient, C. is the-ionic concentration, z. is the 

J J 

ionic charge, F is the Faraday's constant, u. is the ionic mobility, $ 

J 

is the electrical potential, and v is the bulk velocity of the solu- 
ti on . 

Ionic mobilities at infinite dilution can be calculated from the 

O 

limiting values of ionic equivalent conductance, X., by the equation 

J 

u° = X° / |z^| [2.13] 

The ionic diffusion coefficients at infinite dilution can also be cal- 

o 

culated from the values of X. by the Nernst-Einstein equation 

J 

O O O 

D. = RT A./ Iz.l 


[2.14] 
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where R is the gas constant and T 1s the absolute temperature. This 
equation leads to the following equation for the diffusion coefficient 

o 

of a binary electrolyte at infinite dilution, D: 

D° = D° D° (z^ - z.)/(z+ D° - z_Dj [2.15] 

Several experimental methods have been described for the mea- 
surement of the conductances and diffusion coefficients of electro- 
lytes in concentrated solutions (Robinson and Stokes, 1959). These 
can be combined with the transference numbers, which can be obtained 
by cell measurements, to estimate the ionic mobilities and diffusion 
coefficients. In the absence of direct experimental data the dif- 
fusion coefficients can be calculated from the mean activity coef- 
ficients of the electrolyte, by the following equation 

o o d In 

o 

where m is the molality and n and n are the viscosities of the solu- 
tion and the solvent respectively. This equation is a simplified 
form of equations [11-67] of Robinson and Stokes (1959). For more 
accurate calculations additional terms are required to account for the 
electrophoretic and hydration effects. These terms change the values 
by a few percent and certainly are insignificant in the present calcu- 
lations since the diffusion coefficients are reduced drastically be- 
cause of the presence of the matrix. 
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The relation between the conductance and viscosity in concentra- 
ted solutions of strong electrolytes is described by the equation 
(Wishaw and Stokes, 1954) 

where A is the equivalent conductance, X is the electric field, and p 
and q are arbitrary constants. 

In the absence of transference numbers the individual mobilities 
and diffusion coefficients in concentrated solutions can be calcula- 
ted in a simple way (Lin and Winnick,1974) . The ionic diffusion 
coefficients can be calculated with the assumption that the relation 
between D^, D”, and D is the same as at infinite dilution 

D = D^D_(z^ - z_)/(z^D^ - z_D_) [2.18] 

and 

D+ = D° D_/D° [2.19] 

In the same way the ionic mobilities can be calculated from the equiv- 
alent conductance. The ionic mobilities and diffusion coefficients 
can then be calculated in a multicomponent solution utilizing their 
dependence on the viscosity of the solution. The values of U| and 
in multicomponent solutions are calculated by the relations 

Dn, = Dn/n„ [2.20] 
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and 

u„ = u n/Ti„ [2.21] 

m m 

where the subscript m refers to the multicomponent mixture. 

To calculate the effective mobilities and diffusion coefficients 
in a porous matrix it is necessary to multiply their values in ma- 
trix-free solutions by a matrix labyrinth factor, a. The most widely 
used equations for calculating a in a porous medium are the following 
equations by Wheeler (1951) and Mackie and Meares (1955), respecti vely : 


a = e/2 

[2.22] 

a = e^/(2-e) 

[2.23] 


where e is the porosity. 

Electroneutrality is preserved in the solution except in the 
double layer near the electrodes which has a thickness of less than 
100 A. Because there are several anions (CO^I HCO^ and OH ) in the 
solution and only one cation in appreciable concentrations the concen- 
tration of the cation varies across the matrix to preserve solution 
electroneutrality. This gives rise to concentration overpotential 
in the solution which adds to the power losses of the cell. 

Theory of Porous Gas-Diffusion Electrodes 

Electrode reactions involving one or more gaseous components re- 
quire an intimate contact between the gas and electrolyte in the vi- 
cinity of the solid electrode. Because of the capillary forces, this 
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can best be realized in a porous electrode. The performance of gas- 
diffusion electrodes depends largely on the penetration of the gas 
and electrolyte into the electrode pores. 

Several models have been described which simulate the microstruc- 
ture of the electrode pores and the contact betv;een the gas and elec- 
trolyte phases (Bockris and Srinivasan, 1969), but their success de- 
pends on the manner of preparation of the electrode. For instance, 
teflon added to the catalyst particles in the electrode matrix, such 
as those used in the HOC cells, decreases the wetting power of the 
electrode. The situation in this case is probably similar to that 
shown in Fig. 1 (Austin, 1968). The gas-electrolyte area decreases 
as the electrolyte penetrates into the electrode pores towards the 
gas side. For a certain gas pressure there is a limit for the elec- 
trolyte penetration. The capillary equilibrium in the electrode 
pores depends on the surface tension of the electrolyte and the pore 
radius (Markin et al . , 1966). 

The surface tension of the electrolyte depends on the concentra- 
tions and specific nature of the ions (Randles, 1963). Highly hydra- 
ted ions and tetraalkyl ammonium ions are classified as surface-active 
ions because they tend to increase the surface tension with concentra- 
tion at a higher rate than other ions. Electroosmotic forces also in- 
fluence the solvent transport in porous media under the applied elec- 
trical field (Lakshminarayanaiah, 1969), and result in the dependence 
of the electrolyte circulation on the current density distribution in 
the electrode structure (Huddleston and Aylward, 1973). 
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Fig. 1 Illustration of Partially Wetted Electrode 
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Mass transport processes play an important role in the performance 
of gas-diffusion electrodes. Convection and diffusion bring the reac- 
tants from the bulk gas phase to the electrode surface. Inside the 
electrode pores two mechanisms are possible for the transport in the 
gas phase; gas diffusion under a concentration gradient and gas per- 
meation under a total pressure gradient. 

Four types of gas diffusion mechanisms are also possible depend- 
ing on the relative magnitudes of the mean free path of the gas mole- 
cules and the pore diameter (Austin, 1968). In the electrolyte phase 
filling of the electrode pores; diffusion, ionic migration, and convec- 
tion may all restrict. Convective currents result from temperature 
and surface tension gradients in the electrode pores and random fluc- 
tuations of the gas pressure. This results in a relatively uniform 
concentration in the electrolyte throughout the electrode. 

Electrode processes are usually studied in relation to the over- 
potential-current density dependence. Electrode overpotential con- 
sists of three components; activation, concentration, and ohmic. The 
equations derived for the overpotential -current density dependence on 
the basis of simple mathematical models (Bockris and Srinivasan, 1969) 
do not predict the correct values because of the geometric complexity 
of the practical pore systems. Ho\f/ever, they can give insight to the 
desirable properties of the electrodes to obtain the maximum power 
output from a fuel cell. 
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PROCESS DESCRIPTION AND CELL DESIGN 
Process Description 


Carbon dioxide is continuously removed from the cabin air as it 
passes over the cathode of the HDC cell. Oxygen in the air diffuses 
first into the liquid film within the cathode pores and reacts to 
form hydroxyl ions 

O2 + 2H2O + 4 e' 40 H' 

[ 3 . 1 ] 

CO2 diffuses simultaneously and reacts with the OH’ 
carbonate ions 

i ons to form 

CO2 + OH' ^ HCO" 

[ 3 . 2 ] 

HCO" + OH' X + ^2^ 

[ 3 . 3 ] 


Pure hydrogen is fed into the anode side. After diffusing into the 
liquid film at the anode, it forms hydronium ions 

2H2 + 4H2O ^ 4H30'^ + 4 e' [ 3 . 4 ] 

The lowered pH of the anolyte favors a higher bicarbonate concen- 
tration according to the equilibrium [ 2 . 2 ]. These HCO3 ions dissoci- 
ate to CO2 which is removed by the H2 stream 

HCO3 ^ CO2 + OH' [ 3 . 5 ] 
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Because of the low pH at the anode, the HCO^ ions react to form car^- 
bonic acid which also dissociates to CO2 

HCO” + H^O ^ H2CO3 + OH" [3.6] 

H2CO3 ^ CO2 + H2O [3.7] 

These reactions combined result in H2O production, pumping the 
CO2 from a low partial pressure at the cathode side to a high partial 
pressure at the anode side, and the generation of electrical energy. 
The water vapor is carried by the air stream to the WVE cell and the 
CO3/H2 mixture is fed to the Sabatier reactor. The generated elec- 
trical energy is used in the WVE cell. The concentrator cell is 
cooled solely by the flowing air stream. 

Major Requirements 

The main function of the concentrator is to maintain the CO2 par- 
tial pressure in the cabin atmosphere at about 3mm Hg. However, it 
should be flexible enough to accommodate normal atmospheric concentra- 
tions of CO2 (0.25 mm Hg). For one-man requirements, a CO2 removal of 
2.20 Ib/day is needed. With the present design, it is possible to ob- 
tain a removal rate of 2.40 Ib/day at 3mm Hg (Huddleston and Aylward, 
1975c). The concentrator consisted of four identical cells. 

Energy requirements are also important in long-space missions. 

The design should allow no more than 50 yV/hr of long term voltage de- 
cay rate. The voltage decay rate for the combined WVE and HOC cells 
in the present design does not exceed 42 yV/hr. In addition, the con- 
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centrator should show stable performance in the air relative humidity 
range of 35 to 90%. 


Cell Design 

Cell configuration 

The concentrator is shown schematically in Fig. 2. It is based 

on the "cell pair" concept in which back to back hydrogen electrodes 

facilitate the maintainability of the system. The electrodes are 
2 

1 ft each and this relatively large size is chosen to minimize the 
total cell weight by reducing the percentage of the wasteful peripheral 
material. The cathode chamber is made of titanium and it has 96 rec- 
tangular air ducts 0.2 " x 0.19" each. The length of the air path is 
6 inches. 

Hydrogen flows in a serpentine path in a cross flow manner with 
respect to the air stream. The center housing is designed in such a 
way as to minimize the pressure drop in the stream (Huddleston and 
Aylv/ard, 1973). Hydrogen flows at a rate of 600 cm /min with a back 
pressure of 1 psig. The use of pre-bent housing employing spacers of 
rigid material guarantees a uniform matrix thickness at the specified 
Hg pressure (Russell, 1973). 

The cell is provided with a reservoir of porous 'Tissuquartz' 
which absorbs the excess electrolyte during cell flooding at humid con- 
ditions and returns it to the matrix during drying conditions. 
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Fig. 2 Cross Section in Cell Pair 
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Cell electrodes 

The design of the electrodes is of major importance for cell per- 
formance. Both the electrodes in the present design are made of 
platinum precipitated on a metallic screen with teflon binder. The 
teflon binder, being hydrophobic, minimizes flooding problems. The 
amount of catalyst loading and the manufacturing procedure determines 
the electrodes performance. The thickness of the electrodes is of the 
order of 3 to 7 mils. 

Cell matrix 

The matrix material and compression are selected to obtain mini- 
mum internal cell resistance, chemical stability under the operating 
conditions, and maximum CO^ removal efficiency. In the present design, 
the matrix consists of one layer of Tissuquartz (18 mils) sandwiched 
between two layers of fuel cell asbestos (20 mils each) and the whole 
compressed to 25 mils. This matrix set was tested while saturated 
with CS 2 CO 2 solution and a six-month endurance test showed that the 
matrix maintains sufficient compression resistance (60 psig) to pro- 
vide good interface contact between the matrix, electrodes, and 
housings (Huddleston and Aylward, 1973). 

Humidity control 

Previous versions of the device used CS 2 CO 2 electrolyte (Winnick 
et al., 1974). Cell drying occurred when the air relative humidity 
dropped below 60^. At these conditions, CsHCO^ starts to precipi- 
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tate at the anode because of its lov/ solubility. This required ex- 
tensive humidity control and cell cooling which was acceptable in the 
six-man subsystem. 

For a one-man subsystem it is desirable to use an electrolyte 
solution which has lower vapor pressures and thus has higher stabillly 
in low humidity environments. Investigation led to tetramethyl ammoni urn 
carbonate (TfiAC) which enabled the operation at air relative humidities 
in the range of 35-90% (Huddleston and Aylward, 1972). The solubility 
of NMe^HCO^ is about 88 wt% at the cell temperature while the solubili- 
ty of (NMe^)2C02 is 52 wt%. Therefore, NMe^HCO^ precipitation is not 
a problem as it is with CsHCO^. 

The NMe^ ions are somewhat unstable, forming tri methyl amine and 
methyl ion (Huddleston and Aylward, 1972). The latter reacts with 0H~ 
ions to form methyl alcohol. In the HOC cell, however, proper design 
and operation can prevent this decomposition. At the cathode, plati- 
num black is a poor decomposition catalyst in the oxygen atmosphere. 

At the anode, the pH is normally low and does not promote decomposition 
(Huddleston and Aylward, 1972). 

COo removal efficiency 

c. 

According to the overall cell reaction, as described above, two 
moles of CO^ can be transferred for one mole of O 2 and two moles of H 2 
consumed. Since the current produced is proportional to the amounts of 
H 2 and O 2 reacting this corresponds to 100% current efficiency, i.e., 
when one mole of CO 2 is transferred for each two Faradays of electricity. 
Sometimes this efficiency is expressed as the "Transfer Index" which is 
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the equivalent weight ratio of CO^ removed for each pound of O 2 con- 
sumed. A current efficiency of 100% is equivalent to a transfer index 
of 2.75. 

At very high partial pressures of CO^ in the air stream the concen- 
tration of OH" ions in the catholyte drops significantly and the HCO^ 
ions concentration increases. This may lead to the transport of CO 2 as 
HCO^ instead of CO^ and this corresponds to 200% current efficiency 
since now one mole of CO^ is transferred for each one Faraday of elec- 
tricity. In the range of operation expected of the HOC cell the concen- 
tration of CO 2 in the air stream will not exceed 15 mm Hg. The CO^ is 
removed through CO^ ion transport which yields less than 100% current 
efficiency. However, it is possible to obtain current efficiencies 
slightly higher than 100% if the CO 2 concentration exceeds normal levels. 



CHAPTER IV 


MATHEMATICAL MODELING 
Introduction 

Two approaches have been used to predict the performance of the 
concentrator under various conditions. In the first approach, test 
data at selected conditions are used to obtain empirical equations 
relating the performance data, e.g., the CO 2 removal rate to the im- 
portant parameters (Marshall et al., 1974). The second approach is 
based on a fundamental model of the reactions and transport mechanisms 
which take place in the cell (Lin and Winnick, 1974). This approach 
is more useful since the data can be extrapolated with more confidence 
to off-design conditions and can predict the effect of unstudied para- 
meters. Furthermore, some performance characteristics can be predicted 
which are either difficult or impossible to measure. One example is 
the variation of the concentration of the different species across 
the matrix; which can be used to predict the incidence of precipitation 
or, in other words, the limits of cell operation. These concentrations 
can also be useful in qualitatively predicting some of the important 
performance parameters, e.g., the anode overvoltage. 

In this chapter the various reaction mechanisms and prevailing 
equations for the steady-state case are described. These equations are 
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then combined to yield the differential equations which are used in 
design. The numerical methods used in solving these equations are 
also discussed. 

Mass Transfer and Reaction Rates 

The concentrator is considered to be divided into subcells in the 
direction of the air flow and each subcell to be subdivided into five 
zones as shown in Fig. 3. In each zone different reactions and mass 
transfer processes take place. These are examined separately for each 
zone. For convenience, the different species will be represented by 
numbers; 1 for HCO^", 2 for OH", 3 for 002^", 4 for NMe^^, 5 for CO2, 
and 6 for H2O. 

1 . Ai r zone 

In this zone CO2 is transferred from the gas phase to the electrode 
surface and then through the cathode pores to the liquid film inside. 

The flux of CO2 in this zone, Ng, is given by 

"5 = [4.1] 

where h^^ is the overall mass transfer coefficient, and are 
the partial pressures of CO2 in the bulk gas phase and at the gas- 
liquid interface, respectively, R is the gas constant, and T is the 
temperature. The subscript c refers to the cathode region. 
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Fig. 3 Schematic of a Subcell 
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J. 4. 

The relation between Pg^ and the concentration of Cg^ in the 
liquid film is 


*^5c ^ ^5c 


[4.2] 


where H is the Henry's law constant which is a function of the tempera- 
ture and the ionic strength of the solution. The average mass transfer 
driving force can be taken for any subcell as the logarithmic mean of 
the inlet and outlet values 


(VPc)l. 


m. 


(Ai“Ao)/ln(Ai/ A q) 


[4.3] 


4. 

where A is ^nd the subscripts I and 0 refer to the inlet and 

outlet of the subcell. 

The mass balance in the air zone is given by 


[4.4] 

where is the volumetric air flow rate and A is the geometric electrode 
area in the subcell. These equations can be combined together to give a 

4* 

single expression for ''5c 

° ■ '*6 ART/[V^(l-exp{-h5^ A/V^})]]/H [4.5] 

In the same way, the water vapor pressure of the electrolyte is given 
by 

Pec = <‘’6c>I - ”6 ''RT/{Vjl-exp(-h5^ A/V^)]} 


[4.6] 
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2. Cathode reaction zone 

Because of the high pH at the cathode, CO 2 absorption proceeds 
via the basic mechanism described by equation [2.1] and [2.2], The 
combination between CO 2 and OH" ions is the rate-determining step. 

The local rate of reaction, r^, is given by the equation 

r, = k^(CjC+ - Ci/Ki) [4.7] 

where k^- is the rate constant of reaction [ 2 . 1 ] and K.j is its equili- 
brium constant. The second step is instantaneous and the correspond- 
ing equilibrium is established. The equilibrium relation is 

03 /( 0 ^ 02 ) = K.. [4.8] 

The pH at the cathode is usually in the range of 13 to 14 

(Huddleston and Aylward, 1973). The rate constant is in the range of 

8500 sec"^ (mole/il)"^ (Kern, 1960). The characteristic reaction time 

-4 -3 

is thus in the range from 1.2x10 to 1.2x10 sec. The characteristic 
diffusion time is dependent on the electrolyte hydrodynamics which can- 
not be measured within the pores. However, in packed towers the dif- 

-3 -2 

fusion time is usually in the range from 5x10 to 4x10 sec 
(Astarita, 1967). Thus, the condition [2.8] is satisfied and the 
reaction can be considered in the fast-reaction regime. Mathematically, 
this is expressed by the equation 

*^*5 ^ ^c^*^i^5^2^^^^5c " 


[4.9] 
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where S is the total gas-electrolyte interface area per unit of elec- 
trode geometric area and Dg is the diffusivity of COg in the electro- 
lyte. In these equations the equilibrium constants are calculated 
from values at infinite dilution by multiplication by the appropriate 
activity coefficients evaluated at average concentrations. 

3. Bulk electrolyte zone 

2 - - - 

The transport of CO^ , HCO^ , and OH ions in this zone can be 
expressed by equation [2.12]. It is conceivable that all the water 
produced by the fuel cell reaction is removed by the air stream be- 
cause of the high air flow rate compared with the hydrogen flow rate. 
This was actually proved by measuring the dew point of the leaving air 
stream (Lin et al . , 1974). Accordingly, the water flux in the cell 
is given by 

Ng = - 1/2 F [4.10] 

where I is the cell current density. The average molar bulk velocity 
of the solution is given by 

v=^EZ. N. [4.11] 

J J 

where Zj is the mole fraction of the species j in solution, Cj is the 
total number of moles per unit volume, and the summation extends 
over water and all the species in solution. Since the individual 
fluxes are unknown beforehand, an iteration procedure is required. 
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However, water represents about 907o of the total number of moles in 
solution. For simplicity, it is assumed that the other species do 
not contribute to the net bulk velocity of the solution. In that 
case the net bulk velocity is given by 


V = - 1/2 


[4.12] 


In fact, it is found that the assumption of zero bulk velocity gives 
the same results as obtained with the above assumption which means 
that the inclusion of the other species in calculating the net bulk 
velocity is totally unnecessary. 

The individual flux of each component varies with the distance x, 
but the following sums are independent of x at steady state: 


I] + N3 + N5 = 


[4.13] 


1^ + N 2 + 2 N 3 = I/F 


[4.14] 


Here represents the total CO 2 flux. Since the cation does not 
participate in the cell reaction, its flux vanishes at every point 


N4 = 0 


[4.15] 


Substituting for the flux of each component from equation [2.12], 
equations [4.13] to [3.12] give the following differential equations: 


dCc 

C^) 


[ 4 . 16 ] 
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d$ 


dCi . dCp 

^1 W " ^ ^2 d)T “ ^^2^2 dx 


dC. 


[4.17] 


^ I(C,+C2^2C3)/(2FC^) 


dC. 

°4 dT " ^“4^4 f " • V(2 FCt> = ° 

In this zone reactions [2.1] and [2.2] are assumed to be in 
equilibrium. The concentrations C2 and Cg are given by the following 
equations: 

C2 = Cg/CC^K..) [4.19] 

C5 = C^/(C2K.) [4.20] 

Electrical neutrality also results in the following equation: 

^4 " C^+C2+2C3 [4.21] 

4. Anode reaction zone 

In this zone CO2 is released by the decomposition of the HCO3 
ions and H2CO3 by the reverse of reactions [2.1] and [2.3], res- 
pectively. The combined reaction rate constant can be obtained as 
follows: 

'"a ^ *^i^l “ '^i^5^2 •^iii[^2^^3^ " ’^iii^5®6 [4-22] 

where a^ is the water activity, and are the forward and back- 
ward rate constants of reaction [2.1], and and k|- 


are the 
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forward and backward rate constants of reaction [2.3]. The relations 
between these rate constants are 


ki = k./K. 

•^iii " *^iii/^-ii 


[4.23] 

[4.24] 


Since the reaction [2.4] is in equilibrium, [H 2 CO 2 ] is given by 


[H 2 CO 3 ] = C^ag/C 2 K.^ [4.25] 

where a^ is the activity of water and is the equilibrium constant 
of reaction [2.4]. From equations [4.24] and [4.25] we obtain 




[4.26] 


The equilibrium constants and are related to K^. by the re- 
lation 


= S' 


[4.27] 


By substitution of equations [4.23], [4.26], and [4.27] into equation 
[4.22] we obtain the following equation for the local rate of CO 2 
production: 


‘"a ^ (k.j+k^.^.|ag/C2)(C.,/K^.-C2C3’'') 


[4.28] 


Assuming the reaction to take place in the fast-reaction regime 


^5 ” y^5^'^iV'^iii®6^^^^^1^S'^2“^5a^ 




[4.29] 
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5. Hydrogen zone 

The released CO 2 is removed by the hydrogen stream. The rate of 
removal is given by 

"5 - '’BaC’L - PBa’/RT C4.30] 


where hg^ is an overall mass transfer coefficient in the anode zone, 
The equilibrium relation at the gas-liquid interface is 


pL = H cl 


[4.31] 


The mass balance gives the following equation; 


^Ba = [(Psa’l - N5ART/[V^(l-exp{-h5//V^})]]/H 


[4.32] 


where is the hydrogen volumetric flow rate in the anode zone. 


Mathematical Model 


The model describing the performance of the cell can be con- 
structed from the above-mentioned equations. For simplicity, the 
following dimensionless variables are introduced: 

j=l,2,...,5 [4.33] 

j=l,2,...,5 

j=l,2 5 

j=l,2,3,4 


N. = N.d/(C*D*) 

J J 

'd. = D./D* 

c. = c./c* 

vJ J 


u. = U.RT/D* 

J J 


I = Id/(FC*D*) 
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(?5)i = (P5 )i/(HC*) 

K = KC* 

S^d^C*/D* 
kg = k. SgdV/D* 

G - ag k.../(k.C*) 

X = x/d 
i = $F/RT 

f'Sc = 

i'Ba = h 5 a“/(RTD*) 

= V^Hd/{ARTD*) 

V, = V,Hd/(ARTD*) 

a a 

where C* and D* are characteristic constant values of concentration and 
diffusivity and d is the thickness of the cell matrix. 

By substitution and rearrangement in equations [4.7] to [4.32], 
the following two nonlinear differential equations are obtained: 

~ ~ dC-, ~ dC., 

fn (Cl ,CJ — + fp(C, ,C.) 

+ KC1+C3+C5-C1U1/U4 - 2C3U3/u^)/(2Ci-) 


[ 4 . 34 ] 
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^ ~ dC, ^ ^ dCo 
fo(C,,Cj -^+ f4(C,,C.) -^ = - I 
‘ dx ^ ‘ dx 

+ I(Ci+C 2+2C2 - C-jU-j/u^ - CpU2/u^ - 

where 

f-] = Di+D^(C-j “C^/K^- ^ } (U1+2U2C2/C1 )/ (u^ 4 *) 

+ 2D5K..C^/(K.C2^ 

f2 = D4(2c|+l/K..)(u^C^+2u3C3)/(u4i|j) 

+ D3-S5K,,E,2/(',.c32) 

^3 " D^+D^(C^^-C3/K^.)[u^+U2C3/(K^^.C^ ) + 4U3C3C1 


= D^(2C,+1/K..)(u^C^+U 2C3/(K..C^) + 4u3C3]/(u^ 


+ 203 + D2/(K,.^C,) 


and 41 is given by 


i|j = + 2C^C3 + C3/K. . 


[4.35] 
/u4)/(2C.p) 

[4.36] 

[4.37] 

]/(u^ij^) 

[4.38] 

[4.39] 


[4.40] 
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The boundary conditions are 

- C^/C 2 K.} at X = -if [4.41] 

*^5T "" ‘ ^^5ah 

- ^^ 57 /CVg(l-exp{-hgg/V^})]} at x = ig [4.42] 

As it has been assumed, all the water produced by the cell re- 
action goes to the air stream. In dimensionless form, equation [4.6] 
becomes 

~^lc = (Pec^I 

Since the catholyte consists mainly of (NMe^)2C02, its vapor pressure 
data can be used to obtain the catholyte concentration 

C 4 = F{P+^, T) [4.44]m 

The catholyte temperature for the Hamilton Standard design varies 
linearly along the length of the air path and can be calculated from 
the rate of heat generation and the cell dimensions (Lin et al . , 1974). 
The rate of heat generation, Q, is calculated by the equation 


Q = (AHpNg - EI/J)A 


[4.45] 



36 


where AH|^ is the heat of formation of water vapor, E is the cell 
voltage, and J is the heat/energy conversion factor. The last three 
equations provide the additional boundary condition 


+ C 3 /(C^K 2 ) + 2 C 3 = at X = -Js 


[4.46] 


Numerical Method of Solution 


The nonlinear differential equations [4.34] and [4.35] describe 
the variation of and C 3 across the cell matrix. The concentrations 
of other components, C^.C^, and Cg are implicitly included in these 
equations and can be calculated from the equilibrium relations [4.19] 
and [4.20] and the electroneutrality relation [4.21]. The complexity 
of these equations, apart from the high non-linearity of the co- 
efficients, lies in the fact that these equations include the unknown 
parameter, Ngy, which in turn is a function of the boundary conditions, 
These boundary conditions consist of two relations between C-j and C 3 , 
[4.41] and [4.42], which implicitly include the parameter Ngy. An in- 
dependent boundary condition, [4.44], is obtained from the water and 
energy balance on the cell and it relates C,j and C 3 at x = ~h. 

Shooting methods (Roberts and Shipman, 1971) are suitable for 
this kind of two-point boundary value problems. The procedure is to 
assume a value for to calculate C.| and C 3 at x = for this 
value of %TV and then to integrate the differential equation from 
X = -Jg to X = %. The boundary condition at x = can then be calcu- 
lated and compared with the guessed value. Some of the shooting 
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methods use the difference between the calculated and predicted values 
to make corrections for a new guess. However, this cannot be used 
here because of the dependence of the differential equation and 
boundary conditions on the same parameter, Ngj. Instead, the new guess 
is chosen higher or lower according to whether the calculated value at 
X = ^ is higher or lower than the guessed value at x = -h. 

Successful application of shooting methods greatly depends on the 
stability of the numerical integration technique used. Predictor-cor- 
rector methods are of interest because of their relative economy in 
terms of the number of derivative evaluations required for a given 
order algorithm. The method used here is based on a predictor-corrector 
algorithm developed by Crane and Klopfenstein (1965). The algorithm is 
fourth-order and has higher absolute and relative stabilities than 
other predictor-corrector methods. During integration^ the step size 
is increased or decreased in accordance with the estimated error. This 
guarantees a minimum number of calculations for a certain degree of 
accuracy which is determined by the accuracy of the computer. 

Theories of convergence and inherent stability of the differential 
equations (Keller, 1968) cannot be examined for this problem. However, 
numerical experimentation shows that the solution converges as though 
these theories apply. 



CHAPTER V 


EXPERIMENTAL DATA 
Introduction 


This chapter summarizes all the test data and experimental in- 
formation obtained with TMAC electrolyte. These data will be re- 
quired for several reasons. First, they will be used to evaluate 
the parameters which cannot be estimated independently such as the 
mass transfer coefficients and gas-electrolyte interface area. Second, 
they will be compared with the values predicted by the model after the 
required coefficients have been calculated. This will indicate how 
close the model simulates the actual processes taking place in the 
cell. Fundamental property tests were also made on the electrolyte, 
electrodes, and matrix separately. These will help in justifying some 
of the assumptions made here and to explain some of the obtained re- 
sults. 

Performance Data 

Two sets of performance data have been obtained; analytical cell 
data and full-scale cell data. The first series of tests were run to 
obtain detailed information on cell operation with TMAC and to com- 
pare these results to those obtained with CS2CO2 electrolyte. The 
second series was run to test the long-duration performance under 
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design conditions. The design and measurements in each case are des- 
cribed in the following: 

1) Analytical cell test 

The analytical cell was constructed from fuel -cell pair hardware 

2 

to give an electrode area of 1/24 ft with an airflow path length of 
6 inches. The matrix consisted of one layer of Tissuquartz (18 mils) 
sandwiched between two layers of fuel cell asbestos (20 mils each) 
and the whole compressed to 25 mils. The electrodes were PPF anode 
and DS 16-0 cathode (coded names). Air temperature was 70°F. The 
main test results included CO 2 removal efficiency, open-cell voltage, 
cathode and anode overvoltages, anolyte and catholyte pH's and IR 
drop across the matrix. The parameters studied were current density 
(12-30 A/ft^), pQQ (0.25-4 mm Hg), air velocity (5-15 ft/sec) and 
air relative humidity (33.6 - 11 . 1 %). Each test condition was run for 
approximately two days with at least three data reading; taken each 
day. These data were reported by Huddleston and (1973). 

2) Full-scale cell tests 

Several tests have been made on the full cell by Hamilton Standard. 

The most recent lasted for 90 days and involved taking data for operating 

2 

currents and inlet air P^q levels over a limited range (5-18 A/ft 

2 3 

and 1-3 mm Hg). The air flow rate was fixed at 10 ft /min per cell 
which corresponds to an air velocity of 6.6 ft/sec. Mo other variation 
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of test conditions was made except for the decrease of the inlet air 
dew point from 60°F to 40°F in the last stage of the program. CO 2 
removal efficiency and cell voltage were the measured output data. 

The data were recorded twice every day and were reported by Huddleston 
and Aylward (1975 a,b). The variation of the steady-state removal 
efficiency with current density and P^q was summarized in a final 
report (Huddleston and Aylward, 1975 c) to NASA. 

Special Tests 


1) Electrolyte tests 


The first investigation of the TMAC electrolyte (Huddleston and 
Aylward, 1972) was to determine the approximate solubility and maxi- 
mum dew point depression for the carbonate, bicarbonate, and hydroxide. 
The solubility of these compounds at room temperature are 52 wt %, 88 
wt %, and 36 wt %, respectively. The depression of the dew point of 
these saturated solutions are 32°F, 55°F, and 28°F, respectively. The 
electrochemical stability of these compounds decreases in the order: 
NMe^HC 03 > (NMe^) 2 C 03 > NMe^OH. 

The important properties of the three electrolytes were reported 
by Aylward (1974). The density ineasurements were accurate to four 
significant figures; viscosity measurements v/ere accurate to three 
significant figures; specific conductivity were accurate to three 
significant figures; and electrolyte-water vapor pressures were accu- 
rate to +0.5%. 
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The absorption of CO 2 in TMAC in the presence of a catalyst was 
tested. It was demonstrated qualitatively (Huddleston and Aylward, 

1973) that small amounts of the enzyme carbonic anhydrase at very low 
levels are effective in accelerating the absorption reaction in TMAC. 
However, addition of this catalyst to the electrolyte in the analyti- 
cal cell tests showed no significant improvement in the removal effi- 
ciency. This indicates that the CO 2 , OH" reaction is not rate con- 
trolling under normal operating conditions. 

An error analysis study (Russell, 1973) showed that the measured 
CO 2 removal -efficiencies are accurate to +4%. 

2 ) Electrode tests 

The selection of DS 16-0 cathode for the TMAC cells came after a 
series of tests on different cathode designs for cells with CS 2 CO 2 
electrolyte (Huddleston and Aylward, 1973). The significant improve- 
ment of CO 2 removal efficiency with increasing electrode porosity was 
explained by the increase of the total available gas-catholyte inter- 
face area and the ‘increase of the intra-cathode CO 2 gas mass transport 
rate. This was also explained by the increase of the electrode elec- 
tronic resistance with porosity which results in a significant drop-off 
in the current density remote from the current collector screen. This 
results in a decrease of the electroosmotic pumping of the electrolyte 
towards the air side while the capillary force, which is independent 
of position from the current collector screen, returns the electrolyte 
towards the matrix (Huddleston and Aylward, 1973). Floating-electrode 



42 


tests showed that the cathode overvoltage increases with porosity as 
expected. 

The anodes are of the PPF type containing 15 mg of catalyst 
(platinum) per cm . Floating electrode tests showed that they have 
the lowest overvoltage among tested anodes. In the analytic cell 
tests the anode overvoltage obtained is typically one to two hundred 
millivolts greater than that obtained with the floating electrode 
tests. This was explained by the decrease of anolyte concentration 
which was shown to increase the anode overvoltage. The anode polariza- 
tion consists mainly of concentration and ohmic components; the activa- 
tion component is very small at the operating current densities. Char- 
acteristic X-ray examination of the anode showed that over long periods 
teflon migrates to the matrix side while platinum migrates to the H 2 
side. This results in a decrease in the wicking rate of the anode 
since teflon is hydrophobic. This is responsible for the long-term 
power decay rate and may be remedied by current interruption and N 2 
purge (Huddleston and Aylward, 1973). The problem of short-term power 
decay is attributed to the formation of a reduced CO 2 complex layer 
because of the low pH of the anolyte; the anode is reactivated by 10 
sec. N 2 purge at one hour intervals (Aylward, 1975). 

3) Matrix tests 

The effect of matrix thickness was tested for CS 2 CO 2 electrolyte 
on the analytical cell (Huddleston and Aylward, 1973). One, two, and 
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four layers of asbestos were tested at the same compression with thick- 
nesses of 12,5, 25, and 50 mils, respectively. The removal efficiencies 
(based on CO^ transport) at 18 A/ft^ were 66%, 98%, and 106%, res- 
pectively. Increasing the matrix thickness lowers the anolyte pH 
which was found to increase the chemical attack of the asbestos at the 
interface between the matrix and the anode causing the matrix to ad- 
here to the anode. With TMAC electrolyte the anolyte pH does not 
show as strong a dependence on matrix thickness as with CS 2 CO 2 . This 
reduces the problem of matrix sticking, as observed experimentally 
(Huddleston and Aylward, 1973). 

For the matrix materials used in the present cell several tests 
were made to determine the optimum compression, which results in a 
minimum internal cell resistance. The electrical resistance of the 
electrolyte matrix combination was measured using a special micrometer 
cell (Huddleston and Aylward, 1973). The matrix factor, which repre- 
sents the ratio of the conductance of the matrix-electrolyte combina- 
tion to that of the free electrolyte, was found to fit the following 
empirical equation 


" df T + 3 


[5.1] 


where d is the matrix thickness, f is a compression factor equal to 

where d° is the maximum compressed thickness defined as the 
original uncompressed thickness multiplied by (1-e), where e is the 
porosity, t is the matrix tortuosity factor, and 3 is an empirical 
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constant dependent on the matrix material and electrolyte. No equa- 
tion could fit data obtained for composite matrix material. 


CHAPTER VI 


MODEL PARAMETERS 
Introduction 


In this chapter the various parameters required by the model 
are estimated from electrolyte properties, theoretical and empiri- 
cal equations and performance test data. The detailed calculations 
and tables are given in Appendix A. Unless otherwise specified, 
the values used here are at 25“C (77°F). 


Activity Coefficients 

Activity coefficients of electrolytes can be obtained from the 
vapor pressures of their solutions by (Robinson and Stokes, 1959); 

rm 


In = ( 4 ,- 1 ) + 


(4>-l )d In m 


[6-1] 


0 

where the osmotic coefficient ({) at molality m is given by 


4) = (-1000/vmM)ln(P/P°) [6-2] 

Here v is the number of ions in solution per one molecule of 
electrolyte, M is the molecular weight of the solvent, and P and 
P° are the vapor pressures of the solution and solvent, respectively 
Evaluation of the integrand in equation [6-1] requires measure- 
ment of the vapor pressures of dilute solutions. Since accurate 
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measurement at low concentrations is difficult, activity coefficients 
have been tabulated for only a limited number of electrolytes 
(Robinson and Stokes, 1959). Several methods have been described for 
the calculation of activity coefficients where data are available 
only for concentrated solutions. Meissner and Tester (1972) con- 
structed generalized charts for the calculation of the activity co- 
efficients for certain classes of strong electrolytes . Earlier, 
Lietzke and Stoughton (1962) computed activity coefficients from 
osmotic coefficient data using a four-parameter equation containing 
the Debye-Huckel term for the electrostatic contribution and addi- 
tional terms in ionic strength. They later presented a three-para- 
meter equation in which the non-ideality of electrolytic solutions 
is attributed to both the Debye-Huckel interactions and randomized 
fused-salt behavior (Lietzke et al . , 1968). 

Tetraalkyl ammonium salts show significant deviation from the 
Delye-Huckel limiting law in dilute solutions (Pethybridge and 
Prue, 1968; Wen, 1972). Thus, the above-mentioned methods become 
unsuitable. In addition to the electrostatic and simple hydration 
interactions, mutual salting-out and structural hydration interactions 
must be considered. The mutual salting-out effects result from the 
fact that other ions distributed about a given ion occupy a volume in 
the polar solvent medium but with less polarization than an equal 
volume of solvent. In the case of large univalent ions such as NR^, 
these effects are quite significant and can be calculated from the 
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partial molal volumes of the cations and anions (Desnoyers and Conway, 
1964). On the other hand, the structural interactions, which result 
from the compatibility or incompatibility of the hydration of the 
ions (Desnoyers et al., 1969) has been described only qualitatively. 
These interactions are estimated here by comparison with other 
properties of the solution. 

Consider the model of Desnoyers and Conway (1964), wherein the 
activity coefficients are presumed separable into several component 
terms as follows: 

log y = log + log + log + log y^^ [6-3] 

where the subscripts el, h, so, and st refer to the electrostatic, 
hydration, mutual salting-out, and structural interactions respec- 
tively. Equation [6-3] contains one more term, logy^^, than the 
Desnoyers -Conway equation. Equation [6-3] is assumed to apply from 
infinite dilution up to the lowest concentration for which osmotic 
coefficient data are available. Each term can then be estimated us- 
ing appropriate equations, as described below. 

The Debye-Huckel equation was used to calculate the electro- 
static term, y^-j . To account for ion-pair formation, experienced 
even in strong electrolytes, Neff (1970) used the following form of 
the Debye-Huckel equation replacing the ionic concentrations with 
the activities: 
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- 0.5089 1 z+z_ I (Ea.z^.^) 

log Yei = ! 

1 + 0.3291 L/i (Sa.z^.^) 

This equation was tested for NaCl solutions where it was found that 
the distance of closest approach, L, was close to the sum of the 
crystallographic radii of Na^ and Cl“. This equation is here assumed 
applicable for the tetramethyl ammonium electrolytes with the para- 
meter L calculated from the crystallographic radii. 

The hydration term, y^, can be calculated either by the equation 
developed by Stokes and Robinson (1948) or by the following equation 
developed by Glueckauf (1955) 

= (?-v)log(l+0.018xm)/v - (§/v)log(l-0.018^m) 

+ 0.018x(x+?-v)m/[2.303v(l+0.018xm)] [6-5] 

where ^ is the hydration number and x is the ratio between the 

apparent molar volume of the electrolyte and the molar volume of 

water. This equation is preferable to the Robinson-Stokes equation 

since it does not contain the unknown term In a . 

w 

The mutual salting-out term, y^^, can be calculated for an 
electrolyte "MA" from the following equations (Desnoyers and Conway, 
1962): 

for 1:1 electrolyte 


logYso = (7.156 C/4605)[V_/(3L^)+V^/L^+V^/(3L_) + V./Lj 


[ 6 - 6 ] 
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for 1:2 electrolyte 

logYjo = (7.156 C/4605)[4V_/{5L^) + 2F^/L^+VL_+4V_/L_] [6-7] 

where F. and L. are the partial molar volume and 1on1c radius for the 
ion j. To account for ion-pair formation, the same procedure used 1n 
deriving equation [6*4] was also used in modifying equations [6-'6] 
and [6-7], i,e. the concentration was replaced by the activity and 
the ionic radius was taken as the crystallographic radius. 

The final term, can be calculated by difference for electro- 
lytes where activity coefficients are known. According to the 
structural interaction model (Desnoyers et al., 1969) this term re- 
sults from four types of interactions which depend on the structural 
properties of the ions. This term is calculated here for the 
tetramethyl ammonium halides from their tabulated activity coefficients 
(Lindenbaum and Boyd, 1964; Wen et al . , 1S66; Levien, 1965) and use 
these values to estimate similar terms for (NMe^)C 02 , HMe^HCO-,, and 
NMe^OH. The strong relation between the activity coefficients of 
electrolytes on the one hand and the viscosity B-coefficients and 
entropy of hydration on the other hand have been indicated by Gurney 
(1962). Desnoyers and Perron (1972) showed that the term (B- - 

J 

O 

0.0025 V.) indicates whether the ion is 'structure-maker' or 
J 

'structure-breaker' depending upon whether this term is positive or 
negative,respectively . Here B. is the ionic viscosity coefficient in 

J 

the Jones-Dole equation 
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(n-n°/ri°) = 1 + + z_B_)C [6-8] 

This parameter (B. - 0.0025 V.) was used as a basis of comparison. 

J J 

In the case of NMe^HCO^. a value of B was interpolated on the basis 
of the entropy of hydration, s. The detailed calculations and tables 
are shown in Appendix A. 

The calculated mean activity coefficients of the three electro- 
lytes are shown in Fig. 4. The activity coefficients in multicom- 
ponent mixtures can be calculated using the method of Meissner and 
Kusik (1972): 


= l/ 2 (logr^j + Y^^logr^^ + Y^2^ogr^2 ^^3109^43) ^ 6 - 9 ] 

where Y^-j , Y^ 2 ’ Y ^2 ionic strength fractions of the 

corresponding electrolytes, m refers to multicomponent mixtures, and 
is the reduced activity coefficient defined by the equation 


4j 






[ 6 - 10 ] 


Equilibrium and Reaction Rate Constants 

The equilibrium constants at infinite dilution for reactions 
[2-1] and [2-2] are given by the following equations (Harned and 
Owen, 1958]: 

log K.° = -log - (3404. 71/T) + 14.8435 - 0.032786 T [6-11] 

log K.. = -log - (2902.39/T) + 6.490 - 0.02379 T 


[6-12] 
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where is the ionization constant of water, which has been re- 
ported by Robinson and Stokes (1959). The values of can be 
approximated by the following equation: 

- log = 4471. 33/T - 6.0846 + 0.017053 T [6-13] 

The equilibrium constants used in the mathematical model were 
calculated from their values of infinite dilution multiplied by the 
appropriate activity coefficient ratio. This ratio was calculated 
using equations [6-9] with average concentrations estimated for the 
three regions; catholyte, anolyte, and bulk electrolyte. These 
values do not significantly affect the calculated CO 2 flux or the 
calculated HCO^ concentration at the anode. This eliminates the need 
for an iteration procedure to calculate the exact concentrations. The 
details of the calculations are shown in Appendix A. 

The rate constants for reactions [2-1] and [2-3] at infinite 
dilution are given by the equations (Pinsent et al . , 1956): 

log k. = 13.635 - (2895/T) . [6-14] 

log k... = 329.850 - 110.541 logT - (17265. 4/T) [6-15] 

Henry's Law Constants 

Henry's law constants were estimated from their values in water 
(Perry et al . , 1963) multiplied by a factor dependent-on the ionic 
strength of the electrolyte (Van Krevelen and Hoftijzer, 1948). At 
infinite dilution the Henry's law constants are approximated by the 
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equation 

H° = 10^(13.32257 - 0.132516T + 3.198135xl0‘^ T^) [6-16] 

where H is in (nm Hg) il/gmole. Equations [2-10] and [2-11] were used 
to calculate the values of H for the catholyte and anolyte, taking 
i^ = 0.030, i_ = 0.021, and ig = -0.01 (Van Krevelen and Hoftijzer, 
1948). The value of i^ used here is that value reported by these 

"f* ■}' 

authors for NH^ ion since there was no value reported for NMe^. 

Gas-Phase Mass and Heat Transfer Coefficients 

The mass transfer coefficients in the air zone were calculated 
from the limiting slope of Ng against (Pg )j at (Pg )j = 0 (Fig. 5). 
In that case Cn vanishes and equation [4-5] reduces to 

dNg 

h, = - (V /A)ln[l - (ART/V )(gp-L) ] [6-17] 

c 5^ 

From this equation hg can be calculated for each air velocity. Test 
data are available for the analytical cell at air velocity of 10 ft/ 
sec and (Pg )j of 0.25, 1, and 4 mm Hg. The calculated mass transfer 
coefficient in that case is 233 ft/hr. At air velocity of 5 and 15 
ft/sec data are available only at (Pg )j = 0.25 mm Hg. These are 
insufficient to calculate the mass transfer coefficient by the limit- 
ing-slope method. 

For air velocities in the range encountered in the cell, the 
flow is laminar. Two equations are available for mass and heat 
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transfer coefficients in rectangular ducts (Knudsen and Katz, 1958). 

For uniform wall temperature the heat transfer coefficients are 
calculated by; 

Nu = Nujl + (0.003+0. 039 [6-18] 

For uniform heat flux the equation changes slightly to 

Nu = NuJl + (0. 003+0. 019 6 -,/ 62 )RePrdg/ 63 ] [6-19] 

where Nu is Nusselt number, Nu^ is its asymptotic value, 6 -| and 6 ^ 
are the dimensions of the rectangular section, 63 is the length of 
the air duct, d^ is the equivalent diameter of the air duct and Pr 
is the Prandtl number. In the case of mass transfer Nu is replaced 
by Sh and Pr by Sc where Sh is the Sherwood number and Sc is the 
Schmidt number. 

In our case the mass transfer coefficient in the bulk gas phase 
can be calculated by equation [6-18]. The value of h^^ at an air 
velocity of 10 ft/sec is 246 ft/hr calculated at an average cell 
temperature. Comparing these values with the value obtained pre- 
viously from the limiting-slope method (233 ft/hr), it appears that the 
mass transfer resistance in the gas phase within the electrode pores 
is negligibly small. Thus, the overall mass transfer coefficient, 
hg^, is about 94% of the value calculated by equation [6-18]. Assum- 
ing that the overall mass transfer coefficient is 94% of the theoreti- 
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cal value, equation [6-18] can be used to calculate hg^. For the 
present design this reduces to the following equation for hg^ calcu- 
lated at an average cell temperature: 


hg^ = 0.94 {112 + 8.84 V^) 

3 

where hg^ is in ft/hr and is in ft /mm. 


[ 6 - 20 ] 


This equation applies for both the analytical and full-scale cells 
since the dimensions of the air ducts are the same for both designs. 
Other mass transfer coefficients are calculated from equation [6-18] 
while the heat transfer coefficients used in the energy balance calcu- 
lations (Lin et al . , 1974) are calculated from equation [6-19]. 

Transport Properties in the Electrolyte 

The diffusion coefficients of (NMe^)2C02, NMe^HCOg, and NMe^OH 
were calculated from the activity coefficients and viscosity data 
(Fig. 6). The detailed calculations are given in Appendix A. The 
values of the diffusion coefficients are shown in Fig. 7. For each 
electrolyte the diffusion coefficients exhibit a minimum and than a 
shallow maximum as the concentration increases. This complex concen- 
tration dependence is typical of diffusion in electrolytic solutions 
(Bhatia et al . , 1957). From these values the individual ionic co- 
efficients were calculated by the method described in Chapter II. The 
ionic mobilities were also calculated from the equivalent conductance 






C, g mole / il 


Fig. 7 Diffusion Coefficients of {NMe.)pC0os NMe^HCO-,, and 
NMe^OH at 25 ^ ^ ^ j 
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data (Fig. 8) by the method also described in Chapter II. To accommo- 
date computer calculations, the viscosities, ionic diffusion coeffi- 
cients, and ionic mobilities were fit to polynomials. These equations 
are given in Appendix A. Although these equations are empirical and 
have no real theoretical significance, they are suitable for calcu- 
lations because the model is rather insensitive to moderate variations 
in the values of these parameters. 

The diffusivities of CO 2 in the electrolyte were estimated from 
their values in water (Perry et al., 1963) using the following empiri- 
cal equation which was obtained for the diffusion of CO 2 in liquids 
(McManamey and Woollen, 1973): 

D = D°(nVri)°‘^^ [6-21] 

Matrix Factor 

The matrix, as mentioned earlier, consists of two sheets of fuel- 
cell asbestos with an uncompressed thickness of 40 mils and one sheet 
of Tissuquartz between with an uncompressed thickness of 18 mils. The 
matrix is compressed to a final thickness of 25 mils. From Fig. 39 of 
Huddleston and Aylward (1973), the corresponding loading is 128 psig. 
At this loading the same two sheets of fuel cell asbestos can be com- 
pressed only to 23 mils. This indicates that the Tissuquartz layer is 
compressed to approximately 2 mils. The matrix labyrinth factors for 
the asbestos and the Tissuquartz layers at this compression are 0.342 
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and 0.314, respectively, estimated from the conductivity measurements 
(Huddleston and Aylward, 1973). The average matrix factor is then 
about 0.340. The mobilities and diffusion coefficients are multiplied 
by this factor to obtain the effective values in the matrix zone. 

Gas-Electrolyte Interface Area 

Two parameters must be obtained from test data: the air- 

catholyte interface area and the hydrogen -anolyte interface area. 
Experimental methods based on adsorption measurements measure the 
total internal solid surface area (Srinivasan and Gileadi, 1968). 

On the other hand, dilatometry measurements can provide an estimate of 
the electrolyte film area extending inside hydrophilic electrodes which 
provides the active surface for electrode reactions (Moren, 1975). 
Neither of these methods can give an estimate of the true gas-electro- 
lyte interface area in partially wetted electrodes as those used in 
the HOC cells. 

The gas-electrolyte interface area is determined by several 
forces acting on the electrolyte surface; gas pressure, electroosmotic 
forces and capillary forces. The presence of the asbestos matrix in 
the vicinity of the electrode matrix adds another complication be- 
cause of the difference in their porosities and their affinity to the 
electrolyte. Thus, the true gas-electrolyte interface area can only be 
obtained empirically from test data. 
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The calculated N^j is independent of the chosen value of 
but does depend on the chosen value of S^. The latter parameter 
determines the concentration of OH" ions in the catholyte which deter- 
mines, in turn, the rate of CO 2 absorption reaction and the fraction 

2 - 

of the current carried by the CO^ and HCO^ ions in the matrix zone. 

The value of can be chosen so that the calculated CO., flux matches 

the experimental values for both the analytical and full-scale cells. 

In equation [4-9] the term C^/(K^C 2 ) is in the range 10“'*^-10”^^ 

+ -9 -6 

gmole/£ whi le IS in the range 10 -10 gmole/Jl. The equation can 
thus be simplified to the following form: 





[ 6 - 22 ] 


This equation suggests that, with ^5c calculated by equation [4-5], 
the value of Ngy is proportional to (C 2 ) ' for constant values of S^, 
k^, and D^. At a given total ionic strength, the last three parameters 
are constant. Figure 9 shows some experimental data obtained with the 
analytical cell for a constant air relative humidity, or in other words, 
at the same catholyte ionic strength. These values are compared with 
those calculated using the same value of the parameter (S^k^D^). The 
agreement as shown in Fig. 9 is unsatisfactory since this parameter is 


presumably independent of (P 5 (.)];- Better fitting can be obtained if 

1 /2 

Ngj is assumed to be proportional to C 2 rather than (C 2 ) ' and the 
active area changing with the current density as shown in Fig. 10. In 


that case the calculated CO 2 flux fit the experimental data within ex- 
perimental error as shown in Fig. 11. 



C02 REMOVAL EFFICIENCY (%) 






REMOVAL EFFICIENCY {%) 


•- CALC^^-v. i 
f ''n j\ 

z'' \V 

/ T 

'I 

f X\\ \ 








CURRENT DENSITY, A/FT 


Fig. 11 Test Data Compared with Simulation Re- 
sults if the CO 2 Flux is Proportional to 
[OH] and the Air-Catholyte Area is De- 
pendent on Current Density, Analytical 
Cell, Air Inlet Temp. = 70°F, Air Dew 
Point = 50°F 
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To explain this we must reconsider the nature of the chemical 
reaction between CO 2 and OH" ions, which is typical of ion-molecule 
reactions (Frost and Pearson, 1961). The dependence of the rate 
constant on the ionic strength of the solution as given by equation 
[2-6] is resolved into the following equation: 

log k = log k° + + £q|., [6-23] 

where the l terms stand for the dependence of the activity coefficient 
of the corresponding species on the ionic strength. The subscript I 
represents the activated complex which is probably formed as a 
reaction intermediate. It is interesting to note the similarity 
between this equation and equation [2-10] which describes a com- 
pletely different phenomenon, i.e., the solubility of a gas in an 
electrolytic solution. It is the term which we need to reconsider 
here. This term results from the assumption that the activity co- 
efficient of an ion is represented by the following equation: 

logy,- = -0.5809 (1+0.3291 L [6-24] 

J J 

This equation is a modified form of the Debye-Huckel equation which 
is assumed to apply for both the OH" ions and the activated complex. 

In the case of NMe^OH the activity coefficients in their solutions 
increases sharply with concentration which results in some dependence 
of their activity coefficients on C 2 as calculated for the multicom- 
ponent catholyte by equation [6-9]. The activity coefficient of 
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the activated complex is probably of the same order of magnitude as 
that of NMe^HCO^ which does not change much with concentration. The 
result is that the term is dependent to some extent on C 2 

which results in the dependence of k^. on C 2 . The remaining effect 
of the ionic strength on the rate constant is accounted for by the 
term £qq which is represented by the following equation (Frost and 
Pearson, 1961). 


^°9"^C02 " ^C02^ 


[6-25] 


The activity coefficient of CO 2 in an electrolytic solution is also 
represented by the following equation (Harned and Owen, 1958): 


5-CO2 h/h° 


[2-26] 


Combining these equations with equation [2-11] provides an easy method 
to predict the effect of the ionic strength on the reaction rate con-- 
stant . 

The dependence of the gas-electrolyte interface area on the 
apparent current density (current per unit electrode geometric area) 
is less clearly understood. It is known that the equilibrium of the 
electrolyte meniscus within the electrode pores depends upon the 
pressure differential between the gas and electrolyte, the capillary 
forces, and the electroosmotic forces. In regions of high current 
densities close to the current collector screen, the electroosmotic 
forces push the electrolyte towards the air side while the capillary 
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forces force the electrolyte back into the matrix structure in regions 
of low current density remote from the current collector screen. It 
is believed that any factor which increases the drop-off in current 
density remote from the current collector screen, such as increasing 
the cathode electronic resistance by increasing its porosity, signifi- 
cantly increases the COg removal efficiency (Huddleston and Aylward, 
1973). This may explain the change of the air-catholyte interface 
area with current densities as required to simulate the experimental 
data. 

The CO 2 removal efficiency in the full-scale cell tests were 
significantly lower than those obtained with the analytical cell at 
the same conditions. The possible effect of the flow rate and 
pressure and the cell temperature were investigated experimentally. 

It was found that these parameters have little effect on the CO 2 re- 
moval efficiency (Huddleston and Aylward, 1973). The mathematical 
model correctly predicts this lack of sensitivity. 

The difference between the air flow rates in the analytical and 
full-scale cells has little effect on the CO 2 removal efficiency at 
the baseline conditions (P^q = 1-3 mm Hg). The only apparent dif- 
ference between the analytical and full-scale cells is the geometric 
configuration of the cell electrodes as shown in Fig. 12. Both cells 
have the same air path length but the width of the electrodes in the 
full-scale cell is 24 times of that in the analytical cell. It seems 
that this results in a difference in the current density distribution 
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Fig. 12 Electrodes of the Analytical and Full-Scale Cells 



70 


in the cathode with a corresponding change in the balance between 
the electroosmotic and capillary forces which directly affect the 
gas-electrolyte interface. 

The variation of the air-catholyte interface area, as shown in 
Fig. 11 was empirically fit by the following equations using the 
method of Churchill and Usagi (1974): 

For the analytical cell 

W " + (0-06 [6-28] 

for the full-scale cell 

W = “ ■'•24/[1 + (0.06 i)-8]0.5 [6-29] 

These equations represent the relative values of at a given current 

density compared to that at a characteristic current density chosen 

2 

arbitrarily at 18 A/ft with the analytical cell. The absolute 
values of cannot be obtained after making the assumption that Ngj 
is proportional to ^ 2 c- Instead, the parameter (f^5j/C2c^ 
calculated which is directly proportional to S . The variation of S 
with current density can then be accounted for in the mathematical 
model by multiplying the parameter (^^ 57 /^ 2 ^) appropriate ratio 

(Appendix A). 

The magnitude of S cannot be calculated accurately from the test 

d. 

data because the value of and at the anode cannot be determined 
experimentally with sufficient accuracy. However, the calculated CO^ 
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transfer rate and the concentration of C-j at the anode are almost in- 
dependent of the chosen value of S^. Thus, it is possible to use an 

Q 

approximate value of in the mathematical model ignoring any 
possible change with current density or with any other relevant para- 
meter. The concentration of the hydroxyl ions can be estimated from 
the pH measurements of the bicarbonate solutions (Abdel -Salam, 1974) 
after assigning a reasonable value to the activity coefficients of 
the OH" ions. The calculations are shown in detail in Appendix A. 



CHAPTER VII 


MODEL PREDICTIONS 


Introduction 


This chapter compares the performance of the concentrator with 
that of the model over a wide range of the operating variables. It 
also predicts the performance of the cell when one of the design 
parameters, e.g. the matrix thickness, is changed from the used 
value. The limits of stable cell operation in low humidity environ- 
ments can be determined from the concentration distribution obtained 
by the model . 


Effect of P 


CO, 


At a constant current density, increasing P^q in the inlet cabin 
air increases the CO 2 removal rate. This is apparent from Figs. 13 
and 14 which compare the experimental and calculated rates for the 
analytical and full-scale cells. The increase is rather sharp in 
the low P^g region becoming much slower as Pgg rises above the de- 
sign baseline condition (3 mm Hg). 

In the low Pgg region (< 0.25 mm Hg) the process is gas-phase 
mass-transfer controlled and the increase of N^T is approximately 
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Fig. 13 Effect of at the Air Inlet on the CO 2 Re- 
moval Rate in the Analytical Cell (the solid 
lines represent Simulation Results) 



m m IH g 



Fig. 14 Effect of in the Inlet Air on the CO 2 Removal 

Rate In the Full-Scale Cell (the solid lines repre 
sent Simulation Results) 
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linear with P^q as is evident from equation [4-4]. As level 

increases the ionic transport in the matrix zone becomes the con- 

2 _ 

trolling step. In this zone the current is carried by the CO^ , 

HCOg, and OH" ions. Increasing the CO 2 removal rate at a constant 

current density results in a lower fraction of tho current carried by 

OH" ions. A lower OH" flux results in a smaller concentration gradient 

across the matrix which means a lower OH" concentration at the cathode 

since there is no significant change in the anolyte pH. Figure 15 

compares the calculated OH" concentration at the cathode in the 

2 

analytical cell at a current density of 18 A/ft with the pH measure- 
ments obtained at different inlet Pgg . Although these measurements 
are crude for this purpose, not only because of the low accuracy of 
the measurements but also because taking a representative sample of 
the catholyte is hard to achieve, they show the same trend as predicted 
by the model. 


Effect of Current Densitv 


For the same value of P^g , increasing the current density in- 
creases the CO 2 fluxo However, the CO 2 removal efficiency (current 
efficiency), which is proportional to the ratio between and I, has 
a maximum value at a certain optimum current density. The optimum 
current density is higher for higher inlet Pgg . This is shown in 
Figs. 11 and 16 for the analytical and full-scale cells. 
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Fig. 15 Effect of on the Hydroxy] Ion Concentration in the 
Catholyte of the Analytical Cell 
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This dependence of the current efficiency on current density is 
related to the OH” concentration in the catholyte. Increasing I is 
accompanied by increase in 82^, which, although increasing the rate 
of CO2 absorption reaction, also increases the OH" flux in the 
matrix zone. The current efficiency is proportional to the ratio 
between (N-j+N^+Ng) and (N^+N2+2M2). Although these sums increase 
with current density, their ratio has a certain maximum at a certain 


82^ value. This optimum concentration is also dependent on (p5^)j 
The calculated values of 82^ at P^q of 1 mm Hg for the analytical 
cell are compared with the measured pH values ir Fig. 17. 


Effect of Air Velocity 

The effect of air velocity on the CO2 removal rate is shown in 
Fig. 18 for the analytical cell. The calculated values are compared 
with the experimental values. Because of the gas-phase mass-transfer 
limitations in the low P^q region, increasing the air velocity in- 
creases the CO2 removal rate in that region. This is also evident 
from the lowering of the catholyte pH with the increase in air velocity 
(Fig. 19 ). At higher P^g the ionic transport becomes more important, 
and hence, the increase in air velocity has no significant effect in 
that region. 


Effect of Matrix Thickness and Compression 


The matrix thickness and compression are chosen on the basis of 
the internal cell resistance and other mechanical considerations. How- 


mole / 
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Fig. 19 Effect of Air Velocity on the Hydroxyl 
Ion Concentration in the Catholyte of 
the Analytical Cell 
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ever, it is desirable to predict the effect of changing the matrix 
thickness and compression on the CO^ removal rate. Figure 20 shows 
the effect of changing the matrix thickness at the same compression. 
Increasing the matrix thickness increases the CO^ removal rate up 
to a certain limit after which the CO 2 rate decreases. The initial 
increase is attributed to the decrease in back diffusion of the HCO^ 
ions from the anode to the cathode. The effect of increasing the 
matrix thickness on increasing the selectivity of the 0H“ ions as 
charge carriers, which have higher mobilities and diffusion coefficients 
than the other ions, becomes apparent beyond a certain thickness vvhich 
explains the decrease in the CO 2 removal rates afterwards. 

The same reasoning can explain the effect of matrix compression 
which is shown in Fig. 21. Of course, the increase of the CO 2 removal 
rates on increasing the matrix thickness or compression in the region 
around design values is at the expense of the internal cell resistance 
and on the anolyte pH which affects the anode overvoltage. Any change 
in these parameters should also consider these effects on the power re- 
quirements. 

Effect of Air Relative Humidity and Temperature 

The model predicts a slight increase of the CO 2 removal efficiencies 
as the air relative humidity decreases. Test data show no significant 
effect of this parameter beyond the experimental error limits. This 
predicted change may be because the change of the electrolyte surface 
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Fig. 20 Effect of Matrix Thickness on the C0« Removal 
Rate in the Analytical Cell 
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tension with concentration is another factor to be considered in 
determining the gas-electrolyte interface area. Obviously, this 
cannot be accounted for here because of the empirical nature of the 
determination of the gas -electrolyte interface area and the lack of 
quantitative correlation with surface tension. 

Temperature also has little effect on the CO 2 removal efficiency 
over a limited range. The effect of increasing the reaction rate 
constants k^. and the diffusivity of COg in the catholyte with tempera- 
ture seems to be counteracted by the decrease in the selectivity of 
the CO^ as charge carrie''.' as the viscosity of the electrolyte de- 
creases in the matrix zone. 

More important is the change of the electrolyte concentration as 
the difference between the cell temperature and inlet air dew point 
increases. It is important that the concentration of the electrolyte 
at any point does not reach the solubility limit; otherwise the con- 
centrator undergoes what is known as "cell drying" which causes the 
hydrogen to cross over from the anode chamber to the cathode chamber. 
Since the catholyte consists mostly of (NMe^)2C02 and the anolyte of 
NMe^HCO^, the pure component solubilities can be used to check on 
that condition. Figure 22 shows the distribution of the concentration 
of the electrolyte at different points in the concentrator at air 
temperature of 70°F and dew point of 60”F. The cell temperature and 
water content of the air stream at any point along the air path is 
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determined by the current density, air flow rate, and cell voltage. 

As seen from Fig. 22 no precipitation is predicted in the concentrator 
at any point. Still at a lower air humidity, precipitation may com- 
mence at one or the other point in the concentrator. It is clear 
from Fig. 22 that precipitation could commence at the cathode rather 
than the anode because of the low solubility of (NMe^) 2 C 02 compared 
with that of NMe^HCO^. In the example of the computer output shown 
in Appendix B this precipitation is indicated at the cathode. 

Effect of Hg Flow Rate and Anode Active Area 

The effect of changing the H 2 flow rate and the anode active area 
were investigated. No significant change of the CO 2 removal efficiency 
resulted on changing their values by several orders of magnitude. A 
slight decrease in the anolyte calculated pH was obtained on de- 
creasing either the H 2 flow rate or the anode active area. The effect 
of the H 2 flow rate is the same as obtained experimentally (Huddleston 
and Aylward, 1973). 
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Fig. 22 The Concentration Distribution of the CO^, HCO^ 
and NMe^^ in the Analytical Cell at the Follow- 
ing Conditions: = 4 mm Hg, Current Density = 

18 A/ft^, Air Velocity = 10 ft/sec, Air Temp. = 
70°F, Air Dew Point = 60'’F 



CHAPTER VIII 


CONCLUSIONS AND RECOMMENDATIONS 

A mathematical model has been developed for simulating the 
steady-state performance in electrochemical CO^ concentrators 
which utilize (NMe^)2C02 (aq.) electrolyte. This electrolyte 
which accommodates a wide range of air relative humidity is most 
suitable for one-man air revitalization systems. The model is 
based on the solution of coupled non-linear ordinary differential 
equations derived from mass transport and rate equations for the 
processes which take place in the cell. The boundary conditions 
are obtained by solving the mass and energy transport equations. 

A shooting-method is used to solve the differential equations. 

The present model displays two different features from the 
previous model of the CS 2 CO 2 units (Lin and Winnick, 1974). It 
has been found that better fitting of test data is obtained if the 
CO 2 flux across the cell, N^j, is proportional to the concentration 
of the hydroxyl ions in the catholyte, ^2c’ rather than the square 
root of this concentration. It is known that in concentrated 
solutions the value of k^. deF-c:''ds somehow on the activity coeffi- 
cients of the reacting OH" ions. The mean ionic activity coeffi- 
cients of NMe^OH obtained from vapor pressure measurements increase 
much more rapidly with concentration than those of alkali hydroxides. 
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This may explain the dependence of Ng-|- on € 2 ^ rather than 
The second important difference between the two models is that the 
assumption of a constant air-catholyte interface area, S^, at all 
current densities could not be extended to the present model. It 
was necessary to assume that the true interface area decreases with 
current density and that the decrease is more pronounced in the inter- 
mediate region (10-20 ASF). It was also necessary tc use different 
values of Sq for the analytical and full-scale cells to simulate the 
test data. It seems that the relatively open structure of the new 
cathodes used in these cells promotes the role played by the electro- 
osmotic forces in changing the electrolyte equilibrium within the 
cathode pores. These forces are dependent on the current-density 
and this may explain the dependence of on the current density. 

The model has confirmed the importance of the cathode design and 
catholyte pH in determining the CO 2 removal efficiency; other factors 
being of little importance in that respect. The CO 2 removal rate in- 
creases with current density but the current efficiency has a maximum 
value at a certain current density. The optimum current density is 


dependent on the value Pqq • The increase of the COg removal rate 

2 

with current density is more apparent as P^q increases. At low 
values of P^q •25mm) the removal rate is significantly Increased 
as the air velocity increases. However, the increase in the removal 


rate is rather limited at these low levels of P 



To attain the 


atmospheric CO^ concentration level (0.25 mm Hg) inside the cabin 



rather than the baseline design level (3 mm Hg), the number of the 
HDC cells must be approximately increased four to five times to re- 
move the same amount of CO^. 

In addition to calculating the CO 2 transfer rate at any set of 
specified conditions, the model also calculates the concentration 
distribution of all the important molecular species. These concentra- 
tions are difficult to measure experimentally because of the small 
dimensions of the cell. These calculations are important in several 
ways. First they can indicate if there is any precipitation expected 
in the cell under certain conditions. The concentration overpotential 
in the electrolyte can be estimated from the concentration distribu- 
tion of the (MMe^)^ ion. The OH" concentration at the cathode and 
the anode can also be helpful in estimating, in an approximate way, 
the electrode overpotential. From these components it may become 
possible to estimate the cell voltage beforehand and to understand the 
factors which contribute to a low cell voltage. 

The capability of the model lies in the fact that it is derived 
from fundamental mass and energy transport equations with the para- 
meters required by the model either independently measured or assigned 
reasonable values on the basis of previously developed equations. This 
gives us more confidence in extrapolation to off-design conditions than 
any empirical model. The only parameter that could not be obtained 
independently from test data was the air-catholyte interface area. 

The variation of this area with current density was easily identified 
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since the other parameters were fixed. With the growth of our under- 
standing of the nature of porous electrodes and the equilibrium 
between the gas and electrolyte within the electrode pores, a method 
may be developed to estimate this parameter from measurable proper- 
ties such as electrode porosity, pore diameter, and electrolyte 
properties. In that case a completely self-independent model can be 
developed which does not depend in any way on test data. 

Any attempt to improve the removal efficiency should concentrate 
on the cathode design and the factors which increase the gas-electro- 
lyte interface area. The geometric configuration of the electrodes 
may have an important effect in that respect as is shown by the test 
data of the analytical cell compared with the full-scale cell. The 
search for new electrolytes to accommodate wide ranges of air humidity 
should concentrate on higher members of the tetraal kyl ammonium carbo- 
nate series. Increasing the size of the radical is accompanied 
by an increase in y+ for the carbonate according to the model of the 
activity coefficients presented in Chapter VI. It should be expected 
that the higher members will have lower water activity, and thus, a 
greater value of water vapor pressure lowering. Further investigation 
of new electrolytes should concentrate on tetraethyl-, tetrapropyl- 
and tetrabutyl ammonium carbonates. 



NOMENCLATURE 


a activity 

2 

A electrode area, ft 

B viscosity coefficient defined by equation [6-8], £/gmole 
c constant in equation [2-10], )J./gmole 

3 

C mol ar concentration in electrolyte, lb mole/ft 
d matrix thickness, mils 

D diffusion coefficient or molecular diffusivity, ft /hr 

E cell potential, V 

f compression factor in equation [5-1] 

F Faraday's constant, 96,487 C/eq 

g constant in equation [2-6] 

G parameter defined by equation [4-33] 

h mass-transfer coefficient, ft/hr 

3 

H Henry's Law constant of CO 2 in electrolyte, mm ft /lb. mole 

AH|^ heat of formation of water vapor, cal/gmole 

2 

I current density. A/ ft 

J mechanical heat equivalent, 4.184 J/cal 

3 

k reaction rate constant, ft /lb. mole hr 

3 

K equilibrium constant, ft / lb. mole 
Si constant in equation [2-6], VQmole 

O 

L ionic radius, A 

m molality, gmole/kg 

M molecular weight of solvent 



N molar flux, lb mole/ft^ day 

p constant in equation [2-17] 

P partial pressure, mm Hg 

q constant in equation [2-1 7J 

Q rate of heat generation, BTU/hr 

3 

r local rate of reaction, lb mole/ft hr 

3 

R universal gas constant, 555 mm ft /lb mole °R 
s partial molar entropy of hydration, gibbs/gmole 

S gas-electrolyte area per unit of electrode geometric area, 

2 2 
ft^/ft"^ 

T cell temperature, °K 

2 

u electrical mobility, lb mole ft /J sec 

V average molar bulk velocity of bulk solution, ft/sec 

3 

V volumetric gas flowrate, ft /hr 

_ 3 

V partial molar volume of ion, cm /gmole 

X coordinate dimension between electrodes, ft 
X electric field, V/cm 

z ionic charge 

Z mole fraction 

Greek letters 

a matrix labyrinth factor 

3 empirical parameter, equation [5-1] 

V activity coefficient 

r reduced activity coefficient, defined by equation [6-10] 
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A difference between CO 2 partial pressure in the gas phase and 
at gas-electrolyte interface, mm Hg. 

V specific interaction coefficients in equation [A-4], £/gmole 
6^ smaller side of rectangular air channel, ft 

6^ larger side of rectangular air channel, ft 
6^ length of air channel, ft 

e matrix porosity 

Ti viscosity of solution, centi poise 

X ionic equivalent conductance, mhos cm/eq 

p ionic strength of solution, gmole/il 

V number of ions in solution per one molecule of electrolyte 

C hydration number 

T matrix tortuisity 

(J) osmotic coefficient 

$ electrical potential, V 

\p function defined by equation [4-40] 

X ratio between the apparent molar volume of electrolyte and 
molar volume of water 

Superscripts 

t gas-electrolyte interface 

*■ equilibrium value 

backward reaction 
° infinite dilution 

* characteristic value 



dimensionless 


Subscripts 

a anode 
c cathode 
el electrostatic 
G gas 
h hydrati on 
I gas inlet 
j ionic species 
L liquid 
m multicomponent 
so salting out 
St structural 
T total 

0 gas outlet 

1 HCO" 

2 OH" 

3 Cof 

4 NMe^ 

5 CO 2 

6 H2O 

i reaction [2-1] 

ii reaction [2-2] 

iii reaction [2-3] 



iv reaction [2-4] 
+ cati on 
anion 


t activated complex in equation [6-23] 
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APPENDIX A 


DETAILED CALCULATIONS 
Activity Coefficients 

Table 1 contains the activity coefficients of tetramethyl'- 
amtonium halides used in interpolation in the concentration 
range from 0.1m to 1.5m (the iodide forms a saturated solution 
at 0.2715m). 


TABLE 1 

ACTIVITY COEFFICIENTS OF TETRAMETHYLAMMONIUM HALIDES AT 25°C 


m 

NMe^F 

NMe^Cl 

NMe^Br 

NMe^I 

0.1 

0.795 

0.746 

0.7203 

0.6985 

0.2 

0.776 

0.683 

0.6450 

0.6136 

0.3 

0.778 

0.644 

0.5976 


0.4 

0.787 

0.617 

0.5637 


0.5 

0.802 

0.597 

0.5375 


0.6 

0.823 

0.580 

0.5160 


0.7 

0.846 

0.566 

0.4974 


0.8 

0.873 

0.556 

0.4827 


0.9 

0.902 

0.546 

0.4688 


1.0 

0.931 

0.539 

0.4582 


1.2 

0.998 

0.527 



1.4 

1 .077 

0.520 



1.5 



0.4432 
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The electrostatic, hydration, and mutual salting-out terms for 
these halides were calculated using equations [6-4], [6-5], and 
the modified forms of equations [6-6] and [6-7], The hydration 
numbers used were the same as suggested by Glueckauf (1955), 0.9 
for Cl , Br , and I , 1.8 for F~, and 0.0 for NMe^^. The partial 
molal volumes tabulated by Mi Hero (1971) were also used, Table 
2 shows these calculated terms and the term log obtained by 
subtraction from log 

The term (B. - 0.0025 V.) used in interpolation was calculated 
J J 

from the values of B. and V. shown in Table 3 for the ions in 

J J 

question. These values of (-log y^^) were plotted against the 
term (B^ - 0.0025 Vj), for the halides and then used to calculate 
the same terms for (NMe^)2C02> NMe^HCO^, and NMe^OH as shown in 
Fig. 23. In the case of NMe^OH the activity coefficients ex- 
ceeded unity. Thus it was inappropriate to use the activities 
as a measure of ion density. Instead, the Debye-Huckel term was 
calculated using the ionic strength as a measure of the ion 
density with the ion size parameters calculated with the follow- 
ing equation of Stokes and Robinson (1957): 

L. = 0.613 (V. + 18 [A-1] 

u J J 


and 
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TABLE 2 

TERMS IN EQUATION [6-3] FOR THE NMe^-HALIDES 


m 

log Tg, 

log 

log yr-r. 
^ 'so 

'“9 'fst 

NMe^F 

0.1 

-0.0991 

0.0138 

0.0054 

-0.0197 

0.2 

-0.1233 

0.0274 

0.0105 

-0.0247 

0.3 

-0.1390 

0.0408 

0.0159 

-0.0267 

0.4 

-0.1509 

0.0540 

0.0214 

-0.0285 

0.5 

-0.1605 

0.0670 

0.0272 

-0.0294 

0.6 

-0.1689 

0.0798 

0.0336 

-0.0291 

0.7 

-0.1761 

0.0925 

0.0402 

-0.0292 

0.8 

-0.1826 

0.1049 

0.0475 

-0.0288 

0.9 

-0.1885 

0.1171 

0.0552 

-0.0286 

1.0 

-0.1938 

0.1292 

0.0633 

-0.0298 

1.2 

-0.2033 

0.1528 

0.0814 

-0.0318 

1.4 

-0.2117 

0.1758 

0.1024 

-0.0345 

NMe^Cl 

0.1 

-0.0943 

0.0089 

0.0062 

-0.0481 

0.2 

-0.1145 

0.0177 

0.0113 

-0.0801 

0.3 

-0.1268 

0.0262 

0.0160 

-0.1065 

0.4 

-0.1357 

0.0345 

0.0204 

-0.1289 

0.5 

-0.1483 

0.0427 

0.0247 

-0.1431 

0.6 

-0.1540 

0.0507 

0.0288 

-0.1621 

0.7 

-0.1530 

0.0585 

0.0328 

-0.1855 

0.8 

-0.1572 

0.0661 

0.0368 

-0.2006 

0.9 

-0.1608 

0.0736 

0.0407 

-0.2163 

1.0 

-0.1642 

0.0809 

0.0446 

-0.2297 

1.2 

-0.1699 

0.0950 

0.0542 

-0.2557 

1.4 

-0.1749 

0.1085 

0.0603 

-0.2779 
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TABLE 2 - Continued 



m 

log 

log 

log V 3 „ 

log Y 5 t 

NMe^Br 

0.1 

-0.0924 

0.0104 

0.0063 

-0.0668 

0.2 

-0.1114 

0.0206 

0.0113 

-0.1109 

0.3 

-0.1228 

0.0305 

0.0156 

-0.1469 

0.4 

-0.1308 

0.0402 

0.0197 

-0.1781 

0.5 

-0.1371 

0.0497 

0.0234 

-0.2056 

0.6 

-0.1421 

0.0589 

0.0270 

-0.2312 

0.7 

-0.1463 

0.0679 

0.0304 

-0.2553 

0.8 

-0.1500 

0.0767 

0.0337 

-0.2767 

0.9 

-0.1531 

0.0854 

0.0368 

-0.2981 

1.0 

-0.1560 

0.0938 

0.0400 

-0.3167 

1.5 

-0.1590 

0.1332 

0.0580 

-0.3856 

NMe^I 

0.1 

-0.0903 

0.0127 

0.0063 

-0.0845 

0.2 

-0.1136 

0.0250 

0.0113 

-0.1403 



TABLE 3 




IMPORTANT PARAMETERS 

OF THE IONS 



Ion 

0 

B V, 

3 ,-l 

cm mol s. 

gibbs mol""* 

L,A C 

F" 

0.155 

4.2 

3.0 

1.36 1.8 

cr 

0.014 

23.2 

18.5 

1.81 0.9 

Br" 

- 0.021 

30.1 

24.6 

1.95 0.9 

r 

-0.060 

41.6 

31.4 

2.16 0.9 

C032- 

0.171 

6.5 

- 2.0 

1.99 

HC 03 " 

-0.040 

28.8 

28.0 

1.33 

OH" 

0.181. 

1.4 

2.8 

1.33 


0.099 

84.2 


3.47 0.0 



LOG 
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L = + L_ [A»2] 

The hydration numbers of (MMe^)2C02, NMe^HCO^, and MMe^OH 
were chosen so that the experimental osmotic coefficient data 
at the lowest concentrations agree with the value obtained by 
the equation (Robinson and Stokes, 1959) 

m 

4> = 1 + (1/m) m d In Y [A-3] 

0 

The values obtained for these electrolytes are 1.5, 0, and 6.4, 
respectively. The small hydration number obtained for NMe^HCO^ 
may be attributed to a higher degree of association than that ex- 
pected by simple structural arguments. The calculated activity 
coefficient terms for the three electrolytes are shown in Table 4 
while the experimental and calculated osmotic coefficients and 
the derived activity coefficients are shown in Table 5. 

For concentrations below 0.1m, the activity coefficients can 
be calculated by the Guggenheim equation 

log Y^ = -0.5809 |z^z_| + V.^.y [A-4] 

where V.. is the specific interaction coefficient between ions i 

* O 
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TABLE 4 

OSMOTIC AND ACTIVITY COEFFICIENTS OF 



(NMe^jjCOj, 

NMe^HC03, AND 

NMe^OH 

AT 25°C 


m 

d) 

^exp 

*^0310 

Y 

m 

(f) 

^exp 

'^calc 


(NM 


0.1 


0.726 

0.578 

1.202 

1.377 

0.975 


0.2 


0.823 

0.548 

1.603 

1.857 

1 .88 


0.3 


0.877 

0.545 

2.060 

2.417 

4.38 


0.4 


0.924 

0.555 

2.584 

2.974 

11.3 


0.5 


0.967 

0.573 

3.205 

3.434 

28.6- 


0.534 

0.981 

0.981 

0.579 

3.934 

4.027 

90.5 


0.848 

1.108 


0.685 

4.808 

4.326 

230 



le^HC 03 

0.1 


0.811 

0.685 

1.307 

0.604 

0.614 0.286 


0.2 


0.792 

0.589 

1.852 

0.801 

0.313 


0.3 


0.767 

0.526 

2.469 

1.012 

0.375 


0.4 


0.743 

0.479 

3.175 

1 .158 

0.443 


0.5 


0.720 

0.441 

3.989 

1 .323 

0.552 


0.6 


0.700 

0.410 

4.938 

1 .405 

0.648 


0.7 


0.680 

0.383 

6.061 

1 .518 

0.797 


G.8 


0.662 

0.360 

7.407 

1.547 

0.913 


0.9 


0.645 

0.340 

9.054 

1 .619 

1.10 


1.0 


0.629 

0.322 

11.111 

1 .688 

1 .35 


1.2 


0.614 

0.296 






IMe^OH 

0.1 


0.926 

0.863 

1.221 

1.505 

2.21 


0.2 


1 .008 

0.916 

1.939 

1 .702 

3.55 


0.3 


1.113 

1.042 

2.747 

1.982 

6.30 
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TABLE 4 - Continued 


m 

^exp ‘^calc 

T 

1 

m 

(j) 

exp 

^calc ^ 

NM 

le^OH 

0.4 

1.183 

1.166 

3.663 

2.366 

13.0 

0.5 

1.252 

1 . 31 2 

4.710 

2.722 

27.3 

0.5783 

1.320 1.320 

1.463 

5.917 

3.198 

68.7 


TABLE 5 


THE TERM LOG FOR THE HYDROXIDES AT 25°C 



LiOH 

NaOH 

KOH 

CsOH 

NMe^OH 


0.6 A 

0.95 A 

1.33 A 

1.65 A 

3.47 A 

m \ 






0.1 

-0.0694 

-0.0416 

-0.0253 

-0.0091 

-0.0179 

0.2 

-0.1130 

-0.0724 

-0.0436 

-0.0270 

-0.0201 

0.3 

-0.1579 

-0.1038 

-0.0624 

-0.0463 

-0.0211 

0.4 

-0.2033 

-0.1366 

-0.0803 

-0.0638 

-0.0218 

0.5 

-0.2501 

-0.1712 

-0.0975 

-0.0816 

-0.0223 


and j. From the calculated activity coefficients of (NMe^)2C02> 

NMe^HCOg, and MMe^OH the coefficients were found to be 0.439, 

_1 

-0.420, and 0.583 mole , respectively. For comparison the V^.j 
value for NMe^OH is 0.5 +0.1 mole"^ as obtained by cell 
measurements (Prue et al., 1971). 

According to the structural hydration model (Desnoyers 
et al., 1969), the activity coefficients of alkali and tetra- 
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methyl ammonium electrolytes with the same anion vary in a sys- 
tematic manner as the cation size increases. Table 5 shows the 
calculated log term for the hydroxides which exhibit that 
systematic change. 

To test the accuracy of this method, the activity coefficients 

of Na 2 C 02 were calculated by adding the electrostatic hydration, 

mutual salting-out and structural terms, calculated as described 

2 - 

earlier. Since the CO^ and OH ions are very similar in their 
structural effects (Fig. 23), it is reasonable to use the values 
of log y^^ of MaOH in Table 5 for Na2C02. The calculated ac- 
tivity coefficients are compared with the experimental values in 
Table 6. These values agree within + 5 %. 

TABLE 6 

ACTIVITY COEFFICIENTS OF Ma2C03 AT 25°C 


m 

\alc. 

*^exp. 

0.1 

0.489 

0.466 

0.2 

0.406 

0.394 

0.3 

0.357 

0.356 

0.4 

0.319 

0.332 


0.5 


0.288 


0.313 


no 


Characteristic Parameters 


Table 7 contains the characteristic parameters used in calculations. 

Table 7 

CHARACTERISTIC PARAMETERS 


C*, lb mol/Jl 

0.4495 

D*, ft^/hr 

1 X 10 

d*, mils 

30 


7.7 


Equilibrium and Henry's Law Constants 


The concentrations of the different spaces used in calculating 
the average equilibrium constants are shown in Table 8. These con- 
centrations are approximately the same as obtained at the base con- 
ditions (P^g = 3mm Hg, I = 18A/ft , and inlet air dewpoint = 60°F). 
Table 8 also contains the factors which should be multiplied by K? 


Table 8 

AVERAGE CONCENTRATIONS AND EQUILIBRIUM CONSTANTS 




^2 S 

g mole/£ 

K*/K° 
1 1 

K. ./K°. 

TV n 

H/H° 

Cathode Zone 

0.003 

0.5 

1.80 

33 

8.1 

1.686 

Matrix Zone 

- 

- 

- 

20 

4.5 

- 

Anode Zone 

1.8 

10"^ 

0.018 

7.1 

0.85 

1.190 
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Diffusion Coefficients and Mobilities 


Table 9 contains the values of X° and D° for the different ions 
at infinite dilution. Table 10 shows the calculated diffusion co- 
efficients and Table 11 shows the calculated equivalent conductivity. 


table: 9 

EQUIVALENT CONDUCTANCES AND DIFFUSION COEFFICIENTS 
AT INFINITE DILUTION AND 25°C 



Ion 

2 -1 -1 
A°,cm Q equiv 

D°xlO^,cm^/sec 



C 05 - 

69.3 


0.922 



HCO- 

44.5 


1.185 



OH" 

198.3 


5.278 



NMe^ 

44.9 


1.196 





TABLE 

10 




DIFFUSION 

COEFFICIENTS AT 25 °C 



m 

C 

1 + diLny 
dJlnm 

0 /Tl° 

DxlO^,cm^/sec 

D"xl0^,cm^/sec 




(NMe^ 

^2^*^3 


0.1 

0.0982 

0.913 

1.12 

0.887 

0.752 

0.2 

0.193 

0.960 

1.21 

0.863 

0.731 

0.5 

0.461 

1.083 

1.46 

0.807 

0.731 

1.0 

0.846 

1.905 

1.91 

1.085 

0.919 

2.0 

1.482 

4.772 

3.15 

1 .648 

1.397 

3.0 

1.969 

5.718 

7.47 

0.833 

0.706 

4.0 

2.360 

6.552 

12.1 

0.589 

0.499 
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TABLE 10 - Continued 


m 

C 

1 + 

dilnm 

n/n° 

DxlO^, 

2 

cm /sec 

- 5 2 

D xlO ,cm /sec 

NMe^HCOg 

0.1 

0.100 

0.820 

1.07 

0.912 


0.908 

0.2 

0.195 

0.752 

1.08 

0.829 


0.826 

0.5 

0.472 

0.610 

1.18 

0.615 


0.612 

1.0 

0.895 

0.465 

1.37 

0.404 


0.402 

2.0 

1.638 

1.649 

1.69 

1.161 


1.156 

3.0 

2.229 

1.730 

2.25 

0.915 


0.911 

4.0 

2.740 

1.787 

3.09 

0.688 


0.685 

5.0 

3.176 

1.850 

4.04 

0.545 


0.543 

10.0 

4.660 

2.100 

11.8 

0.212 


0.211 

NMe^OH 

0.1 

0.0991 

1.035 

1.03 

1.958 


5.302 

CM 

O 

0.196 

1.240 

1.07 

2.270 


6.147 

0.5 

0.480 

1.571 

1.24 

2.469 


6.686 

1.0 

0.917 

1.685 

1.37 

2.397 


6.491 

2.0 

1.711 

2.510 

2.08 

2.352 


6.369 

3.0 

2.357 

3.333 

2.78 

2.337 


5.329 


TABLE 11 

EQUIVALENT CONDUCTANCES AT 25° C 


c 

^2^^2 

A 

NMe^HCO 

C 

3 

A 

NMe^OH 

C 

A 

0.49 

55.8 

0.35 

63.7 

0.55 

193 

0.74 

46.4 

0.72 

53.5 

1.11 

166 

0.99 

37.9 

1.12 

44.6 

1.67 

137 

1.26 

29.0 

1,50 

38.7 

2.23 

112 

1.52 

22.2 

1.90 

33.3 

2.81 

83.1 

1.69 

17.3 

2.32 

28.8 

3.39 

66.7 

2.06 

11.3 

2.72 

25.4 

3.98 

47.2 
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TABLE 11 - Continued 



The diffusion coefficients and mobilities (as v^ell as the in- 
crease in mobilities for a ten degree Fahrenheit increase) of each 
anion were fit to polynomials in concentration of sixth degree. The 
coefficients of these polynomials are given in Table 11. The vis- 
cosities were approximated by a single polynomial of fourth degree 
such that at any point the viscosity is given by 


n = n(C-, + + 2C3) [A-5] 

The fitting equation is 

n = 0.8937-1.00938 C + 2.37864 - 1.13294 + 0.17625 

CA-6] 

The change of viscosity with temperature is assumed to be proportional 
to that of water. 


Cathode and Anode Active Areas 

Instead of calculating the magnitudes of the cathode active 
areas, the dimensionless parameters (?l'g ) was calculated. This 

parameter is proportional to C^ which should also give the same value 



TABLE 12 


COEFFICIENTS OF THE POLYNOMIALS USED IN CALCULATING 
THE IONIC DIFFUSION COEFFICIENTS AND MOBILITIES 


Parameter 

C° 

C^ 

c2 

C^ 


C5 

0 

cn 

COj^- 

DxlO^, cm^/sec 

0.90718 

-1.10333 

-0.29669 

4.80487 

-4.31020 

1.25519 

-0.10134 

u F^ 

34.395 

- 47.222 

27.566 

5.098 

-14.045 

5.796 

-0.754 

Au F^ 

2.437 

0.380 

-2.072 

0.199 

0.890 

-0.435 

0.060 

HCO 3 - 

5 2 

DxlO , cm /sec 

1.22224 

-3.13108 

4.17542 

-1 .85812 

0.23239 

0.02796 

-0.00587 

u F^ 

44.392 

-45.917 

40.410 

CM 

CO 

1 — 
CM 
1 

6.027 

-0.875 

0.051 

Au F^ 

2.415 

1.739 

-2.139 

0.712 

-0.069 

-0.007 

0.001 

OH" 

5 2 

DxlO , cm /sec 

5.13441 

5.60765 

-6.31471 

1.45324 

1 .28160 

-0.740G1 

0.10685 

up2 

198.110 

-100.195 

59.105 

-17.696 

-3.231 

2.461 

-0.309 

Au F^ 

14.653 

-4.307 

4.263 

- 2.789 

0.285 

0.155 

-0.028 
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of as calculated from the differential equations vjhich represent 
the mass transfer processes in the matrix zone. By trial and error 
the correct value of was obtained in each case from which the 

Q 

value of ) was calculated. As this parameter is pro- 

'll c ^c 

sumably constant at a given current density and total ionic strength, 
an average value was calculated for each current density and air dew 
point. For example, the average value of this dimensionless parameter 
at 18ASF and 60°F for an inlet dew point in the analytical cell is 
2.820 X lo”. The relative change of this parameter with current den- 
sity is the same as that of cathode active area (Fig. 10). 


To calculate the value of we can use the experimental value 
of CO 2 flux obtained for the full-scale cell at a current density of 
18A/ft^, Ppr, of 3mm Hg, and V of 600 cm^/min. Equation [4-32] 

UU2 a 

gives a value of 0.004775 gmole/J?, for Cr . Previous measurement of 
the pH in the NMe^HCOg solutions at a concentration of 1.8 gmole/il 
under constant pressure of CO 2 was 8.21 g (Abdel-Salem, 1974). As- 
suming that the pH of the solution measures the hydrogen ion ac- 


tivity and the activities of hydroxyl ions are the same as those of 


NMe^OH at this low concentration, then the ratio 0^/02^^ at the anode 
is 0.G06508 gmole/il. From equation [4-42] the value of the dimen- 
sionless parameter ['^^^ 5(^2 ^ in that case is 1440 units. In- 

creasing the value of this parameter by several times has only little 


effect on the calculated CO 2 removal rate. Therefore this value is 
used throughout. 



APPENDIX B 


COMPUTER PROGRAM 

Calculation of the heat, water, and CO 2 transport in the cell are 
performed by the attached computer program. The program consists of 
a main program and several subroutines. The calculations for any set 
of conditions proceed as follows: 

(1) The subroutine ENERGY calculates the temperature of each subcell 
and the catholyte concentration. The subroutine LAGRAN is used in 
interpolation from the table of vapor pressure data for (NMe^) 2 C 02 
solutions. (The interpolation is based on the method of Lagrange. 

The calculated temperature, catholyte concentration, and outlet air 
dew point depend on the current density, cell voltage, and air flow 
rate together with the inlet air temperature and dew point. These 
parameters must be specified for each case. 

(Z) An initial guess of the COg flux is made. This is used to cal- 
culate the CO^ concentration in the catholyte which is then used to 
calculate the OH*" concentration in the catholyte. From this informa- 
tion the concentrations of CO^ and HCO^ are obtained at = -1/2. 

(3) With the guessed value of the CO 2 flux the differential equations 
[4-34] and [4-35] are integrated using subroutine CRANE. The deriva- 
tives are evaluated in each step by the subroutine DERIV. The sub- 
routine DIN calculates the ionic mobilities and diffusion coefficients. 

The subroutine DATAI is used to calculate the equilibrium constants at 

2 - - . 

the cell temperature. From the calculated CO^ and HCO 2 concentrations 
at the anode the anode CO 2 flux is calculated and compared with the 
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guessed cathode flux. 

(4) The boundary condition at the anode is satisfied when the cal- 
culated anode flux is within 5% of the guessed value at the cathode. 
This results in an error of no more than 0.01% in the calculated CO2 
flux. Convergence is obtained more rapidly if the initial guess and 
upper and lower limits for the flux are reasonable. For any set of 
operating conditions the initial guess can be obtained approximately 
from Figs. 11 and 16. This reduces the computer time considerably. 

(5) When the boundary conditions are satisfied the concentration 
distribution across the matrix is obtained for all the ionic species. 
If the concentration exceeds the solubility limit at any point a 
warning is issued to that effect. It is sufficient to check the 
catholyte and anolyte and not the entire matrix. 

The program is written in FORTRAN G1 language. The calculations 
were carried out on IBM 360 machine. The attached listing shows cal- 
culations for four cases using different operating conditions. In 
one case precipitation is indicated when a high current density (30A/ 
ft ) is combined with a low air relative humidity. 



C THIS PRGCRAM PERFORMS CALCULATIONS FOR THE MATHEMATICAL MODEL WHICH 
C SIMULATES HEAT, WATER ANC CARBON DIOXIDE TRANSFER IN ELECTROCHEMICAL 
C C02 CONCENTRATOR! ONE MAN AIR REVITALIZATION SVSTEM J 

C THESE CELLS USE TMAC ELECTROLYTE 
C 

C THE INPUT PARAMETERS REQUIRED BY THE PROGRAM ARE AS FOLLOWS 
C yiNPUTy IS A LIST OF PARAMETERS RELATED TO THE CQ2 TRANSFER 
C IC- IDENTIFICATION NUMBER , 1 FOR THE ANALYTICAL CELL AND 2 FCR THE 

C FULL-SCALE CELL 

C PC02A= ANODE SIDE PC02 ♦ INLET VALUE USED HERE* MM HG 
C BLFALL= LOWER BOUND FOR Tl OF ALL SUBCELLS 
C ByFALL=UPPER BOMND FOR TI OF ALL SUBCELLS 
C NTSC= NO. OF TOTAL SU8CELLS 

C CURENT= CURRENT DENSITY OF THE CONCENTRATOR 

C PC02IN= CATHODE (AIR) SIDE INLET PC02, MM HG 

C TI= TRANSFER INDEX* LB. C02 TRANSFERRED PER LS. €2 CONSUME* HERE A GUESSED 
C VALUE FOR THE FIRST SU8CELL IS USED FOR INPUT 

C TRANSFER INDEX OF 2.75 CORRESPCNDS TO 100% CG2 CURRENT EFFICIENCY 
C X0= DISTANCE BETWEEN ELECTRODES* MILS. 

C CFM= AIR FLOW RATE, CUBIC FT. PER MINUTE PER CELL ( FULL-SCALE BASIS) 

C = (AIR VELOCITY, FT/SEC) X I.S2 

C TIN= INLET AIR TEMP. 

C TDEWIN= INLET AIR DEW POINT TEMP. 

C ALL TEMPERATURES ARE IN DEGREE F 
C VOLT= VOLTAGE OF THE CELL, VOLTS 

C VHV=HYDROGEN FLOW RATE, FT3/MI N. CELL ( F ULL- SCALE BASIS) 

C 

C DESCRIPTIONS OF IMPORTANT INTERMEDIATE PARAMETERS 

C XP(N)= X-COOROINATE OF THE N-TH NODE AT WHERE THE VALUES OF CCNGENTRAT IONS 
C ARE SAVED 

C DXP= DISTANCE BETWEEN NODES 

C NSC= THE NO. OF THE SUBCELL BEING SIMULATED 
C FLUX= FLUX OF SUBCELL CC2 TRANSFER 
C BUFLUX= UPPER BOUND FOR FLUX 

C BLFLUX= LOWER BOUND FOR FLUX 



>d § 
8 ® 



C H= INITIAL INTEGRATION STEP FOR CRANE 
C Y( 1 )= FIRST DEPENDENT VARIABLE IN CRANE, Cl 
C Yi2)= SECOND DEPENDENT VARIABLE IN CRANE, C3 

C F{1)= DERIVATIVE OF YU) 

C F(2)= derivative OF YC2) 

c 

C OUTPUT OF THE MAIN PROGRAM - WHAT THE MAIN PROGRAM CALCULATES 
C Cl= BICARBONATE ION CONCENTRATION 
C C2= hYDRCXYL IGN CONCENTRATION 

C C3= CARBONATE ION CONCENTRATION 

C C4=TMA ION CONCENTRATION 

C PC02A-OUTLET PC02 OF ANGCE SIDE 

C C02 REMOVAL RATE ,C02 REMOVAL EFF I CIENC Y, C02 TRANSFER INDEX 
C 

c 

DIMENSION XP(i01),ClCl0l)*C3(101),Fl(l01) ,F2(101) ,CSCATH(40) , 
,TELECI40),XPP<101 ) 

DIMENSION U1RI4) ,01R(5) ,D<5) 

LOGICAL DONT, ENDONX,OKPC 

REALMS X , Y , F , PREC . X H AT , H , FLUX , CU R ENT , FL U X 1 . BL F ALL • BL FL UX , 8 UFLUX , 

4BUF, TFLUX,RFLUX-,BFLUX 

CCMMON yCRANEC/X, Y{ 40 ) ,FI 40 ),PR£C , XHAT »H,N»NB 140) ,MIN,I NIT, DGNT, 
*ENDCNX,OKPC 

COMMON EQK2M,EQK1M, EQK2C ,EOK2A,EQKlA ,EQK1 C ,CKDS , AKD ,HC ,H A , VELC, 

♦VELA, FENRyC,HENRYA,HAK ,FKC,CFM, IC 
COMMON yBLKl/FACMAT ,TVISF ,TEMPD 

COMMON /BLK2/CURENT,FLUX,C4,DTA2,Vi,V2, V3, V4 
COMMON yBLK3yVHV 

N A M EL IST/ 1 NP U1 / 1C , P C O 2A , BLFA L L . B UFA LL , N 7 SC , C U HE NT , F C02 IN,TI,XC, 

♦ CFM, VHV,TIN,TDEWIN,VCLT 

C THE FOLLOWING IS NECESSARY TO AVOIDE THE INTERRUPTION OF CRANE SUBROUTINE 
C WHEN SOME VALUES EXCEED THE LIMITS OF MACHINE CCMPLTATIGN (THESE VALUES 
C DO NOT AFFECT THE RESULTS BUT ARRIS E WHEN SOME TESTS ON THE ERRORS ARE 
C MADE) 

CALL ERRSET (207,300,-1) 





SB 

S lEo- 


CALL ERRSET (2D8,30O»-lJ 
GALL ERRSET (2Gg.300*-l> 

1 READ{5, INFUT* END=99gg } 

IF(IC .EQ. 1)WRITE(6*1111) 

1111 FORMATt *1»^ lOX. • ANALYTICAL CELL AT THE FOLLOWING CONDITIONS J*> 

IF( IC -HO.. 2itfRITE(6. 1112) 

1112 FORMAT! • 1 * ,1 OX. • FULL SCALE CELL AT THE FOLLOWING CONDITIONS 2«) 

WRITE(6, 1113)PC02IN.CURENT, XD.CFM .TIN.TDEWXN 

1113 FORMAT! * 0* »10X.39<1H’^‘)/1<3X»* INLET AIR PCQ2 = • , 3X, F6 . 2. 3X- • MM HG 

1 '/lOX.* CURRENT DENSITY = » .3X .F6- 2 ,3X * * A /FT2 ’XIOX, 

2 • MATRIX THICKNESS=« .4X.F6.2.3X. » MILS */lOX» 

3 » AIR FLOW RATE = • .3X . F6 . 2.3X . • FT3/MIN */10X, 

4 • AIR TEMPERATURE = • , 3X- F6- 2* 3X , • DEG,. F '/lOX- 

5 • AIR DEW POINT = • ,3X . F6 .2, 3X, • DEG . F •) 

WRITE ( 6- 1 1 14) 

1114 FORMAT {» 0^ t lOX, » OTHER CONDITIONS : ») 

WRI TE I 6 . 1 1 1 5 ) VOLT , VH V . PC C2 A 

1115 FORMAT! *0*- lOX. 'CELL VOLTAGE = • , 3X . F6 . 2 .3 X. » VOLT *yiOX. 

2 » HYDROGEN RATE = • *3X ,F6 .2. 3X. • FT3/MIN •✓lOX. 

3 • ANODE PC02IN = * ,3X ,F6 . 2 , 3X . • MM HG * /I OX ,40 ( IH* ) ) 

CALL ENERGY < NTSC, VOLT. CURENT,TIN , TDEW IN ,CFM ,TELEC ,CSCATH) 

C THE FOLLOWING STATEMENTS SPECIFY THE MATRIX FACTOR AND THE FACTOR WHICH 
C SHOULD BE MULTIPLIED BY THE AIR - CATHCLYTE AREA TO ACCOUNT FOR THE EFFECT 
C OF THE CURRENT DENSITY 
C FACMAT= THE MATRIX LABYRINTH FACTOR 
C ACAR =AIR-CATH0LYTE AREA FACTOR 
FACMAT=.34 

IF! IC .EQ. 2) GO TO 111 

ACAfi=3.70-3 .35/ (1 .+ (1 ,/ .G6/CURENT ) ♦’«' 8. 

GO TO 112 

111 ACAR=1 .37-1 .24/! 1 .+ ! l-/.0 6/CURENT)-.«c=«t 6. )*=(«. 5 

112 CONTINUE 
XD=XD/30.. 

DXP=XD/S0. 

XP! l)=-XO/2. 



DO 3 KP=2,51 
XP<KP)=XP (KP- 1 )+DXP 
3 CONTINUE 

C CURENT IS NONDIMENSIGNLIZEO AND THE GUESSED FIRST SUBCELL TI IS 
C CCNVERTED INTO FLUX 

CURENT=0 . 04576>!^CUR£NT 

FLUX=Tiy'5,5*CURENT 

BLFLUX=BLFALL*CURENT/5 .5 

B UFL UX=BUFALL+CURENT/£ • S 

N=2 

NSC=1 

R=0. 

LL=1 

C THE CATHQLYTE TMAC CONCENTRAT I GNS OBTAINED FROM THE SUBROUTINE ENERGY ARE 

C NQN-DINENS IONLI2ED 

4 C4C=CSCATH( NSC1/7.2 

CALL DATA I ( TELEC { NSC) ) 

VELC=VELC*NTSC 

VELA=VELA4NTSC 

C PC02 IS NON-OINENSIONLIZEO 

PC02I N=PC02 IN4HENRY C 
PC02A=PC02AYHENRYA 
WRITE( 6,401 ) SUFLUX.BLFLUX 

401 FORMAT! *0* ,10X, 'UPPER FLUX LIMIT= • , G12 ,62 1 OX . • LOWER FLUX L1MIT= 

l*iG12.6//) 

IT=0 

H=XD/1.0E5 

5 C5C=PCQ21N41 .E6-FLUX/VELC/( l,-EXPt-HC/VELC) )^1.E6 
C5A=PCa2A+FLUX/VELA/ { 1 EXP f-HA/V EL A ) ) 

C CORRECTION TO INCLUDE THE EFFECT OF IONIC STRENGTH ON REACTION RATE CONSTANT 
C C02 DIFFUSIVITY* AND HENRY'S LAW CONSTANT 

CKD=2.820T:ACAR=>I'SQRT ( (5, 1/ VI SCCS ( C4C+7 ,2 ) TTV ISF ) .474c 

1 10.?^4{,Oe3+(3.6-7.2=^'CAC)) )*CKDS 

C2 C= FLUX/ CK 0/ C5 C 
C3C=(C4C-C2C ) /{2.+EQK2C/C2C) 



c 


C 1C=E QK 2C y C 2C 3C 

INITIALIZE Y(l) AND Y(2) FOR DER I V AND CRANE 
7 Y( 1)=EQK2M%C3C/C2C 

Y(2)=C3C 

IF(C1C .LT. 0*) GO TO 24 
K=0 

X=~XDy2. 

XHAT=XD/2 • 

DONT^.FALSE^ 

MIN=3 
PREC=5 . 

FLUX INFLUX 
WRITE (6 *700) FLUX 

700 FORMAT! •-•,10X*47(lH-)yi0X*' THE INITIAL FLUX AT THE CATHODE IS •* 
IG12*6/10X,48( IH-)) 

H=XO/l .0E5 
ENDQNX=.TRUE. 

INIT— 0 
CALL DERI V 
XX=-*5 

WRITE(6*701 )XX»Y{ 1),Y(2) ,C4 

701 F0RMAT(//10X.» X =• , G12 . 6 ,3X* • C 1 C =» . Gl 2*6 , 3X , • C3 C =• , Gl2 .6, 3X. »C4 
1C =•* G12.6) 

XN=X 

C1(I)=C1C 
C3 1 1 ) = C3C 
F1(1)=F(1) 

F2( 1 )=F(2) 

ClN=CiC 
C3N=C3C 
F1N=F(1) 

F2N=F(2) 

KP 1=2 
CONTINUE 

IF{Y( 1 ) > 88*88 ,803 


8 





' 1 ^ 


88 LL^LL+l 

IF(Li_ *NE. 3JGD TO 88 1 
FLUX=TI=!fCURENT*l ,1/5 ,5 
GO TO 5 

881 IF(LL ,.GT . 7)GO TO 4 l 
GO TO 24 

803 IF( Y{ 2) 189. 89, 805 

89 LL=LL+1 

IF(LL ,NE. 31 GO TO 882 
FLUX=8LFALL^CUReNT/5 . 5 
GO TO 5 

882 IF(LL .GT, 7.) GO TO 4 1 
GC TO 26 

805 CONTINUE 

CALL CRANE 

9 XO=XN 
C10=C1N 
C30=C3N 
F10=F1N 
F20=F2N 
XN=X 
C1N=Y{1) 

C3N=Y(21 

FIN=F(1J 

F2N=F<2) 

IF( X.LT.XPIKPl) .OR, KPl .EQ. 51) GG TO 13 
KP2=KP1+1 

10 IF( X.LT.XP(KP2) ) GO TC 11 
IF(KP2 ,EQ. 51) GO TO 11 
KP2=KP2+1 

GO TO 1 0 

11 KP2=KP2-1 

DO 12 KP=KP1,KP2 
SL0P=1 XPXKP)-XO) /(XN-XO) 

EVALUATING AND SAVING VALUES OF Cl AND C3 AT THE NODES 


C 



C3< KP.) = C3 0+( C3N-C30 ) =*SLOF 
Cl( KP)=C10+{ClN-“ClO) *SLQF 
FI ( KP )=F1Q4 ( F IN-F 10 ) ^SLOP 
F2( KP)=F2C+(F2N-F20) 5«cSLaP 

12 CONTINUE 
KPI=KP2+1 

13 CONTINUE 
K=K + 1 

IFIMOOIK^IOOG ) .EQ« 0) WR1TE{6,14) K,X 
14 F0RNiAT(//10X, • INTEGRATION STEP NUMBER =» .G12 ,6 ,SX, *X =»,G12-6) 
IFIK .GT« 2000) GO TO 35 
IF((X+H/2.) .LT. XHAT) GO TO 8 
29 CONTINUE 
XP{51 )=X 
C1(51)=Y(1) 

C3{51)=Y<2) 

XX= . 5 

lfc‘RITE<6,2C.00)XX, Y( 1) ^YI 2) *C4*K 

2000 FORMAT (X/lOX* * X =* » G12 .6,3X . *C1 A = • ,G12. 6t 3X, '»C3A =• *G 12. 6.3X* 'C 
lA =• ,G12.6//10X.* THE INTEGRATION ROUTINE TOOK*, 14, » + 1 STEPS*) 
C1A=Y{ 1 ) 

C3A=Y{2) 

IF(CIA) 24,24,21 

21 IF(C3A) 26,26,22 

22 C2A=EQK2A^C3A/C1 A 

C TEST IF THE ANODE BOUNDARY CONDITION IS SATISFIED 
G=FLUX-1440,4 { EQK 1 A4=C1 A^ 1 .E6/ (C2A41 .£6 )-C5A.) 

CRIT = DABSIG/FLUX) 

IFICRIT ,ET . 0.05) GO TO 5 1 
WRITEI6,2200 ) G.CRIT 

2200 FORMATIXXlOX, * GUESSED CATHODE FLUX - CALC. ANODE FLUX IS*,3X,G12 
16 ,3X. • DIFF./FLUX =• , G12 .6 ) 

IF{G) 26,31,24 
8UFLUX=FLUX 
GO TO 27 


24 



26 


-a' 


BL-FLUX—FLUX 

C HALF INTERVAL TO DETERMINED NEWLY GUESSED FLUX 
27 FLUX=C BLFLUX+BUFLUX )/2. 

IF( DABSCFLUX-FLUXl ) .LT. . 000000 1 ) GO TO 36 
29 rT=iT + l 

IF( IT-50) 5,5,36 

31 continue 
WRITE(6,3101 ) 

3101 FORMAT! • 0* .10X.69 (1 H=) ) 

WRITE! 6,3100) 

3100 FORMAT (//18X , • X» , 13X , • Cl • , 13X , ■ C2« , 13X, *C3« , 1 3X, • C4 • ) 

WRI TE( 6,3102) 

3102 FORMAT! *0 • , lOX, 69! IH- ) ) 

DO 32 IP=1 ,51 

XPP! IP)=XP!IP)/XD 
£GK2 =EGK2M 

IF! IP .EQ, 1)EQK2=EQK2C 
IF!IP .EQ, £1 )EQK2=EQK2A 

C4=C 1 ! IP .) +2 . * C3 ! I P )+ E G K 2 +C3 ! I P )/ C 1! IP ) 

CR1=C 1! IP)*7. 2 

CR3=C3! 1P)*7,2 

CR2=EQK2 *C3(1P)/C1 tlP)¥7 

CR4=C4*7.2 

CALL 0IN!CR1,CR2,CR3,CR4,DIR,UIR) 

IF! IP ,EQ. 1 ) GO TO 311 
IF! MOD! IP, 10) ,NE. 1) GO TO 22 
311 WR1TE!6,3104 )XPP!IP) ,CF1 ,CR2,CR3,CR4 
210 4 F0RMATI//12X, 5!G12. 6,3X) ) 

IF! IF *EQ, 1 .AND, CR4 ,GE. 5 ,4 ) WR I TE ! 6 ,31 1 1 ) 

3111 FORMAT(//10X, * **WARN1NG’*^5«' PRECIPITATION EXPECTED AT CATHODE*) 
IF! IP ,EQ , 51 .AND, CR4 .GE . 7. 3) WRITE! 6,31 12 ) 

3112 FORMAT! //10X, * *^WARNING>f‘* PRECIPITATION EXPECTED AT ANODE* ) 

32 CONTINUE 
WRITE<6,3 103) 

3103 FORMAT! ■ 0* , 10X,69!1 H-) ) 



C R IS THE ACCUMULATED SUBCELL FLUX 

33 R=R+FLUX 

PC02A=(PC02A+FLUX/VELA)/HENRYA 
PC02IN=(PCD21 N-FLUX/VELC) /HENRYC 
WRITE{6 i330{) ) NSC^FLUX 

3300 FQRMAT(y/l OX* • FLUX FOR SUBCELL NUMBER*, 13, • IS *,012.61 
IFINSC .GE. NTSC) GO TO 35 
NSC=NSC+1 
BUF=FLUX^1.05 
IFINSC .GT . 2) GO TO 34 
RFLUX=FLUX 
FLUX=FLUX*0. 98 
GC TO 4 

34 TFLUX=RFLUX 
RFLUX=FLUX 

FLUX=(2 .*FLUX-TFLUX )=*>0 .9 8 
GO TO 4 

35 TI=5 .5+R/CURENTXFLOATtNlSC) 

CCCC=R5i£ .0431 /FLOAT ( NTSC) 

CCCD=CCCC*44. 

WRITE (6, 3502 )CCCC,CCCD 

3502 FORMAT! • 0* , 10X,» C02 REMOVAL RATE - *,G12.6,3X,* LB M0L/FT2.CAY */2 
17X, * = •,G12.e,3X,« LB /FT2.DAY*) 

WRITE (6,3500) TI,PC02A 

3500 format ( // lOX, » TRANSFER INDEX =• ,G1 2. 6,3X, • CD2 ANCDE PRESSURE =*,G 
112 . 6 ) 

EFF=TI*100./2.75 
WR1TE(6* 3S0DEFF 

3501 FORMAT! //lOX** C02 REMOVAL EFFICIENCY = * , 3X, F5 .2 ,3 X, * % *) 

GO TO 40 

36 CONTINUE 

C IF ANY OF THE FOLLOWING STATEMENTS IS PRINTED OUT TRY TO USE NEW INITIAL 
C VALUES For TI 

WRITE (6 ,3600 ) CRIT 

3600 FORMAK//* NON CONVERG ENCE • , 3X , • CRIT =*,614.8) 



GO TO 40 

39 Wf?ITE (6 ,3900 ) 

3900 FQRMAT(/y * TOO MANY R-K STEPS ,STQP C AL.CUi_AT I GNS * ) 

60 TO 40 

41 WKITE(6 ,4100 ) 

4100 FORMAT(//' INITIAL VALUES OUTSIDE RANGE, USE NEW GUESSED VALUES •/• 
1 IN THAT CASE THE USER IS ADVISED TO USE CLOSE VALUES OF TI ,BUFAL 
2L,BLFALL,G0QD INITIAL GUESS CAN BE OBTAINED FROM FINAL REPORT NASA 
39-1 2 £26 ( F £5,1576) • ) 

40 GO TO 1 
9999 STOP 
END 



SUBROUTINE DERIV 
DIMENSION D{S)^U(4) ,FRT<3) 

LOGICAL DONT.£NOONX,OKPC 

REALMS X , Y , F , PREC, X F AT , H . FLUX , CUR ENT 

COMMON XCRANECXX ,Y(40 3 ,F (4 0J . PREC ,X HAT , H-, N, NB (40 J , M IN. INIT ,DONT 
+ENDONX.OKPC 
COMMON SQK2.EQK1M 

COMMON /BLK2/CURENT,FLUX,C4.DTA2»V1 .V2.V3.V4 

CRl-Y (11+7 .2 

CR2=Y( 2>*7,2+EQK2/Y( 1 ) 

CR3=Y(2)+7.2 

CR4=CR1 +CR2HH2 +CR3 

C4=Y( 1 )+2»+Y( 2) «-Y(2) + EQK2/Y( 1 > 

CALL OINICRl ,CR2. CR3, CR4.D.UI 

DO 50 11=1,3 

FRT( II)=U(IIJ+0(4}/U( 4) 

CONTINUE 

W1=D(1 >+FRT(l)+(y(l)-eGK2+Y(2)yY(l ) ) /C4 
W2=2. +Y {2)+(l,-£QK2+Y(2)/Y( I ) ++2 ) +FRT ( 3 )/C4 
W3= { 2 . + Y ( 1 ) +EQK2 ) +F RT ( 1 ) /C4 
W4=D( 3)42 .+FRT( 3)+( l.-Y{ 1) /C4) 

Vl = 1»ll 4W2+D(5 ) +2. +EQK1W/EQK2+Y (1 )/Y (2) 

V2=W3+W4-D(S) +EOK IM/EOKE* ( Y{ 1 1/ Y( 2 ) )+*2 

V3=Wl+2,+V/2+FRT (2) +E0K2+Y(2 )/Y( 1 )/C4 + ( 1 .-EQ K2+Y( 2 ) /Y( U+ + 2) 
V3=V3-D(23+EQK2+Y(2)/Y( D++2 

V4=W3+2.+W4+0( 2) +EQK2/Y( 1)+FRT( 2) +EQK2+ ( 1 .-Y ( 1) /C41 XY( 1 ) 
DTA2=Vl3?cV4-V2+V3 

FLUXX=FLUX-CURENT*(Y{ 13+Y(2>-Y{1 »+U(1 ) /U(4 J-2+Y (2 ) +U ( 3 )XU ( 4 ) )/2 
1 /7.7 

CURNT=CURENT-CURENT+( V (I )+Y(2 )+EGK2/Y (1 )+2, *Y (2J-Y( 1 )*U( I J/U (4 ) 
1 _4,^Y( 21+U(3 )/U{ 4)-Y( 2)+EQK2y Y( 13+U( 2J /U{4) ) /2./7*7 
F (1 )= ( CURNT+V2-FLUXX *V4 )/0TA2 
F ( 2 )= { FLUXK + V 3-CURNT+ V 1 ) XOTA 2 
RETURN 
END 



SUBeOUTINE DlNCCRl ,CR2.,CR3,CR4»D.U) 

c this subroutine evaluated difflsi vities d and mobilities u for given 

C CON cent RAT IONS 

C CR*S ARE concentrations IN MCLE/LITER 

C THE NUMBERS ASSOCIATED WITH THE D»S AND U*S ARE CORRESPOND TO THAT WITH C»S 
DIMENSION D(5 ) ,U 14) .DCl (71 ,OC2(7 ) . DC3 <7 UCl ( 7 ) ,UCT1(7), 

+UC2(7).UC3(7) ,UCT3(7) .VISR(3) ,LCT2{7) 

COMMON /BLKl/FACMAT .TVISF,T£MPD 

DATA DCl/ 1 - 22223653 . 13108479 > 4 , 175 42244^-1 .8581188 4. 

1 .2323907, .02796468,-. 00566855/ 

DATA DC2/ 5. 13441 054 *5. 60764899,-6 .31470563,1 .45324064, 

1 1.281 £9542,-, 740008 12, . 10684605/ 

DATA EC3/ .90717802,-1 . 10332668,- . 29668964 ,4 . 8 0486975, 

1 -4.31019578,1 .25519037 10133989/ 

CATA UCl/44.35165773,-45. 91748832, 40. 40956334 ,-21 .182183598 « 

6.02671557,- .87455289, .05072426/ 

DATA UCT 1/2. 4 1549965, 1 .13364338,-2,13 873683 ,.71171359, 

06922713, -.0069561, .GO 12 12/ 

DATA UC2/198. 1 0978722 , -1 00 , 1 9543755 , 59 . I 0497144, -17 .69603863 , 

* -3.23141665, 2 . 46 1 34 1 03 . 30 9 1 23 1/ 

CATA UCT2/14 .6534768,-4 ,30 70 8816 , 4 . 26 289263 ,- 2. 78870985, 

^ .28474868,. 15513752, -.02804-015/ 

DATA UC3/34.3951439,-47 .22174887, 27.5656625,5.09830616, 

1 -1 4. 04491884, 5. 79626865, -.75395075/ 

DATA UCT3/2. 43698943, .38048952,-2. 07243995, . 1 9862643 , 

1 .89001 14, -.43544337, .05981754/ 

CR1=AM1N1 (CRl , 14. 4.) 

CR2=AMIM(CR2, 14 .4,) 

CR3=AMIM (CR3 ,14.4) 

CR4=AMIN1(CR4, 14.4) 

CR1=AMAX1 (CRl , .1 £-7 ) 

CR2=A MA X 1 { CR2 , . 1 E-7 ) 

CR3-= AM AX 1 ( CR 3 , . 1 E- 7 ) 

CR4= A MAX 1 ( CR4 , . 1 E-7 ) 

VISM=V1SC0S(CR4) 




DO se J=1.3 
D{ JJ=0» 

U<33=0. 

58 CONTINUE 

VISRtl )-VISCOS (CRl )/V ISM+TV ISF 
VISR( 2)=VlSC0S(CR2)yVI£N*TVISF 
V1SR(3)-V ISCOS( 2.^CR3)/VISM^TVISF 
C D(l) IS EVACUATeD AT T=77 F IN PURE SQL'N 

DO 61 1=1,7 

C(1 ) = D(1 )+CCl II)^CR1*4( 3-1 ) 

61 CONTINUE 

C D(2) IS EVALUATED AT T=77 F IN PURE SOL'N 

DO 62 1=1,7 

D ( 2 )=D { 2 ) +DC 2(1) *CR2*i^ ( I - 1 ) 

62 CONTINUE 

C 0(3) IS EVALUATED AT T=77 F IN PURE SQL'N 

DO 63 1=1,7 

C(3 J = D(3 )-^DC3 (I)^CR3 + =4'( I-l ) 

63 CONTINUE 

C D*S ARE EVALUATED AT THE CELL TEMPERATURE IN MIXTURE IN MATRIX 

DO 64 1=1 ,3 

0(1 )=0( I )*FACMAT*VI SR (11/. 258 

64 CONTINUE 

0 (4) =(CR1*1 .2/1 . 1Q4D (1 )+CR2*i .2/5 .25’S‘D ( 2 )+2 . *CR3*1 .2/ .923^ 

=«^C(3) )/CR4 

0(5 J = 1 ,a*Q m93TS/V ISM+TV ISF-tFACMAT/ .258 
C U(l) IS EVALUATED AT THE CELL TEMPERATURE IN PURE SOL'N 

DO 65 1=1,7 

U( 1)=U(1) + (UC1 (I )+UCTl ( I )^TEMPD)=«<CR1*^( I-l ) 

65 CONTINUE 

C U(1 ) IS EVALUATED IN MIXTURE IN MATRIX 

U( I )=U( 1 ) *VISR(1 ) /TVI SF4FACMAT 
C U(2) IS EVALUATED AT T=77 F IN PURE SOL'N 

DC 66 1=1 ,7 

U( 2)=U { 2,) + (UC2( I ) + UCT2( I )* TEMPO ) ♦CR2*^ ( I-l ) 



66 CONTINUE 

C U(2) IS EVALUATED AT THE CELL TEMPERATURE IN MIXTURE IN MATRIX 

U( 2) = UC 2)^VI SRI2) *FACMAT 

C U(3) IS EVALUATED AT THE CELL TEMPERATURE IN PURE SOL* N 

DO 67 1=1 *7 

UI 3)=U(3)+(UC3( I }+UCT3( I ITTEMPDl *CR3^+( I-l ) 

67 CONTINUE 

C U(3) IS EVALUATED IN MIXTURE IN MATRIX 

U{3)=UC3)+VISRI 3) /TVI SF^tFACMAT 

U (4)=(CR1^44 .9/44 .5^U ( 1 )+CR2*44 . 9/ 193 .0 W ( 2 )+2»^CR3*44 .9/ 
+34.7+UC 31 )/CR4 
RETURN 
END 
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SUBROUTINE DATAKTCEEL) 

THIS SUBPROGRAM INITIALIZED THE DIMENSIONLESS PHYSICAL AND CHEMICAL PARATERS 
C NEEDED FOR SIMULATIONS 

C FOR PROCEDURES OF NON-0 1 MENS IQ NL IZATI CN , SEE THE INTERNAL NCTE REFERRED 
C TCELL IS THE ELECTROLYTE TEMPERATURE 
C eOK2C= K2 IN the CATHODE ZONE 

C EQK2M= EQUILIBRIUM CONSTANT K2 IN THE BULK ELECTROLYTE ZONE 
C £QK2A= K2 IN THE ANODE 2CN£ 

C EQK1M= EQUILIBRIUM CQNSTAIsT Kl IN THE BULK ELECTROLYTE ZONE 
C EGK1C= Kl IN THE CATHOCE ZONE 

C EQK1A= Kl IN THE ANODE ZCNE 

C CKDS= SQUARE ROOT OF KCD5 IN THE CATHCOE ZONE AT A GIVEN C4 
C AKD= = SQUARE ROOT OF KAC5 IN TEE ANODE ZONE 
C HA- HYDROGEN ZONE MASS TRANSFER COEFFICIENT FCR CG2 
C HC= AIR ZONE MASS TRANSFER COEFFICIENT FOR CD2 
C VELA= HYDROGEN FLOW RATE 

C VELC= AIR FLOW RATE 

C HENRYC= HENRY CONSTANT FOR THE CATHODE ZONE 

C HENRYA= HENRY CONSTANT FOR THE ANODE ZONE 

C RAK= THE FACTOR B DEFINED IN THE NOTE 

C FKC= TEMPERATURE FACTOR USED FOR CALCULATING PH 

C TVISF= TEMPERATURE EFFECT ON THE VISCOSITY OF ELECTROLYTE 

C TEMPD= TEMPERATURE EFFECT ON SOME D IFFUS I VITIES 

C CKDS= TEMPERATURE EFFECT OF REACTION RATE CONSTANT 
C 

DIMENSION T ( 14 ) ,V ISCQW{ 143 

COMMON EaK2M«EQKlM'.EQK2C»EQK2A.£GKl A,EGK1 C, CK DS . AKD , HC » H A » V ELC * 
^VELA.HENRYC* HENRY A, RAK*PKC,CFM, 1C 
COMMON /BLKl/FACMAT ,TVISF. TEMPO 
COMMON /BLK3/VHV 

C THE DATA OF WATER VISCOSITY AT DIFFERENT TEMPERATURE 

DATA Ty40 *-.45. .50. .55. .60. .65..7 0* .75 . ,77., 60 ..85 . .90 . ,95* ,100 ./ 

DATA VISCOW/ 1 .54 6, 1.43S. 1 .307, 1 *21, 1.122,1, G44, . 9753 ,* 9135 , 

>s‘.8904,.8577, .8071 ,.7612 , .7194, .6812/ 


C TEMPERATURE IN CENTIGRADE 
TK= (TCELL-32 . )y 1 . 8 

HEN=(G. 98 1035+0 ,04210 2^TK+3 * 19 8 1 3SE-4+TK+*2) * 1 , E4 
C TEMPERATURE IN KELVIN 
TK=TK+273.ia 

C EVALUATION OF REACTION RATE CONSTANT 
REACK=10»**< 13. 635-t 2ES5./TK ) J 
TEMPD=(TCELL~75. )/10 . 

PKStf=4471 , 33/TK-6. 0346+ ,. 017053+TK 
PKl=3 404.7iyTK- 14.8435+0. C327£6*TK-FK!* 

FK2=2902.39/TK-6 .4980+0 .0 2379 ^i'TK-FKW 
IFITCELL .GT. T(D) GG TC 71 

TVISF=. 8904/1 {VISC0^I2>-VlSC05«ri II )/(T{ 2)-T( 1) )4={TCELL-Tt 1) 1 + 
*V1SC0W<1)) 

RETURN 

71 DC 72 J=2,14 

1F( TCELL ,GT. T(J)) GG TO 72 

J 1=J-1 

J2=J 

GO TO 73 

72 CONTINUE 
Jl=13 
J2=14 

73 TVISF= .3904/< ( VIS COW (J2I-V IS COW ( J 1 ) J/( T ( J2 )-T ( J 1) TCELL- T( J 1 > ) + 

*VISCOW( J1 ) ) 

AW=0.75 

FKC=PKW+ALOG10 (7.2/ AW 1 
EQKl=10..+4-PKl 
ECK2=10 .5}‘*PK2 
EQKlC=E<3Kl/23a, 

EQK1M=EQK 1/144. 

EGKl A=EQK1/51 . 

EQK2C=EGK2/58. 

EQK2M=EQK2/32. 

E0K2A~EQK2/*. 1 



HEIsiRYC=l ./HEN/1.686/7.2 
HENRYA=1 ./HEN/1. 190/7.2 
HC= { 2 . 39 5 2+ . 1 8 8 5 * CFM ) 4 h EN 
VELC=1. 3654c CFM4HEN 
FA=5.60*HEN 

VELA= .023664HEN 4VHV/0 .067 

CKDS=SQRT {REACK/{ 10. 4^4 ( 13. 635-1 2895. /298. ) ) )) 

RETURN 

END 



FUNCTION VISCOSTCS) 

THIS SUBPROGRAM EVAL.UATES THE ELECTROLYTE VISCOSITY FOR GIVEN CR4- 
VISCOS= .8937 -1 .0093£424:CS+£.37e63884^CS5i=^2-l. 13294407* 

*CS»*3+. 1762541 64=CS**4 
RETURN 
END 
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3U8ROUT INE ENERGY(NTSC , VOLT^CURENT ,71 N,TOEWI N ,CFiV, ,TELEC, CSCATH 5 
C THIS SUBPROGRAM SIMULATES THE HEAT AND WATER TRANSFER IN THE CELL 
C 

C WVP(T) IS THE WATER VAPOR PRESSURE OF PURE HATER AT TEMP. 1(1) 

C SVF(I,J) IS THE VAPOR PRESSURE OF TMAC SOLUTION WITH CR4=24CS(J) AT TEMP. 
C T{ I ) 

C RH07= DENSITY OF AIR, LE/FT*+3 
C CP= HEAT CAPACITY OF AIR, BTU/LB-F 
C AJ~ CONVERSION FACTOR, JOULE/CAL 

C HR= HEAT OF REACTION FOR H2(G) -f 1/202IG) = H20(G), CAL/G-MOLE 

C F= FARADAY'S CONSTANT, A NFS-SEC/G-ECU I V AL ENT 

DUCTS=NO. OF CHANNELS IN AIR COMPARTMENT 
C XAC= LENGTH OF THE AIR CHANNEL, FT 
C YAC= WIDTH OF THE CHANNEL, FT 
C ZAC- HEIGHT OF THE CHANNEL, FT 
C LIMNU= LIMIT NUSSELT NO. 

C AREA= ACTIVE ELECTRODE AREA, FT*’S=2 

C AREAMT= MASS TRANSFER AREA, FT%A2 

C AREAHT= HEAT TRANSFER AREA, FT*12 

C AREACS= AIR CHANNEL CROSS-SECTI QNAL AREA, FTY*2 
C W= AIR MASS FLOW RATE, LS/HR 

C 1= CURRENT, AMPS 

C N 6= MOLAR RATE OF WATER PRGDUCTICN, LE-MOLE/HK 
C CH= EQUIVALENT DIAMETER GF AIR CHANNEL, FT 
C QGEN= RATE OF HEAT GENERATED BY THE CELL, BTU/HR 

C PH20= WATER PARTIAL PRESSURE IN AIR FLCW, MMHG 

C T/»1R(J)= AI R TEMPERATURE IN T HS SUBCELL J 
C TELEC(J}= ELECTROLYTE TEMPERATURE IN THE SUBCELL J 
C MU7= AIR VISCOSITY, LB/HR-FT 
C 057- C02 DIFFUSIVITV N AIR, FT^=«'2/HR 
C D67= VMTER DIFFUSIVilr *N AIR, FT**2/HR 
C SC= SCHMIDT NO. 

C PR= PRANDTL NO. 

C HTC= HEAT TRANSFER COEFFICIENT, 8TU/HR~FT*^2-F 
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C MTC5= MASS TRAM SFER COEFFICieNT FOR CC2, FT/HR 

C MTC-6= MASS TRANSFER COEFFICIENT FOR WATER, FT/HR 

C PE= ELECTRCfLYTE WATER VAFCR PRESSURE, MMhG 

C CSCATH(J)= CATICN CONCENTRAT XiJh4 IN THE CATHOLYTE -GF THE SUBCELL J 
C 

REAE LIMNU,MU7*K7, I ,N6 

OIMENS ICN T (13 J.WVPI 13) ,T ELEC ( 40 ) ,CSC ATH( 40 ) , SVP ( 13.11) ,CS( 11 ). 

1 JT( 2) ,PET( 2) ,CST(2> , X (1 1 ) . TA I R ( 4 1 ) , PH20(41 ) ,RH( 1 1 ) 

DAT A T./'90r.S£»,S0..7S»»7G..E£. »GC.'*S5*»5C. »45.»40'.»35.» 

1 30./ 

DATA WVP/3 6, C7 . 30 . 7 5 . 26 , 27 ,22 . 24 . 1 8 .76 , 15 .77 . 1 3. 29 * 1 1 .09 

1,9.22,7.7,6.30,5.17,4.23/ 

DATA RH07/D.0aC8/. CP/0.24/, AJ/4,186/, HV/106O./, HR/ -57797 .9/ . 

1 R/5S5./, F/96500./ 

DATA RH/1 . , .970,. 949, . 913,. 850, .763, .658, .551 , .424, .324 , . 30/ 

DATA CS/0., . 49, . 74 , . 9 S , 1 . 25 , 1 . 52 , 1 . 79 ,2 , 06 ,2. 33 ,2.61 , 2 .75/ 

DATA SVP/36 .07,30 .75 ,26 .27 , 22 . 24 , 1 8 . 7 6 , 1 5. 77 , 13. 29 , 1 1 . 09 , 

1 54 0 c I.- 3 5. 02, 29. 87,25.55,21 .61,18.23.15.32,12.91, 10.77,540., 

2 34.27,29.20,24.98, 21.13,17. £3. 14.98, 12. 62 , 10 . 54 ,54 0. , 

3 32 .97,28.10, 24.02, 20 .33, 17. 15, 1 4 . 4 1 , 12 .14, 1 0 . 14 , 540 . , 

4 30 . ?0, 26. 18 ,22. 37. 1 £.93 , 15. 96, 13.42.1 1.31 .9.44.540. . 

5 27.66. 23.48, 20.07, 16 .99, 14 .33 , 12.04, 10.14,6.47,540. , 

6 23c8l ,20.28,17.36.14.68 ,12.38, 10 .41 , 8.76,7.32, 540. , 

7 19.67, 16 .98. 14.52, 1 2 . 27 * 1 0. 3 6 , 6. 7 1 , 7 . 34, 6. 1 1 ,540. , 

8 15.33,13 .07,11 .17,9. 45, 7. 98, 6. 70, 5. 65,4. 71, 540., 

9 il ,72,9. 99, 8. 55, 7. 23 ,6. 09 ,5. 13 ,4.32,3 .60 ,540 . , 

1 9 .6,7 .8,6.3, 5. 05, 4. .25, 3 .5. 2. 9, 2. 3, 540. / 

DUCTS=96. 

XAC=6./12 . 

YAC=0.2/12. 

ZAC=G, 19/12. 

AREA =1. 

AREAMT=YAC4XAC40UCTS 

LIMNUi52.89 

AR6AHT-2.04(Y AC+2AC) 4XAC4DLCTS 


^REACS=ZAC=i'YAC 
DH=2.^AREACS/( YAC+ZAC) 

W=CFM*60 .*RH07 

I-CURE^T*AREA 

N6=I/2./Fy453.*36G0. 

W6=N6*18 .02 

QGEN=-I* ( HR/2 ./F+VOLT/ A J ) 4:3600 ./252,. 

TAIR( 1 )=TIN 

CALL LAGRAN(TDEWIN. T .IsVP . 13 ,Ph20IN , 13, 13) 

PH20( 1 )=PH20TN 
TSCN=FLOAT(NTSC) 

DC 821 J=1,NTSC 

TAIRC J+1 )=TATR( J)4-QGEN/CP/Vi/TSCN 
321 CONTINUE 
835 DO 94 JC=1 , NTSC 

TAVE=(TAIR( JC)+TAIR{ JC+1) )/2. 

MU7=1 (TAVE+460. )/492 . >441 .54697 ./< TA VE+665 . ) 4 0. 0 170942. 4 2 
K7=D .01 40 4(TA VE-32 . ) /I 80 .40 .0043 

D57= . l64/( 2.54442)/ 14A. 43600. 4 ( ( TA VE+460. ) /537. ) 441 . 832 
D67= .256/C2 .54 442) / 144 . 43600 . 4 ( {TAVE+460. )/537. ) 442 .334 
R eO=DH4CF M/DUCT S/ARE ACS4RH07/MU7460 . 

FR=CF4MU7/K7 

SC5=MU7/RHQ7/DS7 

SC6=MU7/RH07/D67 

CCNST=(0, 003+0 .0194ZAC/YAC )4RED4PR4DH/XAC4TSCN/FL0AT( JC) 
CaNS5=(0. 003+0. 0394Z AC/ YAC)4RED4SC54CH/XAC+TSCN/FL0AT(JC) 
CONS6-(0 , 003+0. 0394ZAC/YAC)4RED4SCe4DH/XAC4TSCN/FLOAT( JO 
HTC=LIMNU4K7/DH4(1 .+CCNST ) 

MTC5=LIMNU4D574( 1.+C0NS5) /DH 
KTC6=LIMNU4D674 (1 ,+CaNS6 )/DK 
TEL£C{ JC)=TAV£+QGEN/HTC/AREAHT/TSCN 
TELECC JC)=AMAX1{TAIR( JC+1) ,TELEC( JO) 

PH2Q( JC+l ) = (PH20 ( JO/ (TAIR( JO+460 . ) +N6/TSCN4R/CFM/60 . ) 4 
1 {TAIR{ JC+l)+460. ) 

PH2QAV=(PK20 ( JO+PK20 (JC+1 ) )/2. 




O to- 

I® 


PE=FH2CAV+N6*R* CTAV E+46€ . )/AREAMT/MTC6 
PE=AMAXl< PH20AV.PE) 

C USE LAGRANGIAN INTERPOLATION TO DETERMINE WVP AT TELECIJC) 

CALL LAGRAN(TELEC( JC) *T,WVP» 13,WVPTE, 13»13) 

C CALCULATE WATER ACTIVITY, AW, OF THE CATHGLYTE IN THE SUBCELL 
AW=PE/WVPTE 

I FI TELECIJC) ,LT. 90.)GG TO 111 
WVPTE=3E,ll’ScEXP< • 03 1 ( TELEC ( JC)-90,) ) 

AW=PE/WVPTE 

CALL LAGRANIAW,RH,CS,11,CSE,11 ,11 ) 

GC TC 93 
111 CONTINUE 

C determine a set of TIJTIIJ) AND T(JT{2)J SUCH THAT JT ( 1 J = JT ( 2 ) -1 AND THAT 
C T(JTU)), TELECIJC), T(JT(2)) ARE IN DESCENDING ORDER 
DO 85 J=l,12 

IFIT(J) .LT, TELECIJC)) GO TO £6 

85 CONTINUE 
J=J + 1 

86 JT(1)=J-1 
JTI2)-J 

IFIJTIl) .GT. 0) GO TO SS 
PETIl )=WVP(1 ,)*AW 
DO 87 L=l.ll 
X I L)=SVP( 1,L ) 

87 CONTINUE 

CALL LAGRANIPETI 1) ,X,CS,11 ,CSE,11 ,11 ) 

GG TG 93 

88 IFIJTI2) -LE. 12) GO TC 90 
PETI2)=WVP( 12)#AW 

DO 89 L=l,ll 
X( L )=SVPI 12, L) 

89 CONTINUE 

CALL LAGRANI PET 12) ,X,CS,1 I ,CSE, 11,11) 

GO TO 93 

C DETERMINE TFE CATION CONCENTRATIONS OF A CARBONATE SOLUTION WITH WATER 



C ACTIVITY AW AT THE TEMPERATURES T(JT(D) AND T(JT(2)) 

90 DO 92 J=1 ,2 
JTJ=JT( J J 

PET ( J) =WVP f JT J ) «!^AW 
00 91 L=l,ll 
X{L)=^SVP(JTJ»L) 

91 CONTINUE 

CALL LAGRAN(PET( J) ,X,CS,11,CST( J) ,11,111 

92 CONTINUE 
J1=UT( 1) 

J2=JT{2) 

C INTERPOLATE FOR CATION CCNCENT RAT I CN OF THE CATHOLYTE AT TELEC(UC) 
CSE=CST( 1 )+(CST(2)-CST( 1 1 ) /{ T( J2 ) -T ( J1 ) ) * ( TELEC 1 JC) -T ( J1 ) ) 

93 CSCATFC JC)-CS£*2, 

94 CONTINUE 
PH200U=PH20CNTSC+ 1) 

CALL LAGRAN{ PH200U.WVP,T, lo.TDEWOU, 13, 131 
WRITE (6, 94001 TDEWOU 

9400 FORMAT lOX, 'OUTLET AIR D ,P , = • ,3X ,F6 . 2 ,3X , • DEG. F » 1 
WRITE! 6,9401 1 

9401 FORMAT('0*,IOX, • NO. OF SUBCELL TAIRIN TAIRCUT TELEC P 
1H201N PH20OUT CSCATH*} 

00 95 JC=1,NTSC 

WR ITE (6 , 9 40 2 1 JC, TA IR { JC 1 . TA IR { JC4 11, TELEC { JC ) ,PH20 ( JCl , PH2C ( JCl-1 1 
1, CSCATH (JC) 

9402 FORMAT! 16X, 12, 1 1 X, 5( F 6 . 2 ,4> 1 .F6-2) 

95 CONTINUE 
RETURN 
END 



SUBROUTINE L AGR AN ( X A * X , Y , N, ANS* NX ,NY) 

C THIS IS A L.AGRANGIAN INTERPOLATION ROUTINE 
DIMENSION X(NX},Y(NY> 

SUM=0.. 

DO 3 1=1 *N 
PROD=YI U 
DO 2 J=1,N 
A=X ( I 1-X( J ) 

IF(A) 1.2,1 

1 E=(XA-X(J ) ,)/A 
PRCO=PRDDY£ 

2 CONTINUE 

3 SUM=SUM+PROD 
ANS=SUM 
RETURN 

END 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE CRANE CRANOOlO 

IMPLICIT REALMS (A-H.O-ZJ 

HELP IS LOGICAL VARIABLE SET TRUE IK DERI V WHEN 

A DISCONTINUITY OCCURS AND IT IS NECESSARY TO DECREASE THE TIME 
INCREMENT AND RESTART THE INTEGRATION BEFORE CONTINUING. 

THE Y AND F ARRAYS AND X ARE RESET TO THEIR LAST GOOD VALUES 
AND I NIT SET EQUAL TO ZERO BEFORE PROCEDING WITH NEW H. 

ARGUMENTS Or THIS SUBROUTINE ARE ALL CCKTAINED IN THE FOLLOW- CRANG020 
COMMON STATEMENT* LABELED CRANEC. DRIVER PROGRAM MUST CCNFGRMCRANG030 


THE 
IKG 

IN ORDER AND DIMENSIONS WITHIN COMMCN/'CFANEC/.THOUGH NOT NECES- 
SARILY IN THE NAMES OF THE VARIABLES. 

LOGICAL lONT.INDONX.IKPCsiELP 

COMMON /CRANEC/X. Y( 40 ) *F( 40) ,PREC . XHAT ,H . N *NB ( 40) ,MI N «I NIT .IGNT * 
:^INDONX,IKPC 
I EL P= .FALSE. 


CRAN0040 

CRAN0050 

CRAN0060 


X 

Y 

F 

H 

N 

NB 


CRAN0080 

CRAN0090 

CRANOlOO 


MEANING OF THE ARGUMENTS ARE AS FOLLOWS 
THE INDEPENDENT VARIABLE. 

ARRAY OF DEPENDENT VARIABLES. 40 MAXIMUM. 

ARRAY OF DERIVATIVES, FUNCTIONS OF X AND THE Y VECTOR. CRAN012G 

STEP SIZE IN X. MUST BE NON-ZERO. USUALLY MAY BE ALTERED CRAN0130 

1^TERNALLY BY CRANE* AFTER FIRST BEING SET + OR - EXTERNALLY .CRANO 140 
THE NO. OF Y S AND F S IN THE SYSTEM, USLALLY THE NO. ALSO CRAN0150 
OF COUPLED ORDINARY DIFF. EQUNS. CRANO 160 

IF +, HOLD ABSOLUTE TRUNCATION ERROR EST.CRAN0170 
BELOW A SPECIFIED LIMIT =1 0 . ( -FREC) CRAN0180 
IF 0, DITTO FOR RELATIVE TRUNCATION ERH0RCRAN0190 


AN ARRAY OF FLAGS. 


MIN 


INTT 


THE 

THE 

SET 


ESTIMATE. 

IF DISREGARD THE TRUNC. ERROR. EST. 

FOR THE V $ WHOSE FLAGS ARE NEGATIVE. 
MINIMUM NO. OF STEPS BETWEEN OCCASIONS FOR DOUBLING H. 
VALUE OF MIN MAY NOT EXCEED 10, NOR EE LESS THAN 3. 

TO ZERO BEFORE FIRST CALL, CAUSING THE SR TO INITIALIZE 


ITSELF. thereafter* INIT IS BUcMPED INTERNALLY. EACH STEP. 


CRAN0200 
CRAN0210 
CRANO 220 
CRAN0230 
CRANO 240 
CRAN0250 
CRAN0260 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DQNT A FLAG WHICH IF TRUE FORBIDS ANY CHANGE IN H TO BE MADE 
INTERNALLY BY THE SR, (UNLESS £^CC^X=T) 

FRBC A POSITIVE CONST* APPROXIMATELY EQUAL TO THE NO. OF SIGMFI- CRAN0290 
CANT DECIMAL DIGITS OF LOCAL PRECISION* (BASED ON TRUNCAT lONCRANO 300 
ERROR ESTIMATES MADE INTERNALLY. EFFECTIVE ONLY IF D0NT=T.) 

ENDONX A FLAG WHICH IF SET TRUE IN CALLING ROUTINE AND ACCOMPANIED 
BY A VALUE OF XHAT .GT. X AT 1MT=0 FOR h . GT . 0. (OR XHAT 
,LT. X IF H .LT. 0.) WILL ENABLE SR TO RETURN TO CALLING 
PROGRAM WITH X=XHAT AND ENDONX SET FALSE AS A FLAG TO THIS 
EFFECT* DONT WILL BE IGNORED IN THIS INSTANCE. 

XHAT VALUE OF X. IDEPENDENT SEGMENT, AT WHICH A DISCQNTINUI TY 

OCCURS. USED IN CCNJUNCTICN WITH ENOCNX TO HAVE X=XHAT AT END 
OF A STEP. 

OKPC A FLAG TO INDICATE THAT AT LEAST ONE PREDICTOR-CORRECTOR STEP 
WHICH SATISFIES THE ERROR CRITERIA HAS BEEN TAKEN. 

CRAN0320 

THE PARTICULAR PREOICTOR-CCRRECTOR ALGORITHM IS THAT PUBLISHED CRAN0330 
BY CRANE AND KLOPFENSTEIN IN J .A .C.M* VOL 12, PAGES 227-241, APRi LCRAN0340 
1965. IT IS OF FOURTH ORDER, I.E., TRUNCATION ERRORS ARE OF 0RDERCRAN0350 
H^*5. THE ALGORITHM MAS DEVELOPED TO MAXIMIZE THE RANGE OF STEF CRAN0360 
SIZE CONSISTENT WITH ABSOLUTE STABILITY, AND IN A LOOSER SENSE, TOCRAN0370 
HAVE A Gear RANGE OF RELATIVE STABILITY. THESE RANGES, EXPRESSED CRAN0380 
AS H, NORMALIZED BY MULTIPLYING BY PARTIAL DERI V OF Y(=G(X,Y)) WITHCRAN0390 
RESPECT to Y, ARE CRAN0400 

FOR ABSOLUTE, C. GE . HBAR . GE . ( -2 .4809 ) CRAN0410 

FOR RELATIVE, INFIN ITY . GE . HB AR .GE . . 446J CRAN0420 

THE STARTING PROCEDURE IS THE RUNGE-KUTTA VARIANT PUBLISHED BY CRAN0430 
S* GILL, CAMS. PHIL. SOC, PROC. , 47, PS6,<1S51) CRAN0440 

CRAN0450 

TO USE THIS ROUTINE, A SUBROUTINE CALLED DERIV MUST BE FFCV I -CRAN0460 
DEO WHICH USES LABELED COMMGN/CRANECZAND CALCULATES THE F VECTOR. CRAN0470 
THE FIRST CALL OF CRANE WILL CNLY OBTAIN F S.AND SET INIT=1. EACH CRAN0480 
SUBSEQUENT CALL WILL ADVANCE ONE STEP AND UPDATE X,Y AND F, EXCEPT CRANO 490 
THAT AFTER THREE STEPS, CRANE MAY HALVE H AND RETURN TO INITIAL CRANOSOO 
CONDITIONS WITH lNlt=4. (THIS BEHAVIOR WILL RECUR UNTIL PRECIS IQNCRAN0510 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


IS SATISFIED ACCORDING TO AN ERROR ESTIMATE AFTER THREE STEFS.) 


DO NOT 
EXTERNALLY. 


CRAN0S20 
CRAN0530 
CRAN0540 
CR A NO 550 
CRAN0560 
CRAN0570 
CRAN0530 
ENOTNCRAN0590 
CRAN0600 

CHANGE X.Y.F.N OR H EXTERNALLY UNLESS YOU SET INIT=0. CRAN0610 
HOWEVER. N3 AND FREC MAY EE CHANGED AT ANY TIME. SET-CRAN0620 


PRIOR TO FIRST CALL, SET INIT X AND Y AND DEFINE THE VALUES 
OF THE NB VECTOR, ALSO DEFINE VALUES OF N , MI N ,I NI T( -=0 ) , lODKT , AND 
PREC. PREC SHOULD ADVISEDLY SE FROM 2.5 TO 7. (NOT NECESSARILY 
AN INTEGER), THOUGH THERE WILL BE NO EFFECT IF IDCNT.NE.O 


printout and TERVjINATIGN are NOT ACCCMFLIS&EC BO CRANEC 


H FCR hlh 
POSITIVE ON 


C 
C 

c 
c 
c 
c 
c 


TING INIT.LT.O will RESULT IN INHIBITING INCREASE IN 
STEPS, BUT WILL PERMIT DECREASES IN H. INIT WILL BE 
EX ITT IN G CRANE AGAIN. 

DIMENSION Q (40) ,Y1 ( 40 ) , YO 140 ) , FO ( 40 ) , F 1 { 40 ) ,F2t 40 ) ,F3 1 40) ,F4( 40 ) , 
1F5{ 40) ,Y2(40) ,E( 10) ,Y3(40) , Y4 ( 40 ) , Y5 ( 40) 

EGUIVALENCEt Y5 I 1 ) , Q { 1 ) ) , { R , T ) 

Y5 AND Q ARE NOT BOTH NEEDED AT THE SAME INSTANT. 

DATA VNU /I .S2585S6E- Csy 

VNU MAXIMUM ROUNCCFF ERROR ON UNIVAC 11G8 CALCULATED BY 

V=ALPHA^BASE*^{ 1-IT) WHERE ALPHA =1. FOR ROUNDING MACHINE 
AND .5 FOR TRUNCATING MACHINE. BASE IS THE NUMBER BASE FCR 
MACHINE (2 FOR 1108). IT IS NUMBER OF DIGITS IN MANTISSA OF 
the WORD IN MEMORY ( 27 FOR 11 C6 ). 

FIRST executable NEXT 


CRAN0630 

CRAN0640 

CRAN0650 

CRANG660 

CRAN0670 

CRAN0680 

CRAN0690 


:RAN071 G 




IF( INIT)1 ,2*3 
1 INIT=0 

E( MIN)=10 . 

2 IF ( N .GT. 0 .AND. N .LE. 40 ) GC TO 5 
WRITE( 6,200) 

200 FORMATI54HON WAS NOT A POSITIVE I NT EGER .LT . 40 , ON CALLING CRANE.) 
CALL EXIT 

IF ( MIN .LT. 3 ) MIN=3 


CRAN0720 

CRAN0730 

CRAN0740 


CRAN0790 

CRANOaOO 


5 



:-r.' rrnnv«r.v — '^. 


IF ( MIN :.GT. 10 ) MIN=I0 

IF ( INOONX .AND- {XH/iT-XiyH .LT. 0.) GO TO 2040 
beta =-DLOG10( VNU) S 
70 IF ( PFEC-EETA ) 11.13.12 

12 PFEC = BETA 

11 IF t PREC .LT. 2. ) PR£C=4. 

C STATEMENT 11 ASSUMES THAT THE OVERSIGHT OF SPECIFYING NO PRECI- 

C SIGN. OR PRECISION BELO^i 2«0 DECIMALS. SHOULD CALL FOR 4 DECIMALS 

C VIA INTERNAL CORRECTION, ALSO. PREC ABOVE MACHINE ACCURACY IS 

C REDUCED TO MACHINE ACCLRACY.. 

13 PRECN=PREC 

T=l 0. *’!'<-PREC) 

B0UND8=T^16.21S66 
BCUNOA=T* .05 
8GUND=BOUNDA4FLOAT( MI N) 

I FI IN IT U 155, 1155, 16 
1155 CALL DERIV 

IF (lELP) GO TO 2000 
ASSIGN 13 TO KT 
C:**^%*RUNGE-KUTTA SETUP 
17 K=0 

IF ( INCONX .AND. (X HAT-X )/ ( 4 .*h ) .LE . 1. ) H=( XHAT-Xl/4. 

ASSIGN 55 TO IS 

H1-.5^H 

KBC-3 

IKPC = .FALSE. 

GO TO KT, 118,20) 

C IF PROGRAM HAS BACKED UP TO START OF A R-K GO TO 18. 

C IF PROGRAM HAS BACKED UP TO LAST PC THEN GO TO 20. 

C THIS ’a'ILL permit PRINTOUT TO 'START OVER*. 

C ALSO FIRST CALL WITH 1NIT=0 WILL ALSO PASS THRU 18. 

C4**4*EKD R-K SETUP 

C THIS IS THE NORMAL CALL CF DERIV WITH A NEW PC Y VECTOR. HENCE 

C A USEABLE VALUE OF F IF HELP=FALSE. 

15 CALL DERIV 


CRAN0900 

CRAN0910 


CRAN0930 

CRAN0940 

CRAN0950 

CRAN0960 

CRAN0970 

CRAN0980 

CRAN0990 


CRANIOOO 

CRANlOlO 
CRAN1020 
CR AN I 030 


CRANl 050 



ft ft 


CRAN1C60 


IF <IE1_P) GO TO 52 
18 1NIT~1+IN1T 

IF (INOGNX •AN0*DA8S( (XFAT-X)/XHAT ) *LT. VNU ^10.) 

INDGNX=.FALSE* 

19 IF ( IKPC .OR. INIT . £Q . I ) RETURN 

THE ABOVE IS THE ONLY RETURN FROM INTODE. NEXT IS REACHED WHEN CRAN1Q80 
INIT exceeds O. CRAN1090 


3 

IF { PREC - PRECN ) 70.16.70 


16 

IF { fONT .AND. IKPC 1 GO TO 2 13 
IF ( IKPC ) GO TO 212 


c 

20 

K=4 

CRAM 120 

213 

OO 22 I=i.N 

CRANl 130 


F3(I)=F2{I) 

CRANl 140 


F2( I) -FI <1) 

CRANl 150 


Fl{ I )-F0( 1 ) 

CRANl 160 


F0(1)-F(1) 

CRANl 170 


Y3t I)=Y2C 1) 

CRANl 180 


Y2(I1=Y1(1) 

CRANl 190 


Y1(I)=Y0(I) 

CRANl 200 

22 

Y0{I)=Y (I) 

IF ( KBC .EQ. 0 .OR. lONT } GO TO 23 

CRANl 210 

24 

KBC=KBC-1 

CRAN1230 


XN=X+H 

CRAN1240 


X=X+H1 

CRAN1250 


RX=. 2 928 9321 9 

CRAN1260 


ASSIGN 33 TO KS 

CRAM 270 

START OF ONE RK STEP 


37 

DO 29 1-1, N 

CRANl 280 


AK=Hsi«F(I} 

CRAM 290 


GO TO KS. {33.34.301 

CRAN1300 

33 

Y( I )=Y( I )+0.5*AK 

CRAN1310 


Q{I)=:AK 

CRANl 320 



30 

34 

29 


219 

251 

224 

1551 


36 


■C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GC TC 29 

Y( I )=(-Q( I )-Q (I J+AK ) *• 16 6666667 +Y {I ) 
GO TO 29 
R=RX*{ AK-G{ I ) ) 

QCI )={ ( RHhR+R3-RX4AK)+Q( I) 

Y 1 1 )=Y { I ) +R 
CONTINUE 

GO TO (1551, 224, 219, 25l),K 

RX=1 .7071 0678 

GO TO 1551 

ASSIGN 34 TO KS 

GO TO 1551 

ASSIGN 30 TO KS 

X=XN 

CALL OERIV 

IF ( lELP ) GO TO 2020 

THIS CALL OF DERlV IS CURING R-K CNLY 
K=K- 1 


IF(K) 18, 18, 37 
END OF CHE RK STEP 

DELTAX IS DONE. THE COOING IS A BIT 
THE LOOP BEGINNING AT 37 IS EXCTD 
MULTIPLICATIONS AND INDEXING. IT I S 


AS FOLLOWS 


A CC I 
VALUES 


TO S 
OF A, 


WHEREIN 


29 


36 


DO 36 J=1.A 
DO 29 I=1,N+1 
AK( I , J)=DELTAX»F( I ) 

R <1, J )=A( J )#AK( 

Y (I) =Y(I)+R( I,J) 

Q (I) =0(1 )+3.^ R( I , J)-C( J)*AK( I , J) 

X= Y ( N+ 1 J 

CALL DERIV (X,Y,F) 

CONTINUE 

RETURN 

. GILL, C A MB. PHIL. SCC^ 

B , AND C ARE BEST...* 


COMPLEX IN ORDER TO MINIMIZE 
4 TIMES BEFORE ONE STEP IN 
EGUIVALENT TO SLOWER CODE, 

J = 5-K 


PROC. , 47,P96 (1951), THE 


CRAN1330 
CRAN1340 
CRAN1350 
CRAN1360 
CRAN1370 
CRAN1380 
CRAN1390 
CRAN1400 
CRANIA 10 
CRANl 420 
CRAM430 
CRAN1440 
CRANl 450 
CRANl 460 
CRANl 470 


CRANl 480 
CRANl 490 

CRANl 510 
CRANl 500 
CRAN1520 
CRAM530 
CRANl 540 
CRANl 550 
CRAN1560 
CRANl 570 
CRANl 580 
CRANl 590 
CRANl 600 
CRANl 61 0 
CRANl 620 
CRAN1630 
CRANl 640 



c 

A( J>=0.5>1 .-SQRTFI .5 ),1.+SQRTF{.31,1 ./6 . J=l, 4 

CRAN1650 

c 

B{JI=1. ^ 1 .-SQRTF( . 5) , 1 .+ SQRTF( .5) ,1. /3. J=1 ,4 

CRAN1660 

G 

C( J)=0.5,l .-SQRTF( .5) ,1 ,+SQRTF( .5), 0,5 J=l,4 

CRAN1670 

c 


CRAN1680 

C 

c 

FOLLOWING IS REACHED FROM STMT 16 WHEN K0C=0. AFTER PRED-CORR 

CRAN1690 

c 

INITIALIZATION IS COMPLETE. 

CRAN1700 

212 

IF ( INDQNX .AND. ( XHAT-X ) Y ( 4 ) .LT. 1* ) GO 70 IT 

IF ( E(MIN) .GT. BOUNCA ) GO TO 28 


100 

DO 4 1=2, MIN 

CRAM 750 

4 

E(1 )=E(1 )4E< I ) 

CRAN1760 


IF(E( 1 I-B0UND)38.28*28 

CRAN1770 

c 

28 IS USUAL. 38 IS PREPARE TO DOUBLE DELTAX. 

CRAM 780 

28 

DO 39 1=1, N 

CRAN1790 


F5(I)=F4(I) 

CRANiaOO 


F4( I)=F3 ( I) 

CRANI810 


F3{ I )=F2( I ) 

CRAN1820 


F2( X)=F1( 1} 

CRAM 830 


Fill) =F0 < I ) 

CRAN1840 


F0(II=F <I) 

CRAN1850 


Y5( I )=Y4( I) 

CRANl 860 


Y4( I)=Y3C 1) 

CRAN1S70 


Y3I I)=Y2( I ) 

CRAM 880 


Y2( 1 )=Y1 { 1 > 

CRANl 890 


Y1(I)=Y0(I) 

CRAN1900 


Y0(I)=Y Ll'i 

CRANl 910 

39 

Y(I)=1 .547652*Y(I )+2.Q17204^Y2( I )-l .867503=^^71 (I )- .697353*Y3 ( 1 

CRAN1920 


1 ( .gaS50 £124^F( I ) + .89512 1303=!<F2( I 1- . 351 589 £ I ) -Fl £ 1) ) 

CRAN1930 

40 

X=X + H 

CRANl 940 


CALL DERIV 

CRAN1950 


IF ( IEl_P ) GO TO 20 30 

THIS CALL OF DERIV IS FIRST WITH PREDICTED Y VECTOR AND MAY BE 
REACHED WITH LAST STEP OF AN R-K STEP HENCE UNCHECKED OR WITH 
LAST STEP A GOOD PC, 

42 DO 10 1=2, MIN 


C 

c 

c 


CRAN1960 



c 


c 

c 


c 

c 

c 

c 


10 


72 

44 

46 

47 

48 
43 


50 


E( 1-1 )=E< I ) 

E( MIN 1=0 . 

DO 43 1=1 ,N 

T=YO( I )+AK+{ g,*F{ I) + 19.*F0{1 )+F2{ I )-5**FlU ) ) 
FMC=Y(I)-T 

AT THIS POINT, PMC IS PREO I CT OP-CCRRECT OP 
IF ( lONT 1 GO TO 43 
IF< NB(1 ))43 *44,47 


CRAN1970 
CRAM 980 
CRAN1990 
CRAN2000 
CRAN2010 
CRAN2020 

CRAN2040 


IF {Y(I )) 46,43,46 
FMC=PMC/Y(IJ 

IFICABSl PMC)-E(MIN) 143,43,48 
EIMlN )=OABS(PMC) 

Y(1)=T 

HERE, EIMIN) IS THE MAX. T.RUNC. ERROR EST 
WHILE THIS REST OF £ ARE ITS PREDECESSORS . 
IF<E(MIN)-BOUNDB)50 ,5C ,51 
GO TO IS, ( 15., 55) 

15 IS normal and calls ceriv, jumping to 1 

55 IS REACHED FIRST STEF AFTER A SEQUENCE 
WHEN largest ERROR EXCEEDS BOUND, STMT 51 
BY WHETHER THIS IS ON FIRST P-C STEP OR NOT 


CRAN2050 

CRAN2060 

CRAN2070 

CRAN2Oa0 

CRAN2090 

ABS. OR REL., OVER ICRAN2I00 

CRAN21 1 0 
CRAN2120 
CRAN2130 

8 TO STEP INIT+1 ,RETURNCRAN2140 
OF RKG STEPS. ICRC 0510 ) CRAN2150 
ISA BRANCH GOVERNED CRAN2160 

. 52 IS NORMAL. CRAN2170 


51 IF { .NOT. IKPC ) GO TO 59 
C 52 IS REACHED WFEN LARGEST ERROR EXCEEDS BOUND, 

52 X=X-H 

DO 54 1=1, N 
Y(I)=Y0(1) 

54 F(I)=FO(I) 

ASSIGN 20 TO KT 
53 IF { lONT ) GO TO 1000 
H= .625+H 

lELP = .FALSE. 

C NOTE - ALL STEP SIZE REDUCTIONS GC THRU HERE 

C EXCEPT CN ENDONX INITIATED REDUCTION. 

GO TO 17 

RETURN TO 17 SETS UP RKG INI TI ALI ZATI CN 


FNCTN OF PRECI SN.CRAN2190 

CRAN220D 

CRAN2210 

CRAN2220 

CRAN2230 


CRAN2240 


C 


CRAN2250 

CRAN2260 



u y 


55 IKPC = .TRUE, 

ASSIGN 15 TO IS 

GO TO 5 

C 

C ABOVE SEQUENCE COMPLETES I M T I AL IZAT I GN OF PRED-CGRR ROUTINE. 

C NEXT, FAILURE OF FIRST TEST AFTER RKG SEQUENCE CAUSES BACKUP 3 

C STEPS . 

59 X-=X-A.^H 

DC 60 1=1, N 

Yd )=Y3i I ) 

60 FCI3=F3{I) 

61 ASSIGN 16 TO KT 

GO TO 53 

THE FOLLOWING SETS UP DOUBLING OELTX 
REFER TO CARD 1130 FOR RETURN TO RKG INITIALIZATION, 

38 H=2.^H 

DC 64 1=1 ,N 
F2{I)=F3(I) 

F3I I )=F5( I) 

F0( I3=F (I) 

Y2( I )=Y3I 1 ) 

Y3{I)=Y5(I) 

64 Y0(.I) = Yd) 

C 23 IS REACHED AFTER 3 RKG STEPS, I E WHEN KBC=0,(REF, STTMT 223 

23 AK=, 041 6666667 Jf'H 
E(MIN3=10. 

S=H’i<2.03169 
DO 62 1=1 ,N 

C SAME AS STATEMENT 39 

62 YCI 3=1 ,547652’^Y { I 3 + 2 .0 1720 4*Y2( I 3-1 ,667503*7 1 ( I 3- ,697353*Y3( I 34-S* 
1( .985506 1 24^^tF { I 3 + . 8951 21303*F2( I 3 -.351 5 89G7 1 4F3 ( I 3 -F 1 (131 
GO TO 40 

C ABOVE CALLS DERIVATIVE WITH FIRST PREDICTED Y VECTOR. 

C 

2020 K = 3-KBC 


CRAN2280 

CRAN2300 

CRAN23 10 
CRAN2320 
CRAN2330 
CRAN2340 
CRAN2350 
CRAN2360 
CRAN2370 

CRAN2380 

CRAN2400 

CRAN2390 

CRAN2410 

CRAN2420 

CRAN2430 

CRAN2440 

CRAN2450 

CRAN2460 

CRAN2470 

CRAN2480 

CRAN2490 

CRAN2500 

CRAN2510 

CRAN2520 

CRAN2530 

CRAN2S40 

CRAN2550 

CRAN2560 

CRAN2570 



GRIGINAi; PAGE IS 
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X =X N- K 

GO TO (2021 ,2022, 2023)^K 

2021 DO 2024 1=1, N 

F ( I ) = FO ( I ) 

2024 Y(I) = YOCIl 
GO TO 61 

2022 DC 2025 .1=1, N 

F ( 1 ) = F 1 ( I ) 

20 25 Y(I) = YKI) 

GD TO 61 

2023 DO 2026 1=1, N 

F ( I ) = F2 < I ) 

2026 Y(I) = Y2(I) 

GO TO 61 
C 

2030 IF ( IKPC 1 GO TO 52 
GO TO 59 
C 

1000 WRITE (Or 1001) 

1001 FORMAT COTHE USE OF A FIXED STEP SIZE (DONT=TRUE) IS INCOMPATIBLE 
, WITH USE OF HELP IN DERIV AND IS GENERALLY USELESS.*/* RUN ABORTS 
• D. • ) 

CALL EXIT 
C 

2000 WRITE! 6, 200 1 ) 

2001 FORMAT! *0 CRANE WAS CALLED WITH INITIAL VALUES UNACCEPTABLE TO OER 
.IV. CASE ABORTED.*') 

CALL EXIT 
C 

2040 INDONX = .FALSE. 

WRITE(6,2041 ) 

2041 FORMAT! *0 VALUES OF ENOONX, XHAT, X» AND H ARE NOT COMPATIBLE. 

.. RUN AOORTEO.') 

CALL EXIT 

end 


CRAN2580 



ANALYTICAL CELL AT THE FOLLOWING CONDITIONS 





INLET AIR PCU2 - 

4 ,00 

MM HG 

CURRENT DENSITY = 

3C. CO 

A /FT 2 

MATRIX THICKNESS^ 

25. 00 

MILS 

AIR FLEW PATE = 

i r .oc 

FT 3/M IN 

AIR TENPERATUPF = 

7€. CO 

DEG. F 

AIR FEW POINT •= 

4C, OC 

DEG. F 

CTHER CCNDITICNS 
CELL VOLTAGE 

C . 30 

VGLT 

HYDRGGEN RATE - 

C .07 

FT3/MIN 

anodf PCOPIN = 

c , 0 

MM HG 

: SR S: Ti« >S:A * iS S;* A* ***>!'*** 


OUTLET AIR D.P. = 

42.61 

OFG. F 





NC. OF SUP CELL 

TA IP IN 

TA IRQUT 

tf:lfc 

PH20I N 

PH20CUT 

CSCATH 

1 

70 .PC 

72.80 

73 .34 

6.3^ 

6.54 

5.57 

£T* 

72. ec 

75.59 

76,7.3 

6.54 

6.77 

5.4 7 

3 

75.59 

78.39 

79 ,90 

6.77 

7.01 

5.29 


UPPEF FLUX LIMIT=: .f-PBeSC 
LOWEP flux LIiVlT= ,374400 


THE INITIAL FLUX AT THF CATHODE IS .4<?gf^T'' 



ORIGINAU PAGE IS 
OF POOR QUALZn^ 


X 

. SOC-'^ 00 


CIC = 

. 10662 lD-04 

C3C = 

.335602 

C4C = .773765 

X = 

*50Cro'C 


CIA = 

. 491241D-01 

C3A = 

.1276760-06 

C4A = .491244 

THF 

IHTEGFATI CN 

ROUTINE 

TOOK 267 + 1 

STEPS 



GUF5SFO CATHODE 

FLUX - 

CALC. ANODE FLUX IS 

-335.884 

DIFF ./FLUX 

THE 

INITIAL 

FLUX 

AT THE 

CATHODE IS 

.599040 



X 

. 5CCOOO 


C 1C = 

. 136394D~’^4 

C3C = 

,317050 

C4C = ,77376? 

THF 

INITIAL 

FLUX 

AT THE 

CATHODE IS 

.54912<^ 



X =- 

•.5CC00O 


CIC = 

. 16377CD-F4 

CSC = 

.32691 2 

C4C = .773764 

THE 

INITIAL 

FLUX 

AT THE 

cathode IS 

.5491 2C 

. ^ 


X =“ 

, 5CC''CC 


CIC = 

. 16377«^0-F4 

C3C = 

.326912 

C4C = ,773764 





ORIGIiSrAi; PAGE IS 
OP POOR QUALITY- 


THE INITIAL FLUX AT THE CATHCOE IS .5241 60 

X =-,5C0CT0 CIC = .1794i3D-C4 C3C = .331390 C4C ■= .773764 

THE INITIAL FLUX AT THE CATHODE IS .51168C 

X =!-.5CCF.^C CIC = . 187807D-04 C3C •= .333528 C4C = .773764 

X = .SCOGOO CIA = .205002 C3A = ,4&7484D-05 C4A = 

THE INTEGRATICN ROUTINE TOOK 235 + 1 STEPS 

GUESSEr CATHODE FLUX - CALC. ANCOF FLUX IS -159.230 DIFF./FLUX = 311.190 

THE INITIAL FLUX AT THF CATHODE IS .517929 

X =-.5C^0G3 CIC = .18355E5D 

X = .6COOOC CIA = .1262440 

THE INTEGRATICN FOUTINF TOOK 72 


-94 C3C = .332467 C4C = .773764 

-04 C3A = ,136551 C4A = .338095 

+ 1 STEPS 



GUFSSFD CATHODE 

FLUX - CALC, ANODE 

FLU X IS 

.993974 

oiff.xflux = 

THE 

INITIAL FLUX 

AT THE CATHODE IS 

.51480C 

- — — 


X =- 

. ECODfr- 

CIC = .18567CD-04 

C3C = 

.333900 

C4C = *773764 

X = 

* 59 C9 

CIA = .246469 

C3A = 

.3491940-'' 1 

C4A = .316282 

THE 

INTEGPATI CN 

HOUTINE TOCK 130 + 

1 STEPS 



CUES 

SfcD CATHODE 

FLUX - CALC. ANCDF 

FLUX IS 

.956C83 

PI FF. /FLUX = 

THE 

INITIAL FLUX 

AT THE CATHODE IS 

.513249 

. , ,,, 


X =- 

.5C000<^ 

CIC = .1867350-04 

C3C = 

.333264 

C4C = .773765 

X t= 

.5C9099 

CIA = .239719 

C3A = 

. IC7936D-94 

C4A = *239732 

THE 

INTEGHATI CN 

FOOT I ME TOOK 225 

I STEPS 




1 . <?1 743 


1 . 85719 


GU&.SSFC- CATHODE FLUX 


CALC* ANODF FLUX IS 


- 93*891 1 


*^!FF */FLUX 


1B2.93S 



THE INITIAL FLUX AT THE CATHOCE IS .514020 

X =-.50COOC CIC = .lP62n?D-0A C3C = .333132 C4C = .773764 

X = .SCT-'-nc CIA = .262050 C3A = . 25 AS 06 D— 04 C4A = .2621C1 

THE INTEGFATICN FOUTINE TGCK 220 + 1 STEFS 

GUESSEE- CATHODE FLUX - CALC. ANODE FLUX IS -47.50 05 DIFF./f=LUX = 03.''' log 

THE INITIAL FLUX AT THE CATHODE IS .51441C 

X =-,5OOD0O CIC = .1S5936D-04 C3C = .333''66 

X = .5000CO CIA = .275268 C3A = .6572020— ''A 

THE INTEGRATION ROUTINE TOOK 190 F I STEPS 

GUESSED CATHCDH FLUX - CALC. ANGOF FLUX IS -19.5612 

THE INITIAL FLUX AT THE CATHOCE IS ,514605 


C4C = .773764 
C4A = .275419 

01 FF. /FLUX = 36.':265 


u 



X =- 

.SC GOOD 

C 1C = 

. 1 85e?3D-D4 

C3C = 

.333*^33 

C4C 


.773764 

X = 

. SCCCCC 

CIA = 

. 282327 

C3A -= 

.2914710-03 

C4A 


.282910 

THt 

INTEGRATI CH 

ROUTINE 

TOOK 161 T 

1 STEPS 





GUfSSED CATHODE 

J^LUX - 

CALC. ANODE 

FLUX iS 

-3.88692 

DIFF 

./FLUX = 

THE 

INITIAL ELUX 

AT THE: 

CATHODE IS 

,514702 

. — — II ■ 




X =- 

.5C0CCC 

CIC = 

. 1 857360-04 

C3C = 

.333716 

C4C 


.773764 

X = 

.soccoc 

CIA = 

.27296^5 

C3A = 

. 1 0286CD-0 1 

C4A 

= 

.293537 

THE 

INTEGRATION 

ROUTINE 

TOOK 134 + 

1 STEPS 





GUESSED CATHODE 

FLUX - 

CALC. ANOCF 

FLUX IS 

. 357809 

DIFF 

./FLUX = 

THE. 

INITIAL FLUX 

AT the 

CATHODE IS 

.514654 





X =- 

^ cr ^ fv f* 

CIC - 

. 18S77f'D-04 

C3C = 

.333'''24 

C4C 

= 

.773764 

X = 

. S'? C 0 r ^ 

CIA = 

. 282932 

C3A = 

. l?4643r-‘^2 

C4A 


. 28E425 





w, ,K 


THE 

IHTEGF<ATI CH 

ROUTINE 

TOOK 142 + 1 

STEPS 




GUESSED CATHODE 

FLUX - 

CALC, ANODE FLUX IS 

-.1 57754 

oiff 

./FLUX = 

THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

.51 4678 




X =- 

.5CCOOO 

CIC = 

. 1857E.3D--e4 

C3C = 

.3 33^^2 0 

C4C = 

.773764 

X = 

*5CCOOO 

CIA = 

. 27HS8C 

C3A = 

.4699850-62 

C4A = 

. 288680 

THE 

INTEGRATION 

ROUTINE 

TOOK 134 + 1 

STEPS 




CUES 

SED CATHODE 

FLUX - 

CALC. ANODE FLUX IS 

.691 C97 

OIFF 

./FLUX = 

THE 

INITIAL FLUX 

AT THE 

CAT HO HE IS 

• 5 1 4666 

. — 



X =- 

,5CC€C':'- 

CIC = 

. 185761 D-F4 

C3C = 

.333622 

C4C = 

,773764 

X - 

.sccoro 

CIA = 

. 281 7C6 

C3A - 

.2616830-02 

C4A = 

,266740 

THE 

INTEGRATI CN 

ROUTINE 

TGCK 137+1 

STEPS 





. 3fl65?4 


1 . 342T7 



GUESSED CATHODE FLUX - CALC. ANODE FLUX TS .A249D.2 OIFF./FLUX = . B25S88 

THF INITIAL FLUX AT THE CATHODE IS .51466^ 

X =-. 50 nonr CIC = .lb576eD-C4 C3C = .333003 

X = .ECOCOC CIA = .282556 C3A = .168239D-02 

THE INTEGPATICN RCUTINE TOOK 142 + 1 STEPS 

GUESSED CATHODE FLUX - CALC. ANODE FLUX IS ,141105 

THE INITIAL FLUX AT THE CATHODE IS .514657 

X =-.500000 CIC = .18576BD-04 C3C = .333324 C4C = .773764 

X = .500^00 CIA = .2o27<52 C3A = .14203FD-C2 C4A = .285633 

THE INTEGFATIGN FDUTINF TOOK 140 + 1 STEPS 

X Cl C2 C3 C4 


C4C = .773764 
C4A = .285921 

DIFF./FLUX = .274171 



-.sfrc-oo 

. 73794eE-94 

.775425 

2 .39777 

5.57104 

^’f^VslARNI NG** PREC IPX TATI CN EXPECTED AT CATHODE 



- .299999 

. 178827F-03 

.447020 

1 ,8481 1 

4 .14341 

- .999989E-0 1 

,2559C7F-C3 

.278618 

1 ,64838 

3.57563 

.icecci 

. S68852F-03 

. 1 16606 

1 ,53351 

3.18419 

.srocc 1 

♦5151 43 

. 194352E-03 

1 .24278 

3,OC08l 

.SfGCCO 

2.G361 0 

. 1 1 3971F-05 

. 102268E-01 

2 ., 05655 

FLUX ECR SUBCELL 

NUMBER 1 IS 

.514657 



UPPER FLUX L1MIT= 
LOWER FLUX LIWIT= 

.514669 
: .514654 




THE INITIAL FLUX 

AT THE cathode 

IS .504364 

_ 




'^IGINAi; PAGE IS 
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X =-^5coooe 

CIC = 

«21C294D-04 

C3C = 

.327S79 

C4C = ,760133 


X = ^5C-CC0C 

C 1 A = 

, 2C197C 

C3A = 

.40211 5D-05 

C4A = .20 1076 


THE IHTEGPATION 

FDUTINE 

TOOK 205 + 

1 STEPS 




GUESSED CATHGDF 

FLUX - 

CALC. ANODE 

FLUX IS 

-180.414 

DIFF./FLUX = 

357.706 

THE INITIAL FLUX 

AT THE 

CATHODE IS 

,509512 

, ... 



X =-.50^000 

CIC = 

.206246D-C4 

C3C = 

.326998 

C4C = .768133 


X = .50C00C 

CIA = 

.197939 

C3A = 

.8686660-0 1 

C4A = .371716 


THF INTEGftATICN 

FGUTINE 

TOOK 129 + 

1 STEPS 




guessed CATHCDE 

FLUX - 

CALC. ANODE 

FLUX IS 

1 .39459 

DIFF./FLUX = 

2.73710 

THE INITIAL FLUX 

A"^ THE 

CATHODE IS 

.506938 




X Sroorr 

CIC = 

. 2C8259D-04 

C3C = 

.327448 

C4C = .760133 




X = 

.SCOT CC 

CIA = 

.2411 59 

C3A = 

.9903170-05 

C4A - .24117E 


THt 

INTEGRATION 

ROUTINE 

TOOK 222 + 

1 STEPS 




GUESSeC- CATHOPF 

FLUX - 

CALC. ANODE 

FLUX IS 

-1P4.796 

DIFF.XFLUX = 

206. 724 

THE 

INITIAL FLUX 

AT the 

CATHODE IS 

.508225 




X =- 

.5CrrCC 

Cl C = 

. 2G725CD-04 

C3C = 

.327220 

G4C = ,760133 


X = 

. 5C00 3C 

Cl A •= 

.269390 

C3A •= 

.2759960-04 

C4A = .269445 


THE 

INTEGPATICN 

ROUTINE 

TOOK 264 + 

1 STEPS 




CUES 

SEE CATHODE 

FLUX - 

CALC. ANODE 

FLUX IS 

-46.1474 

DIFF./FLUX = 

90. SO 12 

THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

.SC8868 




X =- 

,5crnc(? 

CIC = 

. 2067470-04 

C3C = 

,327109 

r4C. = .760133 



X = ,5CCCC«"' Cl ft = .2^7258 C3A = .163543D-r 3 C4A = ,2875SS 



■■ 12IEKAI] PAGE IS 
JP POOR QUALITY 


THE INTEGR/STTCN ROUTINE TOOK 173 + 1 STEPS 


GUESSED CATHODE FLUX - CALC. ANODE FLUX IS ~1 *7Z?39 DIFF./FLUX = 15.17E6 

THE initial flux AT THE CATHODE IS ,5091 

X =-,5C0D0G CIC = .2C6A96D-04 C3C 

X = ,5C0rC0 CIA = .249855 C3A 

THE INTEGRATION ROUTINE TOOK 126 + 1 ST££=S 

guessed CATHOVDE flux - CALC. ANODE FLUX IS 

THE INITIAL FLUX AT THF CATHODE IS .509029 

X --.S^COCD CIC - . ?f)6622D-<'4 C3C = .32703 1 C4C = .750 133 

X = .500000 CIA = .284352 C3A = .S76662O-0 2 C4A = .296385 

THF INTFGPATICN ROUTINE TOOK 129 + 1 STEPS 

GUFSSHL CATHODF FLUX - CALC. ANODE FLUX IS 1.1474C DIFF./FLUX = 2.25409 


= .327^54 C4C = .760132 

= -3ai662D-01 C4A = .326229 

1.37257 DIFF./FLUX = 2.69560 



THE- INITIAL FLUX AT THE CATHODE IS .5CSD4’9 

X =-.5CC‘CTC4 CIC = .2066E4D-C4 C3C = .327^95 C4C = .769133 

X - .5CCCCC CIA = .289459 C3A ■= .391P25D-93 C4A = .295234 

THE INTEGRATION FOUTINE TOOK 1 49 1 STEPS 

GUESSED CATHODE FLUX - CALC, ANODE FLUX IS -2.46487 DIFF./FLUX = 

THE INITIAL FLUX AT THE CATHODF IS ,598989 

X =-.SCCC9€ CIC = .2966530-04 C3C 

X = .5C99C9 CIA = .289761 C3A 

THE INTEGRATION FOUTINF TOOK 146 + 1 STEPS 

GUESSED CATHODE FLUX - CALC. ANODF FLUX IS 
THF INITIAL FLUX AT THE CATHODF IS .598969 


= .32-7C38 

= .1 l2358D-<?'2 

. 5'4884E-9 


C4C = . 76^133 
C4A •= ,292998 

DIFF./FLUX = 


4, 843*"7 


, 991 936F 



X =- 

•50^0CV 

CIC - 

, ?/'6669D-C4 

C3C = 

. 327102 


CAC = .760133 


X = 

•SfCBCO 

CIA = 

.289823 

C3A = 

.593342^-'' 

3 

C4A = .29101C 


the: 

INTFGRATI ON 

ROUT INC 

TOOK 144 + 

1 STfoS 





GUF.SSFD C4THC1DE 

FLUX - 

CALC. ANODE 

ELUX IS 

-1 . 15821 


DTFF ./FLUX = 

2.27561 

THF 

INITIAL FUJX 

AT THE 

CATHODE IS 

.bCS979 





X =- 

» 5C coc-n 

cir = 

.2^'6661D-04 

C3C = 

. 32799C 


C4C = .76C133 


X = 

.sceorr 

CIA = 

. 2^39888 

C3A = 

.797 14 3D-"' 

3 

C4A = .2Q1483 


THF 

INTEGRATICN 

ROUTINE 

■^COK 143 + 

1 STEPS 





CEJFS 

SET CATHOCT 

FLUX - 

CALC. AMCOF 

FLUX IS 

-.594570 


DIFF./FLUX = 

. 991338 

THE' 

INITIAL FLUX 

AT THE 

CATHODE IS 

.588984 

— 





X =-.«Crr'CO CIC = , 2C66570-34 C3C = C4C = .7B'^133 



X - .Bf'OCvC riA = ■*2e9P7C C'3A = .894847D-03 C4A = .291660 

THF INTEGRATION F^CliTINE TOOK lAC + 1 STFPF 

GUFSFEC CATHDDR FLUX - CALF. ANODE FLUX IS -.P96205 DTFF./FLUX = .581^53 

THF INITIAL FLUX AT THF CATHQDF IS ,5ra986 

X =-.5CC00C CIC = .?''6e56D-‘34 C3C = .327089 C4C = .760133 

X = .5CCCC0 CIA = .289823 C3A = .1^1488D-<'2 C4A = .291853 

THE INTEGRATION ROUT INF TOOK 135 + 1 STpoS 

GUESSED CATHODE FlAJX - CALC. ANCCE FLUX IS -.94B8PSt-01 DIFF./FLUX = .186427 

THF INITIAL FLUX AT THF CATHOCP" IS .508988 

X 5F9'''rC' CIC = .2C6664D-04 CSC = .3270B8 C4C = .760 133 

X = .fiCCCOC CIA = .289808 C3A = .104231D-02 C4A = .291 693 

THF INTFGPATION RDUTINE TODK 134+1 STFPS 




GUESSFC CATHCDE FLUX - CALC. ANODE FLUX IS - . 5536r2F-C 1 DIFF./FLUX = .1CB765 
THE INITIAL FLUX AT THE CATHODE IS .SOBgBB 

X =-.5C0''TC CIC = .206653D-C4 C3C = .3270^8 C4C = .76C133 

X = .ECCDOC CIA = .289785 C3A = .1088420— ?2 C4A = .,2919'58 

THE INTEGRATICN i^QUTINE TOOK 138 + 1 STEPS 

X Cl C2 C3 C4 

-.SCCCOC .82S913F-04 .762733 2.355C4 5,47289 

atm NG=J'4 PPE CIPITATICN E XPECTED AT CATHODE 
-.299999 . 1982EFF-03 .4483-8 1.84439 4.13728 


99Q989E-01 


28F166F-C3 


2772 16 


1 .F5771 


3.59293 



SI mwi. 


3. 1$26Q 
? . ^7541 
2.172 1C 

FLUX FCF SUBCFLL NUMPEF 2 IS .578 9 2S 

UPPEP FLUX LI-VIT-= .50BCI8Q 
LGlvFF FLUX 1.1^17= .5D£'^‘ee 


. 1CC701 


.9726778-03 .769111E-71 1.55240 


.370001 


. 639963 


.916931E-04 1.21768 


. soor 00 


2.CS645 


.947745F-06 .7e22?CE-C2 


THE INITIAL FLUX AT THF CATHODE IS .493253 


X =- 

. SCC^Oo 

CIC = 

.237996D-04 

C3C = 

.317703 

C4C = 

. 734^42 

X = 

. 6 C 0 0 C & 

CIA = 

, 1 8? ''98 

C3A = 

,3704450-05 

C4A = 

.. 1 821 •‘4 

THE 

INTEGPA'^I CN 

routine 

TOOK 244 + 

1 STEPS 




CUES SEX' CATHODE 

FLUX - 

CALC. ANGDE 

FLUX IS 

-222.963 

DIFF 

./^LUX = 452.024 



THE INITIAL FLUX AT THE CATHOCE IS 


501121 


X =■- 

.sr 00 00 

CIC = 

. 23C8 21D-04 

C3C = 

.316463 

C4C = .734042 


X = 

, SCO coo 

CIA = 

.263814 

C3A = 

. 1 85T620-C4 

C4A = .263851 


THE 

INTEGEATIL'N 

ELUTINE 

TOOK 211 + 1 

STEPS 




GUES 

SEP CATHGPE 

FLUX - 

CALC. ANODE FLUX IS 

—66 . 6S44 

DIFF./FLUX = 

133. 4RC. 

THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

,5^5055 




X =- 

.see roc 

CIC = 

.227317D-C4 

C3C = 

.3157<)3 

C4C = .734042 


X =: 

•srocce 

CIA = 

. 166«?77D-n3 

C3A = 

. 1 VH66C 

C4A = .365376 


THE 

INTEGRATI GN 

ROUTINE 

TOOK 65 F 1 

STEPS 




GUES 

SFD CATHODE 

FLUX - 

CALC. ANCCE FLUX IS 

1 .78025 

D IFF. /FLUX = 

3.52467 

THE 

INITIAL FLUX 

AT THE 

CATHOCF IS 

.SCS-RS 

— ■ 




X =-.5C<'C?^n Cic = ^2290&2D~04^ C3C = «31612‘3 

X = .SCOCCC CIA - .24C427 C3A ■= .514951D-C1 

THE INTFGPATICN PGUTINE TGOK 122+1 STEPS 

GUE SSED CATHODE FLUX - CALC. ANODE FLUX IS 1 . 756-'^7 
THE INITIAL FLUX AT THE CATHODE IS .502195 


C4C = .734C42 
C4A = .28295? 

DIFF./ELUX == 56.-^PE7 


X =~.5CCCD0 CIC = .2?’9e€CD-C4 C3C = .316213 C4C = ,734'-<42 

X = .5Crr>"^E CIA - .292416 C3A = .219 134D-D3 C4A = .292P54 



C4C = .73494? 

C4A = .343419 

OIFF./FLUX = 3.49958 



THH AT T ON mUTiNp TOOK 16P + 1 STFPS 


GUESS FP CA~rHODE 

FLUX - 

CALC. ANODE 

FLUX IS 

-5.37424 

DTP1=./FLUX = 1C. 593'' 

THF INITIAL FLUX 

AT THE 

CATHODE IS 

.5G2S42 

. i__ 1 


X =-• 

CIC = 

. 22 92e iO-A4 

C2C = 

.316171 

C4C = .73404? 

X = .5CCC0C 

CIA = 

. 277S02 

C3A - 

.1582840-0 1 

C4A = .309549 

THE IXTEGHATICN 

FDUTINE 

TOOK 128 + 

1 STFPS 



GUESSED CATHDDP' 

FL IJX - 

CALC. ANCDE 

FLUX IS 

1 . 68677 

DIFF./FIUX = 3.35448 

THF INITIAL FLUX 

AT THF 

cathode is 

.502710 



X soonFo 

CIC = 

. 22939 lD-04 

C3C - 

.315102 

C4C = .734D42 

X ■= ..5'^- FA 00 

CIA = 

. 2940 1 5 

C3A = 

. 1 16^5^D-<"'£ 

C4A = .29633F 


THF INTFGPATItK ROUTINE TOOK 131 + 1 STFPS 



GUESSED CATHDOF 

FLUX - 

CALC. ANOOF 

FLUX IS 

.410857 

OIFF./FLUX = .817269 

THE 

INITIAL r-LUX 

AT THE 

CATHODE IS 

.502658 

. __ 


X =- 

. 5C '”500 

cic = 

.229446D-04 

C30 = 

.3162''3 

CAC = .734^42 

X = 

« 5 C C F c 0 

CIA - 

. 293633 

C3A = 

.3BC906D-F3 

C4A = .2Q4395 

THE 

INTFGFATT CN 

FDUTINE 

TOOK 152 + 

1 STEFS 



CUES 

SEC CATHODE 

FLUX - 

CALC. ANODE 

FLUX IS 

-2 . 37277 

OIFF./FLUX = 4.72''44 

THE 

INITIAL FLUX 

AT THF 

CATHODE IS 

.502689 



X =- 

.« f^C'^Tr 

CIC = 

. 2294180-04 

03 C = 

.316197 

r4C = .734042 

X = 

, FDD/) CD 

CIA = 

.2U4C72 

C3A = 

.59 1956P-03 

C4A = .295256 

THE 

INTE-GPATICE4 

FDUTINE 

TOOK 142 + 

1 STEFS 




GUt$?'^r' CATHODF FLUX 


CALC. FLUX TF 


- . -^ 9 1 5FFi 


ni F<- ./FLUX 


1 .79353 



ORIGINAL PAGE IS 
Og J>OOR QUALITY 


THF INITIAL FLUX AT THE CATHOPF IS .502704 


X =-.5CC0FC 

CIC = 

. 2294 0 50-04 

C3C = 

. 31 6195 

C4C = ,734^'42 


X = .sceccc 

CIA = 

. 294152 

C3A = 

.79234 ID— S3 

C4A = .295. '7 


THE INTEGFATICN 

FQUTIME 

TOOK 134 + 

1 STEPS 




GUFSSED CATHODE 

FLUX - 

CALC. ANODE 

FLUX IS 

-.225510 

diff./flux = 

, 44P!=i94 

THR INITIAL FLUX 

AT THE 

CATHODE IS 

.502712 

* — — _ 



X =-.5eoooe 

CIC = 

. 2293980-04 

C3C = 

.316193 

C4C = .734042 


X = .SCfFon 

CIA = 

.2941 17 

C3A = 

.960 66 10-03 

C4A = ,296038 


THE INTFGFATICN 

ROUT INF. 

TOOK 145 + 

1 STEPS 




GUESSED CATHL-r- 

FLUX - 

CALC. ANODE 

FLUX IS 

,125551 

DIFF./FLUX = 

.249808 

THE initial flux 

AT THE 

CATHODE IS 

,502708 

— 





X =-.5CCCr.C CIC = .2P9401D-Q4 C3C = .316194 C4C = .734042 

X = .SCCCr^ CIA = .2941A4 C3A = .8728030-03 C4A = ,295890 

THE INTEGRATICN ROUTINE TOOK 135 + 1 STEPS 

GUESSEO CATHODE FLUX - CALC. ANODE FLUX IS -.407192E-01 DIFF./FLUX = .B099O7E-C1 
THE INITIAL FLUX AT THE CATHODE IS .50271C 

X =-.500000 CIC = ,2293990-04 C3C = ,316194 C4C = ,734042 

X = .500C0C CIA = ,294133 C3A = .91740CD-C3 C4A = .295968 

THE INTEGRATION FCUTINF TOOK 129+1 STEPS 

GUESSED CATHODE FuUX - CALC. ANODE FLUX IS ,475653E-01 riFF./FLUX = . 9461780-*^ 1 

THE INITIAL FLUX AT THE CATHODE IS .5C27C9 

X =-.5f,OOOC CIC = .229400D-04 C3C = .316194 C4C = .734042 

X = .50COOC CIA = ,294138 C3A = ,90210 lD-03 C 4A = ,295942 



THF INTEGPATICN RDUTINi: TOOK 134 1 STEPS 


Pi 

>< ^ 







X 

Cl 

C2 

C3 

1 

1 

1 

1 

1 

1 O 

1 -F 

1 

-.5C<3CCC 

. 91 1273F-P4 

.731743 

2.27660 

5. 2S5C3 

-,299999 

,2176aOE-€3 

.445692 

1 .82748 

4. 10086 

-.999989F-C1 

.319576F-T3 

.275163 

1.65640 

3 . 58828 

.ICCOCl 

e9C4061F-03 

.91 61 62 E- 01 

1 ,56C 17 

3,21285 

,3coeoi 

.722812 

. 8791 97 E- 04 

1 .19705 

3, 11709 

,5ccccn 

2, 1 177<' 

.a5413DF-16 

.64951 3E-0 2 

2.13C78 

FLUX PDF SUPCHLL 

NUMPFP 3 IS 

.502709 



€02 RFNCVAL RATF 

= .2192BeE-01 

= .964856 

LB !VCL/FT2. 
lb /FT2.9AY 

DAY 




TRA^4SFER INDEX = 2.C-364C 
C02 PEMDVAL EFFICIENCY = 




CG2 AisCDE PRESSURE = 176. 35C 


74.12 % 


FULL SCALP: CELL AT THE FOLLOWING CCNDITIGNS : 




INLFT AIR PC02 = 

1 ,CG 

MM HG 

CURRENT DENSITY = 

15.00 

A /FT2 

MATRIX THICKnESS= 

1C. GO 

MILS 

AIR FLOW PATE = 

1C, OR 

FT3/MI N 

AIR TFRPERATURE •= 

7C.0C 

CFG. F 

AIR DEW POINT = 

5C. OC 

DEG. F 


OTHER CONDITIChiS : 

CELL VOLTAGE = C.gn VCLT 

HYDROGEN PATE = C.07 FT 3/M IN 

ANODE RCDOIN = f-.O MM HG 


OUTLET AIR D.F. = 51.15 DEG. F 

Nt . OF SUBCELL TAIPIN TATRDUT TFLFC PH20IK PH200UT CSCATH 

1 7C.CC 74*64 77.16 9.22 9.61 4.59 

UPPER FLUX LIMIT= .P496CC 
LOWER FLUX L. 1MIT= .624^090-0 1 


THE INITIAL FLUX AT THE CATHODE IS ,149760 


X =-.5fOOCr QIC = .346156D-C4 CIC = .290527 C4C = .637927 



DSIGINAD PAGE IS 
lOP POOR QUALITY 


X = .5C-0'?0C CIA = ,443485 C3A = 

THE' IK'TEGPATION ROUTINE TOOK 20€ + 1 STEPS 

GUESSED CATHDOE FLUX - CALC, ANODE FLUX IS 


THF 

INITIAL FLUX 

AT THE 

CATHODE IS 

. 199680 

X -- 

^ c:, ^ /s 

cic ■= 

. 1681 C4D-r4 

C3C ■= 

X - 

. 5C9C00 

CIA = 

.2556480-04 

C3A = 

THE 

INTEGRATION 

CC*UTINE 

TOOK 46+1 

STEPS 

GUESSED CATHODE 

FLUX - 

CALC. ANCDF, FLUX IS 

THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

.1 74720 

X =- 

,5'f*ro9A 

CIC = 

.23B699D-04 

C3C = 

X - 

.5-ro. 

CIA = 

. 347399 

C3A = 

THE: 

INTLGEAT I ON 

TOUT INF 

TOOK 116+1 

STEPS 


5254R5D-C4 


C4A = .44359F 


-67.P489 DIFF,/FLUX ~ 449. C45 


.262691 C4C = .637910 

.236172 C4A = ,534931 

.719293 Dir^E./FLUX = 3.6 9p?3 


,279325 C4C = .63792? 

.9749370-91 C4A - .542369 



GUESSED CAThDDE FLUX - CALC. ANCDH FLUX IS .6''-696C OIFF ./FLUX ~ 3. 4739S 


THE 

INITIAL FLUX 

AT THE 

CATHOCF IS 

.162240 

— 


X =- 


CIC - 

. 2882950-^4 

C3C = 

. 285424 

C4C = .637924 

X = 

. 5^' *"= ? C V 

CIA = 

. 4542C2 

C3A = 

. 1 4''072n-'' 3 

C4A = ,45448'' 

THE 

INTEGRATI CM 

fiOUTINF 

TOCK 181 + 1 

STEPS 



CUES 

ED CATHnnE 

FLUX - 

CALC. ANCDF FLUX IS 

-2 6. 

D IFF. /FLUX = 

THE 

INITflAL flux 

AT THE 

cathcde is 

. 1 684 8'' 

- 


X =- 

.SCOCTT 

CIC = 

.2625780-04 

C3C = 

.28251 7 

C4C = .e>37923 

X = 

■» S C 0 

CIA = 

. 45^364 

C3A = 

.31 2B550-C 2 

C4A = .465021 

THF 

INTFGF.'ATICN 

ROUT INF 

TOOK, 127 + 1 

STEPS 




GUFSSLL CATHODE FLUX - CALC, AHDDL FLUX IS -.6147PC OIFF./FLUX ~ 3.64»62 


171600 


THF INITIAL FLUX AT THF CATHCDE IS 


X =- 

x:; /*• .<05 r- ^ 

* T‘ K.- Jj V U ■ 

CIC = 

. 250421D-OA 

C3C = 

. 2 C 0959 

C4C = .637922 


X = 


CIA = 

. 40 11 77 

C3A = 

.5295920-0 1 

C4A = .5C7D97 


THF 

INTEGRATION 

PGUTI NF 

TOOK 119 + 

1 STEPS 




GULS 

SF.D cathode 

FLUX - 

CALC. ANCDF 

flux is 

. 5631 CO 

diff./p-lux ■= 

3.28147 

THF 

INITIAL FLUX 

AT THE 

CATHODE IS 

. 1 70''40 

— 



X 

. 50 ^•^ 00 

CIC = 

. 2564A4D-04 

C3C = 

. 281 74T 

C4C = .637923 


X 

.5C0C0C 

CIA - 

, 433483 

C3A =r 

.2491 03D-” 1 

C4A = ,4D33C'A 


THF 

INTFGF AT I ON 

FLUTINF 

took 120 -t- 

1 STEPS 




GUF S 

SFD CATHCDE 

<^LUX - 

CALC. AN DDF 

FLUX IS 

.4 758 77 

riFP./FLUX = 

2.79874 

THE 

INITIAL flux 

AT THF 

CATHODE IS 

. 1 6926 F 

— 





X ^- .SrOeCG CIC = . 2’594‘?7D-C4 C3C = -,2«21?4 C4C = .*='•'37923 

X = .5CCCOO CIA = .448B23 C3A = .H9477D-B1 C4A = .47271 B 

THF integration ROUTINE TOOK 122 + 1 STEPS 

GUESSED CATNOC-E FLUX - CALC. ANCCE FLUX IS .3=^4339 D'£FF,/FLUX = 

THE INITIAL FLUX AT TE*E CATHODE IS .16887'' 

X CIC = .261034D-04 C3C = .2B2326 

X ^ .fCOf.CC CIA = .4S5125 C3A = . A 726 1 70-*:' 2 

THE INTEGRATICN ROUTINE TOCK 121 P 1 STEPS 

GUESSED CATHODE FLUX - CALC. ANGDH FLUX IS .5C4498E-C1 

THE lEITIAL FLUX AT THE CATHODF IS ,166676 
V =:-.5fr(^r:C CIC - .26iaC5D-'‘4 C3C = .282421 C 4C = ,637923 


C4C = .637923 
C4A = .468575 

DIFF.xrLUX = 


I .798^6 


.298749 



X = .5C*acr.c 

CIA = 

. 4-E 758C 

C3A = 

.46 7244 0-*' 2 

C4A 

= .466925 

THL INTEGf-if'TICN 

ROUTINE 

TOOK 12F + 

1 STEPS 




GUESSED CATHCDD 

FLUX - 

CALC. ANODE 

FLUX IS 

-,.204022 

D IFF. /FLUX = 1.2CP56 

THE INITIAL FLU>< 

AT THE 

CATHODE IS 

.16877? 




X =--5C€CDC- 

CIC = 

. 2Gi A19D-04 

C3C = 

,282374 

C4C 

= .637923 

X = «5D0C'*S 

CIA = 

. 4SE42S 

C3A ■ 

.5639B4D-C2 

C4A 

= ,467709 

THE IHT>^GFATICN 

PDUTINE 

TOOK 120 + 

1 STEPS 



V 

GUESSED CATHt!DD 

FLUX - 

CALC, ANCDE 

FLUX IS 

-.61 “74HE-C1 

DIFF./FIUX = .361872 

THE INITIAL FLUX 

AT THE 

CATHDDF IS 

. 16S«21 

• — — — 



X =- . EC SC 

CIC = 

, 261227D-«'>4 

C3C - 


f' 4C 

= ,637923 

X = .SfSCCC 

CIA = 

,455793 

C3A - 

.617 '2BO-'’ 2 

C4A 

= .468133 

THE INI E GHAT ION 

KCUTIMF 

TOOK 120 + 

1 cjxcpc; 






" 





X 

C 1 

C2 

C3 

C4 


. 1C377FF~C3 

.527921 

2.93292 

4 . 59296 

- .299998 

. 35f 2<59F-03 

.284884 

2.^4668 

4.37859 

- .999975E-91 

. 376896^-^2 

. 26922 8F-C1 

p.oei'^s 

4. 19?8f 

, l^DC^'2 

1 ,4187c 

,5144665-24 

1 .49698 

4.41281 

.3^00^2 

2,2781 1 

,281 8975 -*'4 

> 94329'' 

4 . 16471 

.FCrCFO 

3,28171 

.346271 5-25 

. 444269f -R1 

3, 37956 


FLUX FLF SUnCFLL NUMREF 1 IS *I6RB?1 

C02 FFI^OVAL FATE = .7?7619F-'^2 LF MCL/FT2,DAY 

= .32C152 L,P /FT2.0AY 

TFAF5FEP INDEX = 1.35273 CG2 ANODF PRESSURE: = SS.5152 





OKIGINAL PAGE IS 
OF POOR gUALmj 


CD2 fFfv'OVAL FFFICIENCY 


49.1^ 



ANi^LYTICAL CELL AT THE FOLLOWING CONDITIONS : 




INLET AIR PC02 ~ 

0.25 

MM HG 

CURRENT DENSITY = 

2f. . CC 

A / FT 2 

MATRIX THICK.NESS- 

25. 00 

MILS 

AIP FLOW RATE = 

25 .CC 

ET3/NIN 

AIR TEMPERATURE ■= 

7 0.cc 

DEG. F 

AIT DEW POINT = 

6C.00 

DEG. F 


OTHER CCNDITICNS : 

CELL VCLTAGF = C , 20 VCLT 

HYDROGEN PATE = C 07 FT3/MIN 

ANODE PC021N = 0.0 MR HG 

S: *:^>k** * ft V 3j« ^ *4= A 4^3}::* 


OUTLET AIR D.P. = 60,50 CEG. F 

NO. OF SUeCFl.L TAIRIN TAIROUT tELFC 

1 7C. i^o 72.47 75.91 

UC’PEP FLUX LIRIT= .137120 

LCEJFf E?UX LIRIT= .4992000-0 1 


THE initial FLUX AT THE CATHODE IS .99B4^0r-” 


= - . 5000 CO 



X 


cic 


37B327D-CS 


C3C 


PH2CIN PH20DUT C SC ATH 

13.2>T 13.51 3.05 


435 ? 1 


CAC 


F34157 



X = .SrorCD CIA = .31P67PD-05 C3A = .69S357D-01 C4A = ,3f^6830 

TC:.: Tr.TEGPATICN POUTINE tCGK 46 + l ST£PS 

GUE:SSEL' CATHODE El UX - CALC. ANODE FLUX IS .3643B2 DIFF./FLUX = 3.C4966 

THE INITIAL FLUX AT THE CATITCDE IS .74E8CCD-r>l 

X =- ,5CT rCC CIC = .133836D-T4 C3C = .214an 

X = <5CC0Cr CIA = ,?25412 C3A = ,^30C6f>6D-'' 5 

THEl INTFGRATintM ROUTINE TOOK 252 -f- 1 STEPS 

guessed cathode flux - CALC. ANODE FLUX IS -304. 1S6 

THE INITIAL FLUX AT THE CATHODE IS .B736''Cn-^l 

X --.5COCOC CIC = .7897720-05 C3C = .lapT^i C4C = .534159 

X - .50rcf'« CIA " .308108 C3A = .221583D-DA C4A =; .308153 

THE INTrCFATICN ROUTINE TOOK 231+1 STEPS 


C4C = .534161 
C4A = .225418 

DIFF./FLUX = 4062.32 



GUESSFD CATHODE FLUX 


CALC. ANODE FLUX IS 


-76 .8655 


DIFF./FLUX 


879. 871 




A ^ 


THE 

INITIAL 

FLUX 

AT THE 

CATHODE IS 

.93600CO-C 1 


X =- 

".5C0C8H 


CIC = 

.5703350-05 

C3C = .169997 

C4C = .534158 

X = 

.SCOT 00 


CIA = 

. 231888D-04 

C3A = .172447 

C4A = .393359 

THE 

INTEGRATION 

ROUT INF 

TOOK 55 F 1 

STEPS 


GUESSFD CATHODE 

FLUX - 

CALC. ANODE FLUX IS .3416C8 

DIFF./FLUX = 

THE 

INITIAL 

FLUX 

AT THE 

CATHODE IS 

.9C4800D-01 


X =- 

.5C0C00 


CIC = 

.6762710-05 

C3C = .180261 

C4C = .534158 

X = 

. 5C00DC 


CIA = 

.213350 

C3A = .105745 

C4A = ,424842 

THE 

INTEGRATION 

ROUTINE 

TOOK 1C9 + 1 

STEPS 



3.64-966 


guessed cathode FLUX 


CALC. ANODE FLUX IS 


.322466 


DIFF./FLUX 


3, 56395 



THE INITIAL FLUX AT THE CATHODE IS .8892000-^1 


X =-.5COOOC 

CIC = 

,7323260-05 

C3C = 

.184837 

C4C = .534158 

X = .500000 

THE I NTEG RATION 

GIA = 
ROUTINE 

,329858 
TOOK 171 + 

C3A = 
1 STEPS 

.9482400-04 

C4A = .33GC47 

GUESSED CATHODE 

FLUX - 

CALC. ANODE 

FLUX IS 

-20.3480 

OIFF./FLUX : 

THE INITIAL FLUX 

AT THE 

CATHODE IS 

.897000D 

-01 


X =-,500000 

CIC = 

.7039060-05 

C3C = 

. 1 82591 

C4C = .534158 

X = .500000 
THE INTEGRATION 

CIA = 
ROUTINE 

,296242 
TOOK 111 + 

C3A = 
1 STEPS 

.3658670-0 1 

C4A = .369416 

GUESSED CATHODE 

FLUX - 

CALC. ANODE 

FLUX IS 

.284160 

DIFF ./FLUX ^ 

THE INITIAL FLUX 

AT THE 

CATHODE IS 

.8931 ODD 

-01 



228, 835 


3.16789 



X =- 

.500e00 

CIC = 

. 71T9C5D-05 

C3C = .183724 

C4C = ,534158 


X = 

.500000 

CIA = 

.335818 

C3A = ,4243640-03 

C4A = .336668 


THE 

INTEGRATI ON 

ROUTINE 

TOOK 

13T + 

1 STEPS 



GUESSED GATHODE 

FLUX - 

CALC. 

ANODE 

FLUX IS -4,45613 

DIFF./FLUX = 

49.8951 

THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

.SS505CO-0 1 



X =- 

.5C00CC 

CIC = 

.710890D-05 

C3C = .183160 

C4C = .534158 


X = 

.5CSQC0 

CIA = 

.324849 

C3A = .1109690-01 

C4A = .347043 


THE 

INTEGRATION 

ROUTINE 

TOOK 

119 + 

1 STEPS 



GUESSED CATHODE 

FLUX - 

CALC. 

ANODE 

FLUX IS .155338 

DIFF./FLUX = 

1.73552 

THE 

INITIAL FLUX 

AT THE 

CATHODF IS 

. 8S4075D-01 




X =-.5CCCC0 CIC = .714394D-05 C3C ~ .183443 C4C = .534158 



X = .5C0CC0 CIA = ^335626 C3A = ,1775660-02 

THE INTF<5KATI0N ROUTINE TOOK 126 + i STEPS 

GUESSED CATHCDE FLUX - CALC. ANODE FLUX IS -.816599 
THE INITIAL FLUX AT THE CATHODF IS .894562D-01 

X =-»5CC0CC CIC = .7126410-05 C3C = .183302 

X = .5G00C0 CIA = .331461 C3A = .5461070-02 

THE INTEGRATION FOUTINF TOOK 119+1 STEPS 

GUESSED CATHODE FLUX - CALC. ANODE FLUX IS -.359644E-01 


THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

,8948060 

-01 



X 

.50 00 CO 

CIC = 

. 71 1765D-05 

C3C - 

.183231 

C4C = 

. 53 4 158 

X = 

• 50 or 00 

CIA = 

. 328285 

C3A = 

.S16082D-02 

C4A = 

.344607 


C4C = ,534158 
C4A = ,3 423 83 

DIFF./FLUX = .4-^2034 


C4A = .539177 


DIFF./FLUX = 9,13345 


THE INTEGRATICN ROUTINE TOOK 117 + 1 STEPS 



GUESSFD CATHGOE FLUX - CALC, ANODE FLUX IS ,38655 1 E-0 1 DIFF./FLUX = 

THE INITIAL FLUX AT THE CATHODE IS .8946840-01 

X a=-.50000C eiC = .7122C3D-05 C3C = .183267 C4C = ,534158 

X = .505^’CO': CIA = .329927 C3A = .6763510-02 C4A = .343454 

THE INTFGFATICN ROUTINE TOOK 123 + 1 STEPS 

GUESSED CATHODE FLUX - CALC, ANODE FLUX IS « 3657938-01 DIFF./FLUX = 

THE INITIAL FLUX AT THE CATHODE IS ,8946230-01 

X =-.5C0CC0 ClC = .712422D-C5 C3C = .183284 C4C = ,534158 

X = .SCFOOO CIA = ,330706 C3A = .6102400-02 C4A = .342911 

THE INTEGPATICN ROUTINE TOOK 117 -J- 1 STEPS 



X Cl C2 C.3 C4 


.990775 


-4CB852 



-.sconce 

-.2 99? 99 
-.999989E-C1 

. irccoi 

.30CCC1 

.5C0C00 


. 26390 3E-C4 
.6 48013E-C'4 
■. 9 1323CF-C4 
.212346E-C3 
.9C32eCF:-01 
2.38108 


1.20660 

. 953987 

.689442 

,310105 

.7566349-03 

.4539965-05 


1 .31965 
1 . 3 I B 1 1 
1.34248 
1 .40404 
1 .4BT2S 
. 439372F-C1 


3.84592 
3.59C28 
3.37449 
3.11 840 
3 . CC559 
2.46896 


FLUX FUR SUPCELL NUMBER 1 IS .8946230-01 

C02 REMOVAL RATE = .38E582E-02 LB MGL/FT2,DAY 

- ,169656 LB /FT2.DAY 

TRANSFER INDEX = .537634 CC2 ANGDF PRESSURE = 31.GC86 

C02 REMOVAL EFFICIFNCY = 19.55 % 



FULL scale cell AT THE FCLLGViflNG CGNCITTONS i 

* « * !fc :fe * ^ ji: 3k ^ ^ ^ -sic jji ^ dt ^ ■* :* •* 

INLET AIR PC02 = 10.00 HG 

CURRENT DENSITY = 1C. 00 A /FT2 

MATRIX THICKNFSS= 25. OG MILS 

AIR FLOW RATE •= 10. CC FT3/MIN 

AIR TEKPERATUPE = 70. CC DEG. F 

AIR DEW POINT = 6C.00 DEG. F 

OTHER CONDITIONS ^ 

CELL VOLTAGE = 0.20 VOLT 

HYDROGEN RATE = C . 07 FT3/M1N 

anode PC02IN = 0.0 MM HG 


OUTLET AIR D.P. = 6C.62 DFG. F 

NO. OF SUBCELL TAIPIN TAIROUT TELEC PH20IN PH200UT CSCATH 
I 7'^.CiO 73.09 74.7a 13.29 13,57 3.74 

UPPER FLUX LIMIT- .232960 
lower flux LIMIT= .166400 


THE INITIAL FLUX AT THE CATHODE IS .199680 


X =-.500000 CIC = .9136780-93 C3C = .258306 C4C = .520307 



X - 

.50-00 00 

CIA = 

.258555 

C3A = 

.38747 10-01 

C4A = 

.366050 

THE 

INTEGRATION 

ROUTINE 

TOOK 94 F 1 

STEPS 




GUESSED CATHODE 

FLUX - 

CALC. ANODF FLUX IS 

.699134 

OIFF 

./FLUX = 3.50127 

THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

.183040 




X =- 

.500000 

CIC = 

.1002820-02 

C3C - 

.256861 

C4C = 

.520347 

X - 

. sfooor 

CIA = 

. 324463 

C3A - 

.6368220-^4 

C4A = 

.324591 

THE 

INTEGFATICN 

ROUTINE 

TOOK 155 + 1 

STEPS 




GUFS 

SFD CATHODE 

FLUX - 

CALC. ANODE FLUX IS 

-28.9639 

DIFF 

./FLUX - 158.238 

THE 

INITIAL FLUX 

AT THE 

CATHODE IS 

.1 9136P- 

. ir- — 



X =*- 

.SCO one 

CIC = 

. 9563120-03 

C3C - 

.258334 

C4C = 

,52 0326 

X = 

.500000 

CIA = 

.329863 

C3A = 

« 18945 CD-0 3 

C4A = 

.3302*2 

THE 

INTEGRATION 

ROUT INF 

TOOK 132 + ] 

STEPS 






GUESSED CATHODE FLUX 


CALCi. ANODE FLUX IS 


-9.59076 


DIFF./FLUX 


50.1 159 


THE 

INITIAL 

FLUX AT THF 

CATHODE IS 

.195520 

— 



X =- 

. 5C0000 

CIC = 

. 934542D-0^ 

5 C3C = 

.258821 

C4C = 

.520316 

X = 
THE 

.600000 CIA = 

INTEGPATICN ROUTINE 

. 328437 
TOOK 97 + 

C3A = 
1 STEPS 

*3836800-02 

C4A = 

.336114 

GUFS 

SFD CATHODE FLUX - 

CALC. ANODE 

FLUX IS 

.218477 

DIFF 

./FLUX = 

THF 

INITI AL 

FLUX AT THE 

CATHODE IS 

.193440 




X =- 

.5CCOOC 

CIC = 

* 945310D-0; 

3 C3C = 

.258828 

C4C = 

.520321 

X = 

.500000 

CIA = 

.331277 

C3A = 

.396723D-03 

C4A = 

.332C7C 


the INTEGPATIGN 9QUTTNE TCCK 110 + 1 STEPS 


1.1 1742 


GUESSFC CATHODE FLUX 


CALC. ANCDE FLUX IS 


-4.24515 


DIF^./FLUX 


21 .9456 



THE INITIAL FLUX AT THE CATHODE IS .154460 

X =-.50norO CIC = . 539897 D -03 C3C = .258824 C4C = .520315 

X = .5COOOO CIA = .331563 C3A = .8524480-03 C4A = .333268 

THE INTEGRATICK ROUTINE TOOK 102 + 1 STEPS 

GUESSED CATHODE FLUX - CALC. ANODE FLUX IS —1.55388 DIFF./FLUX = 8.19559 

THE INITIAL FLUX AT THE CATHODE IS .19SOCO 

X =-.5CCOOO CIC = .9372130-03 C3C = .258822 C4C = ,52^318 

X = .SOODCO CIA = .330912 C3A = .1669C7D-02 C4A = .334250 

THE INTEGFATICN POUTINE TOOK 99 T 1 STEPS 

GUESSED CATHODE FLUX - CALC. ANODE FLUX IS -.455919 DIFF./FLUX = 2.33845 

THE INITIAL FLUX AT THE CATHODE IS 


1 95260 



X =- . 5C0CCC 

CIC = 

. 935876D-03 

CSC = 

.258821 

C4C = .520317 


X = . 5C0A0C 

CIA = 

. 329958 

C3A = 

, 254 158D-02 

C4A = .335041 


THE INTEGPATlCiNi 

ROUTINE 

TOOK 97 + 1 

STEPS 




guessed cathode 

FLUX - 

CALC. ANODE FLUX IS 

-.465968E-’*' 1 

D IFF. /FLUX = 

. 239152 

THE INITIAL FLUX 

AT THE 

CATHODE IS 

.1953 9 9 




X =-.5COOOO f 

CIC - 

. 93529 8D-03 

C3C = 

.25BB21 

C4C = .520317 


X =: .Broooo 

CIA = 

. 329264 

C3A = 

.3139470-02 

C4A = .335543 


THE INTEGRATION 

ROUT INF 

TOOK 97+1 

STEFS 




GUESSED CATHODE 

FLUX - 

CALC, ANCDE FLUX IS 

.1*^2662 

DIFF./FLUX = 

.525423 

THe INITIAL FLUX 

AT THE 

CATHODE IS 

.195325 




X =-.5cr^oo 

CIC = 

.935541 D-03 

C3C = 

.253621 

C4C = .E?''3 17 




X = .5cecr;c CIA = -.329631 C3A = .2825540-02 C4A = .335282 

THE INTEGRATICN ROUTINE TOCK 97+1 STEPS 

GUESSED CATHODE FLUX - CALC. ANGDE FLUX IS . 3? 190 1 1 OIFF./FLUX = 

THE INITIAL FLUX AT THE CATHODE IS .195292 

X ^-.ECCQCC CIC = .9357C8D-03 C3C = .256621 C4C = .520317 

X - .5C0CCC CIA = .329798 C3A = .2630970-02 C4A = .335159 

THE INTEGRATION ROUTINE TOOK 97 + 1 STEPS 



X C 1 C2 C3 C4 


-.5CGC00 .3717C2F-02 ,125161E-01 1.86351 

-.299999 .475508 .1542556-03 1.62076 


3.74326 

3.71718 


. 1 64337 


_,ocC98«E-0 1 


1 . 09363 


534455E-04 


1 .29152 


3.67672 



. 1C CD 01 


1 . 48795 


3GC804E-C4 


. 93899 3 


3.46597 


.300C''l 1,82531 ♦157657E-04 .635374 

.500CCC 2.37454 .192994F-DS . 193?3«'5E-D1 

FLUX FOR SUeCELL NU!«PER 1 IS .19.;292 

C02 REMOVAL RATF = .84171CE-C2 LP MGL/FT2.DAY 

= .37C352 LP /FT2.0AY 

TRANSFER INDEX = 2.34726 C02 ANODE PRESSURE = 67.6904 

C02 REMOVAL EFFICIENCY = 85.36 % 


3 , G970S 
2.41315 




